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MODEL REDUCTION FOR CONTROL SYSTEM DESIGN 


Dale Enns, Ph.D. 

Stanford University, 1684 

In a number of applications areas, the control engineer is faced with 
controlling a physical system for which an analytical model can be derived in the 
form of a very large number of coupled, first order, linear, time invariant 
differential equations. This high order analytical model is an input to the 
controller design process for any design technique. 

This thesis develops an approach and a technique for effectively obtaining 
reduced order mathematical models of a given large order model for the purposes 
of synthesis, analysis and implementation of control systems. 

This approach involves the use of an error criterion which is the H-infinity 
norm of a frequency weighted error between the full and reduced order models. 
The weightings are chosen to take into account the purpose for which the 
reduced order model is intended. 

A previously unknown error bound in the H-infinity norm for reduced order 
models obtained from internally balanced realizations was obtained. This 
motivated further development of the balancing technique to include the 
frequency dependent weightings. This resulted in the frequency weighted 
balanced realization and a new model reduction technique. 
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Two new approaches to designing reduced order controllers were also 
developed. The first involves reducing the order of a high order controller with 
an appropriate weighting. The second involves linear-quadratic-Gaussian 
synthesis based on a reduced order model obtained with an appropriate 
weighting. 

Several numerical examples are used to illustrate the theoretical 
developments of the thesis. These examples include aircraft and large space 
structure problems. The examples clearly illustrate the usefulness of this 


research for practical problems. 
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I. INTRODUCTION 


In a number of applications areas, the control engineer is faced with 
controlling a physical system for which an analytical model can be derived in the 
form of a very large number of coupled, first order, linear, time invariant 
differential equations. This high order analytical model is an input to the 
controller design process for any design technique. 

Two important control law design techniques: linear quadratic gaussian 
(LQG) [Refs. 1-5] and H ^ optimization, [Ref. 5] result in high order control laws 
that are at least as complex (same number of equations) as the differential 
equation model of the system to be controlled. These control laws are typically 
overly complex and simpler designs are sought. 

Because of computational and other practical limitations, the order (number 
of equations) of such a model or controller must be reduced for synthesis, analysis 
and implementation of the control system. Nowhere is this need more clear than 
in the area of flexible vehicle control (aircraft or large space structures) where an 
infinite number of resonant frequencies (hence, equations) characterize the flexible 
vehicle. 

Using model reduction as a part of the controller design process is in effect a 
technique for designing a reduced order controller. The reduced order model 
must characterize the physical system with sufficient fidelity such that 
performance objectives (including, but not limited to stability) for the controlled 
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physical system can be met by designing control laws with the reduced order 
model. 

The intent of this thesis is to develop well understood tools for effectively 
obtaining reduced order mathematical models of a given large order model for the 
purposes of synthesis, analysis and implementation of control systems. 

Related Literature 

There is an enormous amount (e.g. see reference list of Ref. 6) of model 
reduction literature, however, very few of the researchers have studied this 
problem from the perspective of controller design based on the reduced order 
model. Indeed, most present their results without any mention of the controller 
and the impact of the model reduction on stability or other performance 
objectives. 

The most common approach cited is to minimize the integral squared 
impulse response error between the full and reduced order model (e.g. Refs. 7-10). 
For some unknown reason this error criterion was adopted and algorithms for 
solving the minimization problem have been developed. Unfortunately, a closed 
form solution to this problem has not been found and the iterative algorithms 
suffer from such difficulties as: choice of starting guesses, convergence, multiple 
local minima, etc. This approach is somewhat dissatisfying since there is no good 
reason for choosing such an error criterion and its lack of desirable solution 


properties. 
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Another approach to the model reduction problem is to neglect states of the 
internally balanced realization of the full order model [Ref. 11]. This technique 
does not attempt to minimize any error criterion. Rather, it is based on 
transforming the realization of the full order model into balanced coordinates 
where the states of the balanced realization are as controllable as observable in a 
well defined sense. The reduced order model is then obtained by neglecting the 
states which are weakly controllable/observable. In contrast to the previous 
approach, this technique permits a closed form solution involving standard 
matrix software. 

A closely related approach is the Hankel norm model reduction technique, 
[Refs. 12-14]. For any error criterion this is the only model reduction technique 
which results in a closed form solution for the optimal reduced order model. 
Optimal here is in the sense of the Hankel norm which will be defined in a later 
chapter. 

A characteristic common to all of these approaches is that the ultimate use 
for which a reduced order model is sought is not formally a part of the model 
reduction process. For the purposes of this thesis, the ultimate use is for design, 
analysis and implementation of feedback control systems. 

Thesis Contributions 

An approach to the model reduction problem where control design is a 
formal part of the model reduction process was developed. This approach 
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involves the use of an error criterion which is the H ^ norm of a frequency 
weighted error between the full and reduced order models. The motivation for 
choosing this error criterion (to be defined precisely in a later chapter) was the 
relation between model uncertainty and closed loop system behavior. 

The balancing technique was extended to include the frequency dependent 
weightings. The weightings are added in the context of the balancing method i.e. 
controllability/observability aspects of a weighted system. 

A previously unknown error bound (in the H <*, norm) for the internally 
balanced realization was discovered. This bound provides the missing error 
criterion for the balancing technique. This bound was compared to an error 
bound (in the H norm) for the Hankel norm technique which was obtained by 
Glover [Ref. 15). Unfortunately a corresponding, simple error bound (in the H x 
norm) for the weighted case was not obtained. 

An explicit weighting for reducing the order of a controller was developed. 
This results in a straightforward reduced order controller design algorithm, 
however, the full order controller must first be designed with the full order 
model. 

Another approach to the reduced order controller design problem, where a 
reduced order model is used to design the low order controller directly was also 
developed. This approach, however, requires an iterative solution for the reduced 


order model. 
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The theoretical approaches were applied to several example problems to 
illustrate the methodology. A significant design/demonstration example was 
carried out for the ACOSS II large space structure. This example clearly 
demonstrates the effectiveness of the developments of this thesis. This example 
also shows the usefulness of this research for application to practical problems. 

Organization of the Thesis 

The remainder of this thesis is organized into five chapters. Chapter II will 
cover the necessary background material on control system design and model 
reduction. Chapter HI will cover the new model reduction technique. Chapter 
IV will cover control system design with reduced order models. Chapter V will 
present the examples and Chapter VI the conclusions and recommendations for 
further research. 

Chapter II will serve as a reference for the work that follows. In the area of 
feedback controller design some feedback fundamentals as well as the LQG 
synthesis procedure will be discussed. In the area of model reduction the 
internally balanced realization technique will be discussed. No attempt will be 
made at completeness, rather, only those aspects that directly relate to this work 
will be addressed. 

Chapter III will discuss the new model reduction approach and techniques. 
An algorithm for computing the frequency weighted balanced realization will be 
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derived. The properties of the resulting reduced order models will also be 
discussed. 

Chapter IV will develop two new approaches to reduced order controller 
design. Both of these approaches are developed by connecting model reduction 
and controller design together with the relationship between model uncertainty 
and robust design. The first approach involves compensator order reduction and 
the second involves plant order reduction. 

Chapter V will present several design examples in detail. The control 
objectives will be stated and reduced order controllers are designed to meet these 
objectives using the techniques developed in this thesis. 
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H. BACKGROUND MATERIAL 

The purpose of this chapter is to cover background material in the areas of 
model reduction and feedback control system design. Some mathematical 
analysis notation and tools will be covered briefly. 

Notation 

The transpose of the matrix , A will be denoted by A T . The complex 
conjugate of the matrix, A will be denoted by A*. The conjugate transpose of 
the matrix, A will be denoted by A E . Eigenvalues and singular values of the 
matrix, A will be denoted by X[A] and <r[A] respectively. The maximum and 
minimum singular values of the matrix, A will be denoted by <r[A] and c[A] 
respectively. The trace and determinant of the matrix, A will be denoted by 
tr\A) and det[A] respectively. A positive definite matrix, A will be denoted by 
A > 0 and a positive semi-definite matrix, A will be denoted by A >0. Real 
and imaginary parts of the complex matrix A will be denoted by Re[A] and 
lm[A] respectively. That is A = Re[A] + j Im[Aj where j 2 = -1. An n X n 
identity matrix will be denoted by I n or I depending on context. A diagonal 
matrix, A will be denoted by A = diag {a, } where a,- is the (»',») element of the 
matrix A. The notation: block col {A,} and block row {A,} is defined 
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block col{A,} = [A ji 2 A n ] 

A 2 

block row {A,} = 

where A,- are matrices of compatible dimensions. 

The symbol, ~ is read equals by definition. The symbol, V is read for any 
or for all. The symbol, R represents the field of real scalars. The symbol, R n 
represents the space of n X 1 real vectors. The symbol, e is read is an element 
of (e.g. are R n ). The symbol, 5(f) stands for the Dirac delta function (also called 
the unit impulse function). 

Linear Systems Terminology 

Consider the linear system given by 

x = Ax + Bu areR n ueR m 

y — Cx + Du yeR p 

The order of the system is the integer, n. The system has m inputs and p 
outputs. The transfer function of the system is given by C{sl-A)~ x B + D. 
The symbols {A,B,C,D}„ or {A,B,C,D} (when D— 0, {A,B,C} n or {A, B,C}) 
represent a realization (above differential equations) of the transfer function of 
the system. Note that realizations are not unique. The statements (A ,B) 
controllable and (A,C) observable have the standard meaning, [Ref. 34]. The 
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statement (A,B) stabilizable will mean that the uncontrollable modes of the 
system are stable. The statement ( A,C ) detectable will mean that the 
unobservable modes of the system are stable. The realization {A,B ,C,D} will be 
said to be minimal if it is completely controllable and observable. The system is 
said to be asymptotically stable if Re[X[A]] < 0. The system is said to be 
minimum phase if Re[z] < 0 where 2 is a transmission zero of the system. 
Transmission zeros are defined for minimal systems to be the values of z such 
that 

loses rank 

The Laplace transform of the function, /(f) will be denoted by L[f (t)](s), 
or when context permits / (s), and t will be used exclusively to denote time and 
s will be used exclusively to denote the Laplace variable. This is an abuse of 
notation but context will always determine whether the function or the transform 
is meant. 

The function, /(/) which is the inverse Laplace transform of the function, 
F(s ) will be denoted by L -1 [F (s )](t ). That is 

/(<) = L-‘[F(s)l(t) 

F(«) = L|f(t)](s) . 

Singular values will be used throughout this thesis and to make the thesis 
self contained some definitions, identities, inequalities and theorems are tabulated 


in Appendix A. 
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Analysis Tools 

The analysis tools needed for the development of the thesis ideas are several 
norms. These norms will now be defined. 

Absolute value: The absolute value of the complex scalar, z will be denoted by 
| z | and is defined by | z | 2 A zz* . 

Vector norm: The norm of the complex vector, x will be denoted by ||ar|| and is 
defined by ||arj| 2 A x H x. 

Maximum singular value: The maximum singular value (sometimes called the 
spectral norm) of the complex matrix, A will be denoted by <?\A ] and is defined 
by 0 2 {A] A X max [A F A). 

L 2 norm: The L 2 norm of the complex matrix valued function, f (x) of a single 
real variable, x will be denoted by ||/(z )||2 and is defined by 

OO __ 

ll/(*)llf ~ / o tr[/ H (x)f(x)\dx. 

L oj,: The L ^ norm of the complex matrix valued function, f (x) of a single real 

variable, x will be denoted by ||/(a:)|| 0 o an d is defined by 

ll/OOIloo = sup a[f{x)). 

2 

H ^ norm: The H <*, norm of the complex matrix valued function, f(s) of a 
single complex variable, s which is analytic in the closed right half plane will be 
denoted by ||/(s)|| 00 and is defined by 11/(5)1(00 = sup cr[f{ju)]. 
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Hankel norm: The Hankel norm of the complex matrix valued function, f (s) of 
a single complex variable, s which is analytic in the closed right half plane will 
be denoted by ||/(s)||jy and is defined by ||/(s)j|# where / («) is 

assumed to have a state space realization: f[s) = C[sI-A)~ l B and U, Y satisfy 
the Lyapunov equations: 

AU + UA t + BB t = 0 
A T Y + YA + C T C = 0 . 

Note that this definition is equivalent to other definitions of the Hankel norm 
given in the literature. A finite dimensional realization of f(s) was assumed as a 
convenience for this thesis. 

General Framework for Control Systems Analysis and Synthesis 

The purpose of this section is to briefly review a new framework for analysis 
and synthesis of control systems developed by Doyle, [Ref. 5j. This review will 
state the assumptions which imply the use of the specific analysis and synthesis 
tools to be used in this thesis. The point is that other assumptions could just as 
well be made and they in turn would require other analysis and synthesis tools. 


A general control system is depicted in Fig. II. 1. The diagram already 
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COMMANDS 
DISTURBANCES 
SENSOR NOISE 
INITIAL CONDITIONS 



J RESPONSES 
| ERRORS 


Fig. H.l General Control System 


makes the assumption that the actual physical system can be represented as a 
nominal system plus perturbations. This assumption is necessary since it is 
desirable for the nominal model (typically linear, time invariant ordinary 
differential equations) of the system to be as simple as possible (otherwise we 
cannot do much analysis, let alone synthesis). This is also the motivation for 
model reduction. 

Typically the input signals are given for a control problem. That is their 
magnitude, energy, frequency content, ... are specified in some way. The output 
signals on the other hand are functions of the inputs and the rest of the system. 
Typically some desirable output signal magnitude, energy, frequency content, ... 
are specified. These specifications mathematically take the form of frequency 


weighted norms. 
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The perturbations denoted by A can typically be bounded in some way. In 
general, these bounds mathematically take the form of frequency weighted norms. 
A simple example is the system in Fig. II. 2 where the time constant of the 
nominal linear system 

x(t) — a:r(£) + u(£) 
y{t) = x{t) 

is uncertain but bounded above and below with the absolute value norm, i.e. 

I A I < — I a I where a == constant 

i 10 1 



Fig. II.2 Simple Uncertain System 


Another less trivial example is an actuator with a finite bandwidth and rate 


limit shown in Fig. II.3. This system can be nominally modeled as unity 




Fig. II.3 Finite Bandwidth Actuator with a Rate Limit 


plus a perturbation as shown in Fig. II.4. In this case the uncertainty is bounded 
using the absolute value and norms and is expressed [Ref. 16] 

IMOIloo <na- => 
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A(S) 


4 » 


1 — 


Fig. II.4 Actuator Nominal Model Plus a Perturbation 

Note that the bound for A(s) depends on the signal level. 

Without elaborating any further here (see Ref. 5 for more details) Doyle’s 
Table will be introduced now and is shown as Table II. 1 . This table lists out 
various assumptions and the analysis and synthesis tools that result. The 
assumptions that will be used for the remainder of this thesis are found in the 
second and third rows of Table 11.1. 

The second row assumes the input signals have bounded L 2 norms and that 
the output specifications are expressed in terms of L 2 norms. Actually the 
frequency content of these signals is also of interest but without loss of generality 
the weightings for this frequency content can be absorbed into the nominal 
model. The perturbations are assumed to be zero. These assumptions result in 
the use of Bode plots of singular values as the analysis tool and for this thesis 


loop shaping as a synthesis tool. 




Table II. 1 Assumptions and Corresponding Analysis and Synthesis Tools 
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The third row assumes zero inputs and guaranteed stability as the only 
output specification. However, a specific class of perturbations is permitted. 
These assumptions also result in Bode plots of singular values and loop shaping 
as analysis and synthesis tools respectively. 

Note that for these two rows performance and stability robustness are 
treated separately. This is, in fact, one of the weaknesses of loop shaping as a 
design tool and singular values as an analysis tool. Singular values can be used 
in general but may be too conservative. In general, the weakness of loop shaping 
is that we are typically concerned with properties of more than one of the 
feedback system’s loops simultaneously. 

This weakness is eliminated with the use of the structured singular value 
[Refs. 17, 18]. The structured singular value permits the non-conservative 
evaluation of robust performance. That is, it answers the question of whether or 
not the performance specification is met in the face of the perturbations (not 
necessarily small). This is not otherwise possible. 

Although the structured singular value and H <*, synthesis techniques [Ref. 
19] are very promising tools, they were considered beyond the scope of this thesis. 

Why Bode Plots of Singular Values? 

As already indicated in Table 33.1, using Bode plots of singular values is a 
consequence of the assumptions of the second and third rows of Table II. 1. The 
purpose of this section is to develop these consequences in more detail. 



- 18 - 


Performance 

Let (7(s) be the asymptotically stable transfer function from the input, u to 
the output, y, i.e. 

y(«) = <?(*)«(*) 

Then the following two statements are equivalent {Ref. 20]: 

IMOII2 < 1 for a11 IMOII2 < 1 

and 

Halloo < 1 . 

In other words, the H ^ norm of the transfer function (or the L <*, norm of 
the transfer function’s frequency response) measures the worst case gain of the 
system. The idea is that for an arbitrary input, «(/), with unit L 2 norm, the 
transfer function’s norm gives the worst case L 2 norm of the output, y(t). 

Typically frequency content is also important and the utility of this 
approach to measuring performance is increased by introducing the weighting 
transfer functions W’ (s) and ^(s) respectively. In this case the following two 
statements are equivalent [Ref. 20]: 

||L' i |W y (s)y(s)](t)|| 2 < 1 for all J|L - J [W u (s)u ( s )] (t)|| 2 < 1 


and 
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||Hy«)(7(.)W^IL< 1 • 

The Bode plot of the singular values of G(ju>) or W y (joj) G(ju)W~ 1 (juj) is 
just a graphical means of assessing performance (because the magnitude of the 
highest peak of the frequency response is the H ^ norm). Often the weightings 
are not explicitly constructed, rather, the singular values of G{ju) are plotted 
and the effect of certain weightings can be readily assessed by the shape of the 
frequency response. 

Stability Robustness 

The stability robustness of a MIMO control system can be determined with 
Bode plots of singular values with the aid of the following theorem. Consider the 
perturbed feedback system shown in Fig. II. 5 where G(s) is the nominal open 
loop transfer function and A(s) is the 



Fig. H.5 Perturbed MEMO Feedback System 


transfer function of the additive perturbation. 
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Stability Robustness Theorem: 

Assume that the nominal closed loop system is stable i.e. (?(«)[/ + G'(s)] -1 
has no poles in the closed right half plane. Then for any stable A($) such that 

||L( S )A(„)J?-'(«)IL< i 

if 

IlflWi/ + (?(»))-■ /,-■(. Illco < 1 

the perturbed closed loop system is stable [Ref. 18]. 

Note that this test evaluates closed loop stability of the perturbed system for 
a large class of perturbations with a single test. 

Actually the Stability Robustness Theorem stated above is a special case of 
a more general theorem in Reference 18 but it will be adequate for the purposes 
of this thesis. 

Three stronger requirements for robust stability are: 

s\I + G'(;o?)] > a[A(;u;)j \/oj 

s\I + G-'ljoj)} > 5lG _1 (.;w)A(/u>)] Vw 

z\I + G~ { (ju)] > W[A{juj) G -1 (;w)] \/uj 

i.e. satisfaction of any of these three requirements implies that the condition of 
the stability robustness theorem is satisfied. A proof of this statement is given in 
Appendix B. 


The advantage of these more conservative requirements is that bounds for 
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the perturbation by itself or the perturbation normalized by the nominal are 
most times easily obtained or estimated and then can be readily compared to 
something that depends only on the nominal system. 

The Bode plot of singular values is again just a graphical means for checking 
whether the control system has stability robustness with respect to a given set of 
perturbations. For example consider the set of perturbations 

(A(s) such that a(G“ 1 ( i ;'w)A( < ;'w)] < l(w) \/w} 

where l(ui) is a known function of frequency. Then if the Bode plot of 
<z\I + G~\ju)\ lies above the Bode plot of l(u>) the perturbed closed loop system 
is guaranteed to be stable for any perturbation in this set. 

Closed Loop Responses 

The purpose of this section is to derive expressions for a MIMO control 
system’s closed loop responses. Consider the MIMO feedback control system 
show-n in Fig. II. 6 where G(s) is the plant transfer function; K(s) is the 
compensator transfer function; y c is the commanded value for the output, y; u is 
the plant input; d is the disturbance and n is the sensor noise. 
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Fig. H.6 MIMO Feedback Control System 


Employing standard feedback algebra results in the following closed loop 
responses where for notational convenience, the Laplace variable notation, s has 
been suppressed: • 

y = d + GK[y e - (n+ y)} 

= (/ + GK)- 1 \d + GK(y e -n)) 

e = y c ~ y 

= -(/+ GK)~ l d + [I- (/+ GK)~ l GK\ y t + (/+ GK)~ l GK n 
= (/+ GK)' 1 (y e -d) + (/+ GK)' 1 GK n 
u — K\y c - (n+ d+ Gii)\ 

= (1+ KG)- 1 K(y c -n-d) . 


Desensitization 

One important goal of feedback is to accomplish desensitization. 
Desensitization means making a system or system component insensitive to 
changes in operating conditions. A classic example is the feedback amplifier 
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where some of its component’s dynamic characteristics vary widely with 
temperature, however, with feedback this variation can be made negligible. 

The purpose of this section is to quantify the effect of feedback in 
accomplishing this desensitization. Note that this will not be a stability 
robustness analysis (which requires another analysis), here, stability is assumed. 

To quantify the desensitizing effect of feedback in the multivariable case, the 
comparison sensitivity approach will be used [Refs. 21, 22). In order that we 
compare “apples to apples,” the relative sizes of perturbations to two transfer 
functions (open and closed loop) which are nominally the same will be compared. 
To facilitate this comparison, constant pre- or post-compensation will be used 
such that the nominal closed loop transfer function will be identical to the 
nominal open loop transfer function. 

Output Sensitivity 

For the effects of a perturbation on the output consider the open loop 
system in Fig. II. 7 and the closed loop system in Fig. D.8 where (?(«) is 
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Fig. II.8 Closed Loop System With Perturbation and Post-compensation 

the nominal plant transfer function, K(s) and A (s) are the transfer functions of 
the compensator and perturbation respectively. 

For notation al convenience, the Laplace variable notation, s will be 
dropped. Let the perturbed plant be denoted by 

g p J a + a 
= U+P.)-' a 
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where 

P, A - AG; 1 

is a normalized (from the right) version of A(s) (i.e. % error). 

Then the transfer functions from y to y for Figs. 13.7 and n.8 are given by 

T 0 i = G p K = open loop transfer function 

T d = (/+ GK){I+ GpK)' 1 G p K = closed loop transfer function 

It can be verified that when A = 0 (i.e. nominal condition) the two transfer 
functions: T o! and T c i are the same. What is of interest here is the effect of the 
perturbation on T ot and T c i. The effect of the perturbation on the open loop 
transfer function is given by 

T., = (/+/>.)-* GK . 

To obtain an analogous expression for the closed loop transfer function, 
substitute for G p and perform the following algebra: 

T a = (/+ GK)\I + (/+P,)- 1 GK}-' (/+FJ- 1 GK 
= (1+ GK)}(I+P,)-' (I+P.+ GK)}-' (/+/>„)-■ GK 
= (/+ GK)(I+ GK+ P,)-' GK 
= (/+ GK) {[/ + P,(I+ GK)~' } (/+ GK) }-> GK 
= \I + P,[I+ GK)-' ]-' GK 
= [/+ (P.)a I' 1 GK 


where 
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(P.)d M PAI+GKr 1 . 

At this point the effect of feedback on the perturbation is readily seen. The 
closed loop perturbation is essentially a scaled version of the open loop 
perturbation. In terms of the spectral norm we have 

WM = *\PAi+ cat 1 ] 

< V\P.) *[(/+ GKY' ] 

= z [/+ gk] ' 

That is if q\I+ GK] > 1 holds at some frequency ( s=ju> ) then the use of 
feedback accomplishes desensitization at that frequency. Bode plots of singular 
values again can be used as a graphical means of assessing a feedback system’s 
desensitization properties. 

Input Sensitivity 

For the effects of the perturbation on the input consider the open loop 
system in Fig. II.9 and the closed loop system in Fig. ELIO. 
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Fig. II.9 Open Loop System With Perturbation 



Fig. n.10 Closed Loop System With Perturbation and Pre-compensation 


Let the perturbation be normalized from the left, then 


o r = 0(1+ py 

where 

Pi =£ -G - 1 A 


Performing the same steps as for the output results in the closed loop 
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perturbation being given by 

(Pih £ ( I+KGY'P i ■ 

At this point the effect of feedback on the perturbation is readily seen. The 
closed loop perturbation is essentially a scaled version of the open loop 
perturbation. In terms of the spectral norm we have 

w,-y = ffi(/+ kg )-' p,i 

< y|(/+ KG)-' 1 5|P, j 

= a.[l+KG\ • 

That is if g\I+KG] > I holds at some frequency (s—joj) then the use of 
feedback accomplishes desensitization at that frequency. Bode plots of singular 
values again can be used as a graphical means of assessing a feedback system’s 
desensitization properties. 

Feedback System Transfer Functions 

The purpose of this section is to define some special transfer functions which 
will be used for assessing a feedback system’s performance, desensitization and 
stability robustness properties. The following terminology will now be attached 


to some of the transfer functions from above. Let 
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L 0 (s) A G(s)K(s) = output loop transfer function 
L,(s) A K{s)G(s ) = input loop transfer function 
/+ (S'(s)A'(s) = output return difference 
7 + /v(s)G'(s) = input return difference 

,?,,(«) A [I + G ? (s)7r(s)]“ 1 = output sensitivity 
^•(s) [/ 4- ^(s)G(s)]' 1 = input sensitivity 

H 0 (s) A [I + G(s)/^(s)]" 1 G(s)K(s) — output closed loop 
H { {s) A /i(s)G'(s)[7 + A(s)G(s)] -1 = input closed loop 
/ 4- [G'(s)/C(s)] _1 = H~ l (s) = output inverse-return difference 
/ -f [K (s )G'(s )] -1 = Hf\s) = input inverse-return difference 

The feedback properties are summarized in Table H.2. 


Typical Control System Design Problem 

The design problem assumes that the external inputs (commands, 
disturbances, sensor noise) are specified in some way (weighted L 2 norm bounds). 
An assumption regarding the uncertainty of the nominal model (frequency 
response error bound) is also required. It is also assumed that the objective of 
the design problem (small error between the output and commanded output) is 
specified in some way (weighted L 2 norm bounds again). Another objective may 
be that control energy is minimized in meeting the small error objective. 
Satisfying the control energy objective will be treated as a secondary objective 
and will only be mentioned briefly in a later section. 


Mathematically the control design problem takes the following form. 
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Table 13.2 Summary of Feedback System Properties 


Closed Loop Responses: 


yis) 

-i 

H.(«) S.(s) - H,(s ) 

«• 

— 

S.(>) -SJs) H.(») 

L«(«) 

• 

HMG-‘(s) -HMG-'(s) - H;(s)G-'(s ) 



Stability Robustness: 

> a[(? _1 (;'u>)A(;a;) j X/u> 

] > ff[A(jw)G- l {ju) ] \/w 

Desensitization: 

frjb’ilrfO'w)] < »[w«) ]»[«(»«)] V" 

s ^ p ,)Ai u ) ] < ] sjs.ljw) ] Vw 
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Assume that weightings have been obtained such that the magnitude and 
frequency content of the commands, disturbances and sensor noise for the control 
problem are described by 

ll L_t |W,,(s)y c (s)] (t)|| 2 < 1 
||L- 1 [W d (s)d(s)](t)|| 2 <l 
||L-'|W„( S )n(s)](t)|| 2 <l . 

Likewise assume that a weighting has been obtained such that the 
acceptable magnitude and frequency content of the error (i.e. e A y c -y) is 
described by 

l|L -1 [W e (s )e(s )] (t)|| 2 < 1 . 

Assume that the modeling process has produced a nominal model of the 
plant, <7(s) (transfer function corresponding to linear, time invariant, ordinary 
differential equations) and one or both of the input and/or output scalar, 
frequency dependent, multiplicative uncertainty bounds and/or / 0 (w). These 
bounds describe the uncertainty between the nominal model, G(s) and the true 
system Gtruci 5 ) 35 follows. Although the true system is not precisely known we 
will assume it belongs to one or both of the sets: 

{Gtru.( s ) : 5 , [C"‘(yw) |G lnle (ju) - G(juj)] ] < l,(u) } 

{G lr Js): <i{ [G true (i W ) - C(iw)] G-‘(ju>)i < } . 

Given the weightings W^s), W d (s), IT b (s) and ITJs), the nominal model, 
G(s) and its uncertainty bounds, /,-(w) and l 0 (u) the control system design 
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problem is to find the compensator transfer function such that the performance 
objective is satisfied for the feedback system of Fig. 11.11. 



Fig. n.ll Feedback Control System Block Diagram 


Although Doyle’s structured singular value and H^/L^ synthesis method 
[Refs. 5, 17, 18, 19] solve this problem, the theoretical background required is far 
beyond the scope of this work, hence, this thesis will only discuss an approximate 
solution. In addition to the substantial theoretical understanding required, the 
synthesis method requires significantly more numerical computations to 
obtain the compensator transfer function, K(s). 

A major simplifying assumption is to require only that the performance 
objective be satisfied for the nominal system (second row of Table II. 1) although 
stability is required for the true system (3rd row of Table II. 1). The other 
simplification is that no formal optimization will be attempted. 
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Performance and Uncertainty Requirements 

To meet the performance objective for each one of the command, 
disturbance and sensor noise individually requires that 

||iv,(,)|/ +/,,(.)]-» W->(.)IU <1 
iiw',M/+ i.wr 1 hw)Moo < 1 
II "'.(•P + i • 

To satisfy the output stability robustness requirement, requires that 

<l\I + > '» • 

To gain a clearer picture of these requirements, note that they are satisfied if 

< S.1 1 + Lo(ju)} 

< e\l + L 0 (joj)] 
vW 1 + LAJ")}- 1 MJu)) < 

v{\I + L 0 {ju)Y l L 0 (joj)} < / 0 "V) 

where w £ (w) £ <?[W e {ju)] 

Mu ) J c ( w *(;<*>)1 * = 


It will be helpful to choose the simple performance objective 

IWOlb J IMO - »(Olb < ‘ 

=> W'.i,) = | 

and plot — w~ l {u>), — €tr n (w) and /~ 1 (cu) versus frequency. Such a plot 

is shown in Fig. 0.12 where a log-log scale has been used. 
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Satisfying the requirements for the commands and disturbances requires that 
a plot of sl[I 4- L 0 (ju ) )] lie above the command and disturbance curves in Fig. 
11.12. Satisfying the requirements for the sensor noise and uncertainty requires 
that a plot of a[ [I + L 0 {ju)\~ l L 0 {ju) ] lie below the sensor noise and 
uncertainty curves in Fig. 11.12. 



Fig. H.12 Typical Magnitude and Frequency Content of Commands, 
Disturbances, Sensor Noise and Uncertainty 
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Large and Small Loop Gain Approximations 

Although the requirements are in terms of q\I + L 0 (ju)} and 
<r[ [I + L 0 {ju)\~ l L 0 (ju>)} they can be well approximated by requirements on 
large and small loop gains as follows. Note that if &\I + L 0 (ju)] is large (i.e. 
>> 1) then &[L 0 (ju))] is also large in fact 

z[/+ L 0 (ju)} = a\ L 0 {ju)] 

when 

dMM I » 1 

Also note that if a\\I + L 0 (j cj)] -1 L 0 {jui)} is small (i.e. << 1) then a \L 0 (ju))} is 
also small in fact 

a[[I + L 0 (ju)}~ 1 L 0 (ju)\ = o\L 0 (ju)\ 

when 

« 1 

These observations suggest that Fig. 11.12 can be used to determine 
graphically the suitability of any L 0 (s) and thus of any compensator, K(s) (since 
L 0 (s) ~ G(s)K(s)) by simply sketching the singular values of L 0 [joj). Such a 
plot is shown in Fig. 11.13. 
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Fig. 11.13 Loop Shape Constraints and Trial Design 

The secondary objective of minimizing control energy can be approximately 
satisfied by having the c[L 0 {j(X>)} line follow as closely as possible to the 


command and disturbance limits. 
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Crossover Requirements 

The large and small loop gain approximations are not valid when the 
singular values of the loop transfer function are near unity (i.e. crossover) and 
thus additional constraints are required for crossover. These constraints are that 
the minimum singular values of the return difference and the inverse return 
difference don’t get too much smaller than unity near crossover. This is such 
that performance and stability robustness are not compromised for frequencies in 
the crossover region. Recall that when the minimum singular value of the return 
difference is less than unity, disturbance and command responses of the error and 
sensitivity are amplified compared to open loop. When the minimum singular 
value O'f the inverse-return difference is much less than unity the closed loop 
system could be unstable for small perturbations. 

SISO Interpretation 

The crossover requirements for SISO systems have historically [Ref.23] been 
expressed in terms of gain and phase margins. The relationships between the 
gain and phase margins and the magnitudes of the return difference and the 
inverse-return difference will now be derived. 

Phase Margin 

By the definition of the phase margin, PM of the loop transfer function, 
L(s); L{ju ) = - (cos PM + j sin PM) for some frequency, u. In this case the 
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magnitude of the return difference and the inverse-return difference are the same 
which can be seen by 

|i+ iW«)| _ u + . 

In terms of the phase margin they are given by 

|1+ L(ju) | = |1+ L~ l (ju) | = \/(l - cos PM) 2 + sin 2 PM 

= V2 >/ 1 — cos PM 

Thus clearly a small phase margin (say j PM | < 30 deg.) implies a small return 
difference and inverse-return difference magnitude. 

Gain Margin 

By the definition of the gain margin, GM of the loop transfer function, 
L(s ); L(joj) = - 10 GM for some frequency, oj. In this case the magnitude of the 
return difference is given by 

|1 + L{ju))\ = |1 - 10 GM | 

and the magnitude of the inverse-return difference is given by 

|1 + L~\ju)\ = |1 - 10~ GM | 

Thus clearly a small gain margin (say | GM | < 0.3) implies a small return 
difference and inverse-return difference magnitude. Historically the gain margin, 
GM as defined above is multiplied by 20 to give its value in dB. That is 
GM = 1 and GM = 20 dB are equivalent. 
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LQG Loop Shaping 

The process of obtaining satisfactory (in terms of Fig. 11.13) Bode plots of 
the singular values of the loop transfer function is called loop shaping. The LQG 
design methodology can be a remarkably effective tool for achieving the loop 
shaping demanded by Fig. 11.13. A detailed description of the manner in which 
LQG can be used to solve multivariable control problems is given in Reference 
24. Some of the properties of LQG loops which make the LQG methodology 
effective for loop shaping will be briefly summarized below. In addition to LQG 
properties, some algorithms for loop shaping and two simple examples will also be 
discussed. 

LQ Regulator 

The linear quadratic (LQ) regulator problem assumes that a model of the 
system 

x = Ax + Bu R" U€ R m 

and a performance index 

i r°° 

J — ~ J 0 {x T Qx + u t Ru} dt Q > 0 , R > 0 

have been specified where the objective is to minimize J. This results [Refs. 1-5] 
in the control law u — -K c x where K e is obtained from the positive semi- 
definite solution, P e of the Riccati equation: 
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A t P c + P e A 4- Q - P c BR- 1 B T P e = 0 

as K c — R~ 1 B t P c where it is assumed that ( A,B ) is stabilizable. 

It was shown in Reference 25 that there is no loss of generality in taking 
Q — H t H where H is an m X n matrix. That is for any Q,R the full state 
feedback matrix, K c could also have been obtained for some H where 
Q = H T H . Thus the designer’s only input to the LQ regulator problem will be 
taken to be H. To insure a stable regulator H is always taken to be such that 
(A ,H) is detectable. 

The LQ regulator’s block diagram is shown in Fig. n.14. It is well known 



Fig. H.14 LQ Regulator Block Diagram 

[Ref. 1-5] that if ( A ,H ) is detectable the regulator is stable. Other theoretical 
properties of interest here can be derived from the Riccati equation as follows. 


For notational convenience let 
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4>{s) A (s/-^4) 1 
L,-(s) J K e (sI-A)~ 1 B 
VK e (s) £ H{sI-A)~ 1 B . 

Note that L,-(s) is just the regulator’s loop transfer function and lT c (s) is the 
transfer function corresponding to the weighting in the LQ performance index 
which the designer can specify. Starting with the Riccati equation we have 

A t P , + P,A + H t H - P,BB r P, = 0 
(~juBA T )P, + P,\jul-A) - H T H + KjK t = 0 
PJ+ <t>"P, - 4> b H t 1U + <fr H KjK'<fr = 0 
B T P,4,B + B t < t> B P c B - B T <f> 3 H T H<t>B + B T <t> B KjK'4,B = 0 
K C *B+ B t 4, b KJ- W B (ju) W'Uu)+ LfUuVAM = 0 
[/+ L,Uu)) = /+ W%ju)W'U») (II. LQ) . 

Two important properties of the LQ regulator can be determined from this 
final equation and they are 

LQ Property 1: g\I + L { (ju)\ >1 \fu 

LQ Property 2: a[L { (ju)] = a[W e (joj)\ 

when e[W e (ju))] » 1 

The first can be seen by observing that the right hand side of (II.LQ) has 
eigenvalues which are all greater than or equal to unity (since 
W^(ju) W e (ju) > 0 Vw) together with the definition of singular values 
applied to the left hand side. The second can be seen be observing that the 
identity matrices are negligible when fz[ W e {ju)\ is large (i.e. > > 1 ). 
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The first property is important since it implies that the LQ regulator’s input 
return difference has a minimum singular value which is larger than or equal to 
unity for any frequency. This in turn implies that the regulator has desirable 
disturbance rejection and desensitization properties (i.e. both are equal to or 
better than open loop for any frequency!). 

The second property is important since the command and disturbance 
requirements of Fig. 11.13 are for large loop gain. Thus the regulator loop 
transfer function’s high gain characteristics can be specified a priori by choosing 
the LQ weighting matrix, H such that U^fs) = H(sI-A)~ 1 B has the required 
high gain characteristic. 

To achieve the high gain performance requirement with H{sI-A)~ x B it may 
also be necessary to append additional dynamics. For example, to achieve zero 
steady state errors may require additional integrators in the plant (i.e. integral 
control). This is equivalent to frequency dependent weighting [Ref. 26]. 

It is also well known [Ref. 27, 28] that the LQ regulator has certain 
guaranteed stability robustness properties. This is also a consequence of 
q\I +■ L { (ju)] > 1 \/uj. It is the multivariable generalization of avoiding the 

-1 critical point on the Nyquist diagram for SISO systems. It implies that LQ 
regulator loops provide reasonable transition or “crossover” between the low and 
high frequency regions shown in Fig. 11.13. 
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Finally we note that at high frequencies the LQ loop transfer function 
approaches [Ref. 24] 

»IMM ] S - H\HB\ 

OJ 

when HB is full rank. This shows how the high frequency roll-off characteristics 
are related to H. This is a relatively slow attenuation rate and is the price the 
regulator pays for its excellent return difference properties. We recognize that no 

physical system can maintain a — characteristic indefinitely [Ref. 29], This is 

not a concern since L,(s) is only a design function and will have to be 
approximated by the full state loop transfer recovery procedure. 

Full State Loop Transfer Recovery 

The next step of the design process is to provide estimates of the states by 
processing the output measurements with a Kalman filter. This procedure is also 
well known [Refs. 1-5]. It involves the model of the system: 

i = Ax + Bu - f £ art R n u«R m 
y — Cx 4- t] ye R p 

where £ and r) are uncorrelated white noise processes with spectral intensities 
given by 


£[£(*)£ r (»)] = where M > 0 

E\rj{t )r} T {r)] = 6{t-T)N where N > 0 
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the estimation equations 

i = Ax + Bu + Kf(y-y) 
y = Cx 

and the performance index 

J = lim tr £[(z(0 -*(<)) WO -*(0) r ] 

( — >00 

where J is to be minimized by choice of Kj . 

The solution for the constant matrix, Kj is well known [Refs. 1-5] and is 
obtained from the positive semi-definite solution, Pj of the Riccati equation: 

APf + PfA T + M - Pf C t N- 1 CPj = 0 

as Kf = P/ C T N" 1 where it is assumed that ( A,C ) is observable. 

Using the estimates of the state in the control law i.e. « = -K e x results in 
the compensator transfer function 

K(s) = K c {sl-A + BK e + KfCY'Kf 

It is also well known that this compensator results in a stable closed loop system 
for any noise parameters: M,N. For the purposes of this thesis these noise 
intensities will be treated as design parameters which can be manipulated by the 
designer and not as some sacrosanct noise intensities for the system. By duality 
with the LQ results [Ref. 25] there is no loss of generality in taking M — TF 7 
and N = I where f is an n X m matrix. 
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For full state loop transfer recovery, the following value for T will be used: 

r = qB 

where q is a scalar design parameter. Then as q becomes large (assuming 
G(s) A C(sI-A)~ 1 B is minimum phase) it has been shown in Reference 30 that 
the filter gain behaves in such a way that 

lim K{s)G(s) = K e {sI-AY x B 

q-+ oo 

where the convergence is pointwise in s. 

This design procedure essentially inverts the plant from the left i.e. 

lim K(a) = K e (sI-A)~ l B G~ l (s) . 

^ “■♦OO 

This inversion, intuitively is why (^(s) is not allowed to have zeros in the right 
half plane. In practice the recovery procedure is effective as long as G(s) has no 
right half plane zeros with magnitudes in frequency ranges where high loop gain 
is required. The limitations on the achievable performance of feedback systems 
because of non-minimum phase zeros is discussed in Ref. 31. 

It has been suggested [Refs. 32,33] that an improved recovery is obtained by 
using colored rather than white process noise. This procedure however does not 
recover the LQ loop transfer function i.e. 

lim K{s)G(s) ^ K e {sI-A)~ x B 

q — oo 


when colored noise is used. The reason the authors of References 32 and 33 
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concluded that an improvement was obtained is that although their LQ loop 
transfer function, K c (sI-A)~ 1 B did not satisfy the requirements of Fig. 11.13, 
their LQG loop transfer function K(s)G(s) (obtained with colored noise) did 
satisfy the requirements of Fig. 11.13. The point is that the desirable /tT(s)G'(s) 
loop transfer function could have been recovered with white noise if the LQ loop 
transfer function had in fact been desirable. 

SISO Double Integrator Example 

The design methodology will now be illustrated with a simple example. 
Assume the desired loop shape must satisfy the requirements shown in Fig. 13.15 

where (7(s) = — . The state space matrices are 
s 2 



Fig. H.15 Desired Loop Shape For Double Integrator Plant 
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given by 

a ro n d roi p r, 

A ~ Lo oj B ~ LU 0 - l 1 °J * 

The first step is to choose H such that H(sI~A)~ l B satisfies the low frequency 
requirement. This is obtained by choosing 

H = K 2 0] 

as can be seen from the plot in Fig. 11.16. 

Solving the LQ Riccati equation gives 

K e — v/2 o, 0 ] and 

the LQ loop transfer function is given by 

1 

K e (sI-A)-'B = -f 

That it satisfies the requirements of Fig. 11.15, can be verified by examining the 
plot of | K c {jujl-A )~ l B | shown in Fig. 11.16. LQ Property 2 can also be verified 



by comparing 
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Fig. 11.16 LQ Loop Transfer Function for Double Integrator 


| H(ju>I-A )~ l B | and | K e (ju>I-A)~ l B | for frequencies ( v « u> 0 (i.e. where 
| H(jul-A )~^B | » 1). 

Solving the filter Riccati equation gives 


K f = 


V2q 

Q 


Computing the compensator transfer function results in 
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s + 


w 




K{s) = s/2u 0 sfq (u 0 +\/q) 


y/2 (u 0 + s/q ) 


p 4- s/2 (w,+ >/q)a + (w 0 + s/q f j 


Finally the input loop transfer function is given by 


A»C(s) 


U* q 

[. 

OJ 0 s/q 

1 

(u> 0 + y/j f 

s/2 (u 0 + s/q ) 

6 2 


O 

3 




[ n/2 (w 0 + s/q ) j 



K + y/q f 

[s 2 + s/2 (w # + s/q )* + (u; tf + s/q ) 2 ] 

and its magnitude for s — ju: is sketched in Fig. II. 17 for several values of q. 
The LQ loop transfer function is also sketched in Fig. 13.17 for comparison. 
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Fig. H.17 LQG Recovery for Double Integrator 

From the expression for K(s)G(s) above it can be seen that 

lim *(«)<?(«) = 4 .. 

?—oo s 6 U 0 

y/2 

The convergence of K(s)G(s) to K e (sI-A)~ 1 B as q goes to infinity occurs 
pointwise on s not uniformly on s. This means that for finite q the 
approximation A"(s )G(s) = K e (sl-A )~ l B will be valid only over a restricted 
frequency range. Outside that range the approximation can be quite poor. This 
is evident from the plots in Fig. H.17. Practically of course the errors at high 


« + 


u 

w o 

n/2 


= KAsI-AY'B 
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frequency cause little concern provided j K(joj) G(ju) | stays small. 

Designing a full state LQ regulator and then recovering the full state design 
with the Kalman filter has just been discussed. This approach was followed 
because it is the usual sequence one considers for LQG design. This process 
allows the designer to shape the input loop transfer function itT(s)G(s). 

A dual procedure for shaping the output loop transfer function G(s)A(s) is 
to design the filter first and recover with the regulator. The equations for this 
alternate procedure are mathematical “duals” of the ones given above. The 
subtle differences between the two procedures are discussed in Reference 31. 

The filter loop transfer function, L 0 (s) = C(sI-A)~ l Kj shown in block 
diagram form in Fig. II. 18, enjoys the same properties as the LQ loop transfer 
function. These properties are summarized in Table 13.3. 



Fig. H.18 Filter Loop Transfer Function 





Table II.3 Regulator and Filter Loop Transfer Function Properties 
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The duality of these two design procedures is easy to see in the SISO double 
integrator example given above by interchanging the roles of a and q. That is 
since they appear symmetrically, let u> 0 go to infinity while holding q fixed to 
give the desired loop shape (either /f(s)G(s) or G(s)K(s) since they are the 
same for SISO). 

The LQG input and output loop shaping procedures are summarized in 
Table 11.4. 

Advanced Loop Shaping 

The LQG loop shaping procedures just discussed required that the plant 
have the dynamics of the desired loop shape i.e. the poles of either K e (sI-A)~ 1 B 
or C(sl-A)~ l Kj are the same as those of the plant. Actually, as already alluded 
to, dynamics may be appended if they are not already present in the plant. 
Advanced loop shaping just formalizes this process. 

Assume the high gain characteristics of the desired loop shape are given by 
C[(sl-Ai)~ l B l i.e. A[, B t , C t are all three specified by the designer to satisfy the 
high gain requirements of Fig. 11.13. The plant transfer function will still be 
given by C(sI-A)~ l B . This desired loop shape can be approximately obtained 
for either the input or output loop transfer function with the appropriate one of 
the following two dual procedures. 

The procedure for the input will now be discussed and the other will follow 
by duality. The augmented plant state space realization is given by 



Table II. 4 LQG Input and Output Loop Shaping Procedure 


Input 


append additional dynamics 
if necessary and choose H 
or T such that the indicat- 
ed transfer function has 
satisfactory high gain 
characteristics 


solve^the Riccati and 
gain equations 


to obtain the ideal 
loop transfer function 


choose a scalar q and 
solve the Riccati and 
gain equations 



Output 


C(sl-A) _1 r 


AP f + p f A T + rr T - p f c T CP f = o 


K f = P f C' 


C(sl-A) 1 K f 


A ! P o + P A + q 2 C T C - P BB T P 
c c ^ c c 


K = B P r 
c c 


then if G(s) = C(sI-A)"'b 
is minimum phase 


K(s)G(s) -> K c (sI-A)"'B 


as q -v 
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A . = l A 0 

10 A, 

Q = [C 0] 


Bi 


B 

Bi 


Note that this does not change the input-output transmission of the plant i.e. 

G(s) = C(sI-AY l B = C { [sI~A { )- 1 Bi . 


Now to obtain an LQ loop transfer function which has the high gain 
characteristics of Ci(sl-A [)~ 1 B t let 

Hi = 10 C, 1 

and solve the regulator Riccati equation for P e > 0: 

Ajp c + F'A c + ElHi - P,BjE?F e = o 

to obtain the full state gain 

K = B?p e . 

When A is stable K e is given by 

K = [0 K e ) 

with K c — B[ T P e where P c > 0 is the solution to the (lower order) Riccati 
equation: 

A( T P c -I- P c A t + C t T C[ - P e B[ BfP e = 0 

Note that when the plant is stable it is not involved in the computations so far. 
That is only A h B h C t (typically of lower dimension than A,B,C) have to be 
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manipulated. 

Next the full state feedback recovery procedure is used as follows. Let 

r i = q B; 

where q is the scalar design parameter used for recovery and solve the filter 
Riccati equation for Pj > 0: 

A { P f + P f A, T + T { T, T - P f CfC { F f = 0 

to obtain the filter gain 

K f = PjCj . 

The compensator transfer function is obtained as usual from 
K(s) = K^I-Ai + E { R C + K '/C,)- 1 K, . 

Assuming G(s) is minimum phase and stable we have 

lim K(s)G{8) = K'isI-AtY 1 Bt 

q-+oo 

where cr[K e (j B{] = cr\C[(jujI-Ai)~ l B t ] for cu such that 

(T[C[[jLoI-Ai)- 1 B[] » 1 which was the design objective. 

The comments regarding the minimum phase assumption in a previous 
section covering ordinary full state feedback recovery apply here as well. That is 
the desirable properties of the LQ loop function K e {sl-A l )~ x B[ will be recovered 
for frequencies other than those near the magnitudes of right half plane zeros. 
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Using the procedure for an unstable system results in a loss of accuracy in 
the approximation of the desired loop shape. This is consistent with the fact 
that [Ref. 29] arbitrary loop shapes are not possible for unstable systems. Thus 
the procedure is still applicable to unstable systems but the choice of A { , B h C t 
must be made with attention to the limits to achievable performance. 

SISO First Order Example 

A simple example will now be discussed to help clarify the steps of this 
procedure. Let the plant be given by 

<?(«) = *7~T’ A = -1> B = 2, C = 1 . 

S *7* X 

Let the desired loop transfer function be given by 

C^sl-Ai)-' B t = — Ai = 0, B l = 100, C { = 1 . 

8 

The plant is stable and K e can be verified to be given by 

K, = 10 K,) = [0 1] . 

This results in the desired loop transfer function being achieved exactly by the 
LQ loop transfer function i.e. 

K c (sl-A { )- 1 B t - — = Cfal-Atr'Bt . 

s 

The filter gain can be verified to be given by 
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K, = 


v/I+4? 5 -1 

200 q 2 
\/T+4q 2 


from which it can be seen that 


lim Kt = 

g —*co 


2 ? ' 
Uoo 9 . 


Using this limiting expression for Kj to compute K(s) it can be seen that 


lim K(s)G{s) 


— lim 


100q(s + 1) 


oo s 2 4- (101+2g)s + 1 + 1 

200 


lim 

ii+ U2L+ 2 

? I 9 

100 

s 


s + — 
9 


As already alluded to, this procedure has a dual for shaping the output loop 
transfer function G(s)K(s). Both procedures are summarized in Table H.5. 

The design procedures just discussed are very powerful for shaping either the 
input or output open loop transfer functions. However, in most design problems 
both loops must have desirable properties. Thus, one major weakness of the 
design procedures is that only the properties of one open loop transfer function 
may be formally tailored. Of course for SISO problems the input and output 
loop transfer functions are the same since (7(s) and K(s) commute. 




H, = O’ V 


xlp +P X+hTHj-P b.bJp„ = 0 

1C Cl 11 C11C 


^jAP^ J +TT J -f f C J C P. = 0 
of fo 00 foof 


solve the Riccati 
and gain equation 
to obtain either 

Regulator Loop Transfer: T< (sl-ft. 

or _ i 

Filter Loop Transfer: C n (sl-A n ) Kf 


P c >0 


t 


K = b]p„ 

c 1C 


A stable 


P f > 0 


I 


K f ■ Vo 

A stable 


A I P c +P cV C I C r P c B £ B I P c • 0 


P c >0 K c =[0K c ] K c =B{P c 


A £ P f +P f A £ +B £ B £" P f C .e C £ P f = 0 


P f > 0 K f = 


K f = PfCf 


tn 

vo 


solve the Riccati 
and gain equation 


A. V P f X{ +q 2 Bi BT-P f cT C .P f .o 


A T P +P ft +q 2 C T C -P B B T P =0 
oc C 0 ^ 00 cooc 


Pf > 0 


Kf - P f C { 


P > 0 
c — 


K = b]"p 
c o c 


Form K(s) 


K(s)=K c (sI-A i +B i K c +K f C.)' 1 K f 


“1 17 


K(s) = K (sl-ft +B K +K-CJ"'K- 
' ' c 0 0 c f o' f 


then if G(s) 
is minimum 
phase 


lim K(s)G(s) 

q-x» 


KjsI-A.J-'B. 
,-l 


lim G(s)K(s) = C 0 (sI-A 0 ) _l K f 


= K c (sl-A^)~ B^ for A stable 


q-*» 


= C^(sl-A^) - ^ Kf for A stable 
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Motivation For Model Reduction 

As we have seen a powerful design technique ( LQG loop shaping) exists to 
formally solve the control problem. However, this design technique requires the 
manipulation of matrices with dimensions greater than or equal to that of the 
plant model. Also it results in a compensator transfer function which has order 
equal to or greater than that of the plant. 

For many problems the order of the plant is so high, as to prohibit a 
successful design of the compensator (either by cost of computations or accuracy 
of computations). Even when the design process can be carried out it may be 
desirable to simplify the compensator (i.e. reduce its order) for implementation 
purposes. This could be to save implementation computer time and memory, 
provide ease of checkout and verification of the implementation, eliminate 
excessive gain scheduling, or for any of many other practical reasons. 

Thus we are led to search for reduced order models to simplify the design, 
analysis and implementation of the control system. Many methods of obtaining a 
reduced order model exist [Refs. 6-15] and the fundamental ideas of the 
internally balanced realization method will be discussed next. 

Internally Balanced Realizations 

Balanced realizations will be used throughout the sequel. When the 
adjective internally is used, it represents the essentially unique realization of an 
asymptotically stable MIMO transfer function defined by B. C. Moore [Ref. 11]. 
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A geometric interpretation of the balanced realization, an algorithm for 
computing the realization and some of its properties will now be discussed. 

A state space realization of the system to be balanced is given by 

x = Fx + Gu xdt N ucR 1 " 
y — IIx ye R p 

It is assumed that the system is asymptotically stable, i.e. /?e[X[F]] < 0. The 
minimal order of the system will be taken to be n where n < N. That is the 
given system may be uncontrollable and/or unobservable. 


Geometric Interpretation 

With reference to the block diagram shown in Fig. 11.19 ask the following 

two 


u(t) 




x(t) 

r — | 


G 


(Sl-F)' 1 


JLT 


Fig. 11.19 Block Diagram for Internally Balanced Realization Discussion 

(dual) questions. What set of points in the x-state space could be part of the 
zero initial condition response for some input, u(t) such that ||«(/))| 2 < 1? and 
what set of points in the x-state space as initial conditions could produce an 
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output, y(t) such that ||y(?)|j 2 A 1 with zero external input? 

Moore showed that these two sets are in general different ellipsoids. The 
following general theorem was proved by Moore, [Ref. 11]. 


Grammian- Ellipsoid Theorem 

Let F(t) be an impulse response matrix of some asymptotically stable, 
linear, time invariant system. Let 


G f A j o F(t)F T {t) dt > 0 


S c A x(t) = j g F{t-r)u(T)dT, \/t, ||«(/)|| 2 < lj 

5, A |ar(o): y{t) = / 0 F T (t-T)x{o)6{T)dT, \ft , ||y(/)|| 2 = 1 


the grammian, controllable set and observable set for F(t) respectively. 
Furthermore let the eigenvalue decomposition of the grammian be given by 

G f = V T V — !,'£ = diag {cr, }, V = block col[t--] . 

Then 

S c = ellipsoid with semi-axes given by 

and 

S g = ellipsoid with semi-axes given by — — v { 

For the first question let F(t) = e Ft G , then the answer is the controllable 
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set, S c . For the second question let F(t) — e fTt H T , then the answer is the 
observable set, S 0 . These grammians have been given special names [Refs. 11, 
34, 35]: 

r 00 r 

U = J ^ e Ft GG T e Ft dt = controllability grammian 
Y = J 0 e FTt H T He Ft dt = observability grammian 

It is well known [Refs. 11, 34, 35] that these grammians can be computed from 
the Lyapunov equations 

FU + UF T + GG t = 0 (D.la) 

f t y+ yf + h t h = o (n.ib) . 

These results are now summarized. The lengths and directions of the 
controllability ellipsoid semi-axes are the square roots of the eigenvalues and the 
eigenvectors of the controllability grammian respectively. The lengths and 
directions of the observability ellipsoid semi-axes are the reciprocals of the square 
roots of the eigenvalues and the eigenvectors of the observability grammian 
respectively. 

It is well known [Refs. 34, 35] that a zero eigenvalue of the controllability 
grammian implies {F\G} is uncontrollable and that a zero eigenvalue of the 
observability grammian implies {F,H} is unobservable. The association of these 
grammian eigenvalues with the lengths of the ellipsoid semi-axes is intuitively 
pleasing in that: an uncontrollable direction would clearly correspond to an 
ellipsoid semi-axis of zero length and an unobservable direction would clearly 
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correspond to an ellipsoid semi-axis of infinite length. Hence the eigenvalues of 
these grammians provide scalar measures of how controllable or observable a 
given grammian eigenvector direction is. 

Moore also showed that there exists a realization of the system obtained 
with a change of variables by a similarity transformation, T, such that the axes 
of the ellipsoids are the same (balanced) for the new state variables (say z, where 
x = Tz). For model reduction, this balancing is the key idea because it provides 
a basis for the n dimensional x-state space where the direction of a given basis 
vector is as controllable as observable in a well defined sense. The lengths of the 
ellipsoid’s semi-axes provide a scalar measure of how controllable and observable 
a given basis vector direction is. Finally the reduced order model is obtained by 
neglecting the weakly controllable/observable states of the system. 

Algorithm For Computing the Internally Balanced Realization 

There are in fact several algorithms for computing the internally balanced 
realization. Most suffer from numerical difficulties when they are applied to non- 
trivial problems. The difficulties arise due to the squaring up nature of the 
problem i.e. the G and H matrices are squared to compute the grammians and as 
will be seen the controllability and observability grammians are then multiplied 
together. This is an important research area and progress has been made by 
Laub at USC however these results were not used for this thesis. 
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Tliie algorithm for computing the internally balanced realization that will be 
presented here is slightly more general than Moore’s in that {F,G,H} may be 
uncontrollable and/or unobservable. The objective of any balancing algorithm is 
to find a minimal (i.e. controllable and observable) realization, {A ,B, C) of the 
given possibly non-minimal realization, {F,G,H} such that the controllability 
and observability grammians for {A,B,C} are equal and diagonal. 

To this end consider how the original grammians change due to a similarity 
transformation. Starting with (0.1 ab) we have 

T~ x FTT- 1 UT~ t + T~ 1 UT~ t T t F t T~ t + T~ 1 GG t T- t = 0 
t t f t t -t t t Y t + t t YTT~ x FT + T t H t HT = 0 

which can be rewritten as 

FU + UF t + GG t = 0 
F t Y + YF + H t H = 0 

where 

F A T~ l FT 
G A T~ l G 
H A HT 

0 A t~ x ut- t 
Y A T t YT 


Now our objective can be stated as follows: find T such that 
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0 = T-'UT~ t = block diag {£,£,} 

Y= T t YT = block diag {£, £„} 

where 

E — diag {crj, cr 2 , ex n } 

E u E y = diag <rf +2 , a%) 

with (7 X > cr 2 > • • • > <r„> cr n+1 = cr n + 2 = • • • = a N — 0. The 

uncontrollable and/or unobservable modes result in E u E y = 0. 

The similarity transformation, T which accomplishes the objective is 
obtained from the eigenvector decomposition 

UY = TAT- 1 

That is the columns of T are eigenvectors (eigenvector nonuniqueness discussed 
below) of UY corresponding to A = diag{\,-} the real diagonal eigenvalue 
matrix of UY (assuming without loss of generality that the X/ s have been 
ordered such that Xj > X 2 > • • • > \ N ). The fact that UY has a real 
diagonal Jordan form is a consequence of the fact that U and Y are both positive 
semi-definite. It is also true that A is positive semi-definite. The proofs of these 
facts are not well known and are rather long, hence they are contained in 
Appendix C. 

Eigenvectors are not unique since a scalar times an eigenvector is also an 
eigenvector and when the eigenvalues are not distinct linear combinations of 
eigenvectors corresponding to a repeated eigenvalue are also eigenvectors. 
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Specific eigenvectors can always be chosen such that 

T~ 1 UT~ t = block diag{E, E u } 
T t YT = block diag {E, E„} 


where 


E = diag {cr j, (t 2 , <r n } 
E u E y = diag a% +2 , 


with 

> ^2 ' ' ' > > <T n+ 1 ~ ° r n+2 == ' ' ' ~ a N = 0 

The details of this choice of eigenvectors is messy and hence is also contained in 

App. C, 

Fortunately any choice of eigenvectors will work for the purpose of model 
reduction. The eigenvector scaling just leads to a choice of scale for the 
individual balanced state variables {the reduced order model transfer function is 
independent of this scaling!). The complication, due to repeated eigenvalues of 
UY, is eliminated by either retaining or eliminating all the balanced states 
corresponding to the repeated eigenvalue of UY in the reduced order model 
(there is no justification for doing anything else!). 

The square roots of the eigenvalues of UY are the singular values of the 
balanced grammian. That is 

a i = \A7 « = 2, N 
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Let T n be the first n columns of T (i.e. those columns which correspond to non- 
zero balanced grammian singular values), likewise let S n be the first n rows of 
T ~ l , then the n th order minimal internally balanced realization of {F,G,H} is 
given by {A,B,C} = {S n FT n , S n G, HTJ. That is 

AY, + Y aT + BbT = 0 (II. 2a) 

A T Y + T, A + C T C = 0 (H.2b) 

where 

S = dia S (r 2. <L,} > 0 . 

Moore showed that the internally balanced realization is essentially unique 
when the balanced grammian singular values are distinct. Essentially unique is 
taken to mean unique up to a change in sign of a state variable. 

The algorithm for transforming a given asymptotically stable, possibly non- 
minimal realization {F,G,H} into an internally balanced, minimal realization 
{A,B,C} is summarized in Table II. 6. 

Simple Example 

A simple example will now be used to fix ideas. Consider the following 4th 
order non-minimal realization with one controllable and observable mode, one 
controllable but unobservable mode, one uncontrollable but observable mode and 


one uncontrollable and unobservable mode. 
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Table II. 6 Internally Balanced Realization Algorithm 

Given: F,G,H with Re [X[F]] <0 

Solve for U and Y from 

FU + UF t + GG t = 0 
F t Y + YF + H t H = 0 . 


Solve for eigenvalues and eigenvectors of UY i.e. 

UY = TKT~ l . 

Partition T, A and T~ l such that E > 0 

'S n | 

&N- n . 

Compute A, B, C with 

A = S n FT n , B = S n G , C = HT n 


— \r T 
L n 


1 

N ~ n J 


Y 2 0 

T-l — 

0 0 

i — 


Then AY, + YA T + BB T = 0 

a t y+ ya + C T C = 


0 
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-1 

0 

0 

O' 


T 

0 

0 

-2 

0 

0 

-3 

0 

0 

G = 

1 

0 

.0 

0 

0 

-4. 


0 


II = [1 0 1 0] 


The two grammians can be verified to be given by 


— — 00 

2 3 

i i 0 0 


— o — o 

2 4 

U = 

3 4 

0 0 0 0 

0 0 0 0 

J 

r = 

0 0 0 0 

0 1 0 

4 6 

0 0 0 0 


The similarity transformation or the eigenvector matrix of UY can be verified to 
be given by 


T = 


1 0 -T 0 

-10 0 
3 

0 0 10 
0 0 0 1 


The grammian singular values can be verified to be given by 

E = ~ 

^ 2 

E« = dia s °> °) 

E y = dia s {°> ^ °) • 

Finally a minimal realization is given by 

A = -1 B = 1 C = 1 . 
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Some Properties of the Internally Balanced Realization 

The internally balanced realization has several fascinating properties. The 
properties of interest here are stability, controllability and observability of 
reduced order models obtained from the balanced realization. 

As mentioned earlier a reduced order model is obtained from the internally 
balanced realization by neglecting the weakly controllable/observable states of 
the balanced realization. The singular values of the balanced grammian provide 
a measure for determining how controllable/observable a given state direction is. 
The idea is that the states corresponding to the smallest singular values can be 
neglected. This is expressed by the following. Choose the order of the reduced 
order model, r such that a r > <r r+1 then let = diag {<Tj, a 2 , ..., c r ) and 
]T] o = diag {<T r+1 , <r r+ 2 , ..., <7 n }. Partition A, B, C compatibly as 

A _ -^11-^12 

A 21 A 22 

c = [C, cy 

then {A n B lt C 1 } r is the r th order reduced order model of {A , B, C} n . 

Immediately from (II. 2 a,b) it can be seen that the realization of the reduced 
order model is internally balanced i.e. 

•^■11^1 "h = 0 

A y Sj + SjAjj + = 0 

Moore showed that the reduced order model is genericallv asymptotically 


B 


B i 
B 2 
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stable, controllable and observable. Pernebo and Silverman [Ref. 36] proved the 
following stronger result. The condition that < 7 r 7 ^ cr r + j implies that the reduced 
order model will be asymptotically stable, controllable and observable. 
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m. FREQUENCY WEIGHTED MODEL REDUCTION 

One of the main points of this thesis is that model reduction and control 
system design are not independent of each other. That is, given that a control 
system will be designed, analyzed and implemented based on the reduced order 
model, the technique used for model reduction must be cognizant of this fact. 
The point is that the reduction will introduce error and the criterion used to 
define a small error must reflect the purposes for which the reduced order model 
is intended. 

Model Fidelity With Respect to Control System Design 

The overriding concern in control design is stability of the closed loop 
system. The next concern is that the performance objectives are satisfied. In 
Chapter II, it was shown that for good performance and robustness a loop shape 
must be obtained which has high gain (typically at low frequency) to satisfy 
disturbance attenuation, desensitization and command response specifications and 
low gain (always at high frequency) to satisfy sensor noise and robustness 
specifications. Furthermore the loop must accomplish the transition (crossover) 
between these two regions in a stable manner. Therefore the model reduction 
error criterion must assess how much the use of a reduced order model can affect 
the desired loop shape. 

In the low frequency range, the feedback compensator should provide 
adequate loop gain, so the high gain of the loop shape will not be seriously 
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degraded by some low frequency model error. In the high frequency range the 
loop must be rolled off such that stability robustness to uncertainty in the full 
order model is maintained. Thus an accurate approximation of the full order 
model is not necessary in the high frequency range. This leaves the mid 
frequency or crossover region as the only critical region for accurately modeling 
the plant dynamics. 

To assure stability it is also clear from the Nyquist stability criterion (and 
its MIMO generalization) that the reduced order model must not neglect unstable 
poles. Thus another requirement for model fidelity is that the reduced order 
model have the same number of unstable poles as that of the full order model. 

These comments can be interpreted graphically for the SISO Nyquist plot 
shown in Fig. III.l for the loop shape requirements: 

\L(ju)\ > R w < U[ 

1 1 + L (ju) | > 1 - r < u> < u> h 

| L(ju>) | < r w A < uj 

where L(s) — G(s)K(s) — K(s)G(s) and G(s) and are the plant and 

compensator transfer functions respectively. The solid line is the locus of 
G r (ju)K(ju>) versus ui i.e. the loop is analyzed and designed with the reduced 
order model transfer function G r (s). 

It can be readily seen from Fig. III.l that if G r (joj)K(ju>) satisfies the loop 
shape requirements and the error between G(ju) and G r (jui) is such that the 
locus of G(joj)K (ju) versus u> lies in the shaded region, then the loop shape 
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Figure 1 1 1.1 Allowable Reduced Order Model Error Bound 
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requirements will be satisfied for G(ju)K(ju) as well. Thus C r (s) is an 
acceptable reduced order model and the frequency dependent width of the shaded 
region is the frequency dependent allowable model reduction error. 

Motivation For Frequency Weighted Model Reduction Error Criterion 

The above comments have already indicated that for control design the 
model reduction error criterion must take into account the frequency dependence 
of the allowable error. This rather intuitive discussion will be made more formal 
now. 

Given that the model reduction method will not neglect unstable poles and 
that a stable closed loop, system can be designed for the reduced order model, the 
stability of the full order model closed loop system can be assessed with the 
stability robustness theorem of Chapter II. 

Consider the full order model closed loop system in block diagram form 
shown in Fig. III. 2 where G(s), G r (s) and K(s) are the transfer functions for the 
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full order model, reduced order model and the compensator respectively. Note 
that, the block diagram is really just a rearranged version of the compensator 
and the full order model in a standard feedback configuration (i.e. G r (s) cancels 
out). Applying the stability robustness theorem to this configuration implies that 
the full order model closed loop system is stable if either 

or 

||^( 5 )[G( S )-G r ( S )]|| 00 <l 

where 

WM A *•(»)[/+ G,(s)K(s )]-' 

KM MU + «■(•) ■ 

It is of interest to compare the magnitude of these weighting transfer 
functions with the intuitive discussion given above. For simplicity consider the 
SISO special case. Let 

W(s) = W { (s) = W 0 (s) 

= K(.)|l + G f («)K(.)]-' 

= Gr,(«)K(«>[l + C? r («)A-(-)|-' . 

Now then G r (jw)K(ju>)[ 1 + G r (ju)K(ju>)]~ i (the closed loop system frequency 
response) will typically have a magnitude as shown in Fig. HI. 3 (i.e. 
| G r (ju>)K(juj) | » 1 for small w and | G r (ju>)K(ju>) | « 1 for large w). 




Fig. HI. 3 Typical Closed Loop System 


Next the reduced order model transfer function, C r (s) will commonly have a 
magnitude plot like one of those shown in Fig. HI. 4. 



Fig. EH. 4 Two Common Reduced Order Models 
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Using Fig. HI. 3 the magnitudes of the weightings corresponding to the two 
plants in Fig. III. 4 can be computed and are sketched in Fig. III. 5. 


log IW (jco) I 



Fig. HI. 5 Two Common Weightings 

It can be seen that these weightings indeed emphasize the crossover region 
and place less weight on low and high frequency error depending on the reduced 
order model DC gain. Thus the intuitive discussion is consistent with the formal 
stability analysis. 

Definition of Frequency Weighted H ^ Model Reduction Error Criterion 

It has been shown that for stability, the frequency dependence of the model 
reduction error is important. That is, errors should be small in some frequency 
ranges (crossover) and can be larger in other ranges (low and high). This 
motivates the use of a weighted error criterion. 
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Stability also requires that an inequality be satisfied for every frequency, 
hence, it is the magnitude of the highest peak in the weighted error that is 
important (as opposed to some integral squared error criterion). This motivates 
the use of the norm. 

To obtain a formal definition of the error criterion let G(s), H^s) and 
^(s) be the given full order model, input weighting and output weighting 
transfer functions respectively. Then the scalar model reduction error, for 
the reduced order model transfer function, (? r (s) is defined to be given by 

£» & ||»,(.)[(?(.) -<?.(*)] w ( (.)ll„ • 

The purpose for which the reduced order model is intended (i.e. control 
design) is reflected by the choice of the weighting transfer functions. The 
remainder of this chapter will develop the model reduction technique assuming 
that the weightings are given. The choice of weightings is important and will be 
discussed further in Chapter IV. 

Model Reduction Problem 

The model reduction problem will be defined to be the following. Given an 
nth order state space realization (P, G, H} n of the full order model transfer 
function, G'(s) = H (si -F)~ l G and transfer functions for the input and output 
weightings, W'-(s) and W^s) respectively, find an r tk order ( r specified and 
r < n) state space realization {A n , J5 1; C j} r of the reduced order model 
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G r (s) == C X {sl-A n )~ x Sj such that is minimized and G(s) and G>(s) have 
the same number of unstable poles. 

In mathematical notation the problem is 

Given: {F , G, H} n , Mft*), W c (s), r < n 
Find: {A n , B h C x } r 

such that F and A u have the same number of closed right half plane eigenvalues 
and 

||lV„( s )[ff( s /-F)-‘C - C^I-A n )-' B£ »;-(.)ll„ 

is minimized. 

This problem appears to be intractable, however, the SISO special case with 
unity weightings and r = n-1 has been solved with the Hankel norm technique 
[Ref. 14], although the reduced order model is also required to have a constant 
feedthru term, D j, i.e. G r (s) = CrfsI-A jj ) -1 B x + J9j. 

Internal Balancing As An Approximate Solution 

Consider the unity weightings special case and reduced order models 
obtained from an internally balanced realization. This also requires that the full 
order model is asymptotically stable. For this special case the model reduction 


error criterion, can be bounded. 
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Internally Balanced Realization Model Reduction Error Bound 

Assume that the asymptotically stable, minimal realization {A, B, C} n of 
the full order model transfer function, G(s) = C{sI-A)~ 1 B , is internally 
balanced, i.e. 

AE + EA T + BB t = 0 
A T E + EA + C T C = 0 

where 

E — diag {cr,- } 
and 


ffj > <r 2 > • • • > a n 

Let the matrices A, B, C, be partitioned compatibly as 



C = [C, c 2 1 





B 



where the dimensions of A u and Ej are r X r. Assuming that a r+l > a r 
implies that the reduced order model G r (s) — Cj(sI~A jj) - 1 B j is asymptotically 
stable and thus satisfies the model fidelity requirement that G(s ) and G’ r (s) have 
the same number of unstable poles (zero in this case). 
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Then the model reduction error criterion is bounded: 

||C(.)-G r (*)IL<«rlE2l 

and for r = n-1 the bound is tight, i.e. 

I|G(«) )IL = 2cr n . 

The proof follows. Let 

m J (*i-A n r 1 

A(s) =r $I~A<v£ — A 2 i<f)(s )A 12 
B(s) J A 2 ^( S )B,+ B 2 
C(s) ^ C^(s)A l2 + C 2 , 

then 


<?(«)- G r («) 


C^Z-A)" 1 B - C i <l)(s)B 1 
[^i ^2] 


(si-Ajj) 

~A 12 

-1 

~A 21 

sI-A 22 





Using the inverse of a block matrix formula [Ref. 34] 
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G(s)~ <?,(«) = [C x C 2 ] 

<p{s) + <t>{s)A n A- 1 (s) j 4 2i «^(s) ^(«)^ 12 A -1 («) B x 
A _1 (fi)^ 21 0(s) A -1 (s) B 2 

- C 1 <P(s)B 1 

— C x 4>{s)A X 2& 1 ( s M 21^( s )'®1 + 1 ( s M21 < ^( 5 ) B x 

+ C x <p[s)A 12 A _1 (s)S 2 + C' 2 A" 1 (s)B 2 
= C(5)A- 1 (s)5(s) 

by algebra and the definitions of B(s) and C(s). 

From the definition of the maximum singular value we have 

a^G(juj) - G r (ju>)^ = - G r {ju)y^G(ju) - G r (;w)] J ] 

Substituting for G(ju) - G r (ju) in terms of B(ju), C(ju>) and A(;w) from 
above, we have 

u.)-C r (j«)] = XAfi [cU U )A-Hj U )BUw)B B Uw)A-> 1 Uu,)C B uA i , 

Using the fact that X max [;4i3] = X max [&4] for any matrices A, B where the 
products are defined [Ref. 37], we have 

»[<?(/«) - GAM] = >'^[A- , (jw)B(; W )B s (i W )A- ff (j W )d fl (ju,)C'(; U )] . 

Expressions for B(ju)B E (ju>) and C H (jui)C(j u) are obtained by using the 
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partitioned form of the internally balanced grammian equations 

A n Ej + Mn + == 0 

A joSo d* ^-'1^21 d" B \B% == 0 
J d-22“-'2 d" 1^2^22 d" B 2 B 2 ~~ 0 

A 11 ^1 d" + C'j^Cj == 0 

-^21^2 d- BiyA 12 + CfC t == 0 

^2^2 d B’jA 22 + C'JC'g == 0 

where the first three of the above equations were obtained from 
A S + TiA T 4- BB T = 0 and the last three were obtained from 
A t E + EA + C T C = 0. 

~ ~ TY 

An expression for B(ju)B (ju) is obtained as follows. By the definition of 
B[s), we have 

B(j<jj)B n {ju) = A 2l (f>{ju)B 1 Bj' ^ H {ju)A^ + A 2 ^{ju)B x Bl 
d- B 2 Bi <f> B (juj) A 21 + B 2 B 2 

Substituting for B^Bf, B X B 2 and B 2 B 2 from the partitioned form of the 


grammian equations, w'e have 
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+ A<n<j){ju)) ^4 i 2 E 2 4- E^^ j 

+ ^4 21 5U j — f- £2-4 {j w)A.2i + -4 22 E 2 d“ So-^oo^ 

— _ ^2i [ < ^(i a, )(-^n^i+ S^4 T\j<f> B (joj) + <j>(ju)Yi X 

+ Ej^ ff (;o;)j A Jj 

+ ^21^(i U )A 12+ ^ 22J ^2 + E 2 ^4 \2<!> B {jw)A2\ + ^oojj 

after rearranging terms. Substituting jwl - A(;u>) for A^U^A^ + 4 22 from 
the definition of A(s), we have 

B{ju)B H (ju>) — -|-4 2 i<^(iw) jU jjEj + Si Ah 

+ Si <f>~ H {joj) + ^'{ju^^ijujAZi 
+ J/w/-A(/w)jE 2 + E 2 [-> ul - A- ff (;w)j| 

where the first term was also further factored. Noting that the first term is zero 
by the definition of <^(s), we have 

B{ju)B n (ju) = A(;'w)E 2 + E 2 A ff (jw) 

The expression for C H (ju)C(juj) is obtained analogously and is given by 
C H (ju)C(ju) = E 2 A(;w) + A*(;w)E 2 . 
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These expressions for B(ju)B 3 (ju) and C B (joj)C(joj) are then substituted into 
the final expression for <r\G(ju>) - G r (j a;)] above to obtain 

?[G(,w)-G r (>u)] = Xi£[A- , (;u)[A(y«)E 2 + E 2 A*(y«)]A-®(H 

[s 2 A(j W ) + A fl (>V)£ 2 ]] 

= X*fi[[s, + A">(;a.)E 2 A fl (j W )] 

[A-'f/uJEjAt/w) + S 2 ]] 

after algebraic simplification. 

Consider order reduction by one state (i.e. r == n-1) then E 2 = c n and 

a[G{ju)~ G n _ x {ju)\ == cr n \Z[T+ A(j'w)][l + A-\ju)} 

where A (ju>) A A H (ju)A~ l (ju ;) is an “all pass” i.e. |^4(yw)| = 1 \fu. Thus 
after rearranging terms and making use of the “all pass” property of A(ju) we 
have 

v[G[ju)- G n _i(jw)] — o n 

= U + A(ju>) | 

By the triangle inequality we have 

a[G(yw) - < cr n [l + | A(;w)|] = 2 a n 

This completes the proof of the bound for r = n-1. 

The remainder of the proof is achieved by using the order reduction by one 
state result and by noting that {A n B j, C x } k obtained by the k th order 
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partitioning is internally balanced with balanced grammian, 
Ej = diag {a h <r 2 , cr k }. 

Let E k (s) A G k+1 (s)~ G k (s) for A: = 1, 2, n-1 and let G n (s) A G(s). 

Then 


*[£*(;«)] < 2 (r h+i 

since G'^s) is a reduced order model obtained from the internally balanced 
realization of G k+ !(s) and the above bound for order reduction by one state 
holds. 

Noting that 


G(s)-G r (s) = £ E k (s) 

k—r 

by the definition of E k (s), we have 


»r, 


*[G{jv)- G r {jw ) j = a 


n-1 


£ Ek(j“) 




Using the triangle inequality we have 


- G r (joj ) j < £ a[E k {ju ) j 


1 V 1 sT/r. 


* — r 


But a[E k (ju))\ < 2a k + j from above, hence 


v[g{}») ~ G r (/«)] < 2 £ <r* +1 = 2tr[E 2 ] 


by the definition of the trace operator and E 2 . This completes the proof of the 
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bound: 

5 [<?(>«)- c r (;w)] < 2ir[E 2 ] . 

Note that when r = n - 1 by the definition of A (s) we have 

A(0) := A 21 A n ^!2 _ -^22 

(the inverse of An exists since A n has eigenvalues with strictly negative real 
parts). But, from the determinant of block matrices formula [Ref. 34] 

det A 
det A u 

Neither det A or det A n is zero since both A and A n have eigenvalues with 
strictly negative real parts, thus we have 



0 < | A(0) | < oo 


and hence 


^<°> - m = 1 


Using this in the equality from above for u> = 0 we have 


»[<?«>) - c..,(0)] 


= »„|1+ A (0)| 

= 2<t„ 


Thus the bound is tight for order reduction by one state and is in fact achieved 
at DC, i.e. oj — 0. 
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Tightness of the Bound For Two Special Cases 

For the following two special cases the model reduction error criterion, E 
and the bound, 2tr[£ 2 ] can be computed with a limiting process and then 
compared. The first example has poles and zeros which alternate along the 
negative real axis and the bound is tight in this case. The second example has 
poles and zeros which alternate along the ju axis and the bound is not tight in 
this case. 


Alternating Real Poles and Zeros Example 

Consider the SISO system with transfer function given by 


G(s) = 


n n k 

£ “Hr 

* = 1 5 + f 


v > o 


For s real, the derivative of G(s) with respect to s is negative so the poles and 
zeros must alternate along the real axis [Ref. 38]. The results in Reference 39 
show that 


lim G r (s) 

p —►00 

lim E 2 

p — ►00 

lim HG'(s) - G T ( s ) 1 1 qo 

p— ►OO 


= E 


* + V 

1 


2 
n-r 


h-r 


Applying the bound and the above result for E 2 to this example results in 
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lim 1 1 O' ( s ) — G'r ( s ) 1 1 oo < 2 lim fr[E 2 ] == n-r 


p—+oo 


p —*oo 


Thus the actual error is equal to the bound in this limiting case. 


Alternating Imaginary Poles and Zeros Example 


Consider the SISO system with transfer function given by 


G(s) = £ 




*■=1 S z + 2 $u k s + 

where f > 0 and 0 < Wj < w 2 < • • • < u n . Note that f > 0 is required such 
that the transfer function is asymptotically stable, but 


lim <?<*) = E ' T U ‘~ 


has alternating poles and zeros along the imaginary axis (again, proved in Ref. 
38). With the further assumption that the w/s are sufficiently separated, i.e. 


fmax(o;,-,a;y) 
I W,- - Uj | 


« 1 V i,j — 1, n , ij£j 


the results of References 40 and 41 show that 


lim C r (s) 
f-o 

lim £ 2 


r uj; 

lira £ — — 

f-*° ft — i s s + 2 $u k 8 + 

1 


= — /, 
4f 


2(n-r) 


Substituting for C r (s) and by the definition of the H ^ norm we have 



sup 

w 


“k 


lim ^ — — 

f-*° *« r+ i wf - or + 


lim — 
f-*o 2 f 


Applying the bound and the above result for this example results in 


lim \\G(s) - G’ r (*)|| 00 < 2 lim tr^ 2 ] = “ (»-r) 


Thus the actual error is a factor of 2(n-r) smaller than the bound in this limiting 


case. 


The results for these two cases are indicated graphically in Fig. HI. 6. 



V. 


V-/ 


V-/ 
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Fig. HI. 6 Tightness of the Error Bound For Two Special Cases 
On the Question of Optimality 

As already alluded to, the Hankel norm technique does provide a reduced 
order model for r = n-1 and SISO which does achieve the minimal E ^ (with 
unity weightings). As was also noted this technique does require that the reduced 
order model have a finite D term, whereas the balancing technique uses a zero D 
term. With this extra degree of freedom it is expected that the £’ 0O resulting 
from the Hankel technique would be less than or equal to that resulting from the 


balancing technique. 
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It is of interest to compare these two results for E^: for the Hankel 
technique [Ref. 14], E <*, = tr n and as proved above for the balancing technique, 
j = 2 a n . Thus the balancing technique yields a solution which has an error 
twice that of the Hankel technique which has an extra degree of freedom (i.e. the 
D term). 

In the more general case of any r < n and MIMO but with unity weightings 
Glover [Ref. 15] states an algorithm for model reduction based on the Hankel 
technique for which 

£oc ^ tr j^2 j 

where again it involves a finite D term. 

It should also be noted that if the D terms for these techniques are set to 
zero and the reduced order model is otherwise unchanged, the following bound 
holds in this case 

£» < ■ 

Thus short of knowing the optimal solution for the unweighted, zero D case, the 
internally balanced realization model reduction technique will be regarded as an 
attractive technique for performing model reduction. The computations for the 
Hankel technique start with a balanced realization, hence, obtaining the reduced 
order model from the internally balanced realization is computationally less 
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ex pensive and satisfies the same error bound as the Hank el technique for the 
same assumptions on the reduced order model. 

Frequency Weighted Balanced Realizations 

We have seen that the internally balanced realization is an attractive model 
reduction technique for unity weightings. An extension to the balancing 
technique to include weightings will now be developed. 

This extension is motivated by the geometric interpretation of the internally 
balanced realization given in Chapter 13. The n th order, minimal, asymptotically 
stable system to be balanced with respect to the asymptotically stable input and 
output weightings and W 0 (s) satisfies the given state equations 

x = Fx -f Gu zeR" ueR m 
y — Hx yeRP 


With reference to the block diagram shown in Fig. 131.7, ask the following 



Fig. HI.7 Block Diagram For Frequency Weighted Balanced 


Realization Discussion 
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two (dual) questions. What set of points in the x-state space could be part of the 
zero initial condition response for some weighted input, fi(t) such that 
IIMOII 2 ^ 1 and zero initial conditions for the states of the input weighting? 
and what set of points in the x-state space as initial conditions could produce a 
weighted output, rj(t) such that ||»7(< )|| 2 < 1 with zero input and zero initial 
conditions for the states of the output weighting? 

These questions are again answered by appealing to the grammian-ellipsoid 
theorem of Chapter D. First assume asymptotically stable state space 
realizations for the input and output weightings are given, i.e. 

W,(s) = H^sl-F,)-' G t + D, 

M» = HM-F,)-' G.+ D, . 

For the first question, let F(t) — H i e F,t G ,• where 


F GH- 


GD- 

0 Fi 

w — 

G { 


Hi = [/ 0 ] , 

then the answer is the controllable set, S e . 


M P £ MW 

For the second question, let F(t) = H c e ‘ G 0 where 


F a = 


F 0 
G„H F„ 


* _ m 

- Loj 


H 0 = \D 0 H H 0 ) 


The grammians corresponding to these questions are given by 
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_ w 

U A J o H, e Flt G { G? e F ‘ l tij dt 

/ oo 

0 aj e F -‘ nJn c e F ‘‘ G,dl 

and thus can also be computed from the cascaded grammians 


U A 

1 

f o e F,t GiG? e F '\ 

Y A 

< 

f o e F ° l HjH 0 e F ° l 


which can be computed from the Lyapunov equations 


F { U + UF? + G { Gj = 0 
FjY + YF. + HjH 0 = 0 . 

Partition U and Y such that their upper left block is n X n as 


U = 


Un 

U 2 i 


ul 

U 22 


Y = 


I'll ^12 

Yn Y 22 


then 

U = H { UH? = U n 
Y = GjYG „ = Y n 


(HI. la) 
(m.ib) 


Thus again, the two sets are in general different ellipsoids with n semi-axes. 
The lengths and directions of the weighted controllability ellipsoid semi- axes are 
the square roots of the eigenvalues and eigenvectors of the weighted 
controllability grammian U, respectively. The lengths and directions of the 
weighted observability ellipsoid semi-axes are the reciprocals of the square roots 



- 98 - 


of the eigenvalues and the eigenvectors of the weighted observability grammian, 
Y respectively. 

The next step is to obtain a realization of the system such that the two 
ellipsoid axes are lined up (balanced). To obtain a new realization with this 
property a change of variables with a similarity transformation will be used i.e. 
x = Tz. 

A possibly subtle point should be noted here: this is not the same as 
internally balancing the cascaded systems. That is, the n X n grammians, U 
and Y are balanced not U and Y which in general will not even have the same 
dimensions. 

Algorithm For Computing the Frequency Weighted Balanced 
Realization 

The objective of the balancing algorithm is to obtain a realization 
{A, B, C}„ of G(s) = H(sI-F)~ 1 G = C{ S I-A)~ l B such that U and Y, the 
weighted grammians are equal and diagonal. Towards this end, consider how the 
original grammians change due to a similarity transformation on the full order 
model state variables. 

Starting with III. la in partitioned form together with U = U n we have 
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F ghA u U&] U u! 1 ] l F t o' 
0 F { U n C/ 2 2 C/ 21 ^22 H?G t F? 



Introducing the transformation, T, we have 



which can be rewritten 



Similarly III. lb can be rewritten 



where 
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A A T~ l FT B A T-'G 

C J HT 

^ ^21] A T-'UT~ t T~ l Uli 
U 21 U 22 \ ~ [U 21 T~ t U 22 

}V ^ 12 ] A t t yt T t Y 12 
y f 2 y 22 j = y[ 2 t y 22 ‘ 

Now the objective is to find T such that 

0 = T~ x \JT~ t = E 
Y = T t YT = E 

where E = diag {a 1 , <x 2 , fr n } 

with > <r 2 > • • • > (T n 

The similarity transformation, T which accomplishes the objective is the 
eigenvector matrix obtained from the eigenvector decomposition 

UY = TAT' 1 

as for the internally balanced case discussed in Chapter II. Noting that 

U > 0 => U > 0 

and 

F > 0 => Y > 0 , 

the result in Appendix C shows that A > 0, thus letting 

<r { A y/Xj i = 1, 2, n 
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where A — diag {X, } accomplishes the objective (assuming without loss of 
generality that the X,’s have been ordered such that Xj > X 2 > • • • > X n ). 

The frequency weighted balanced realization {A,B,C} is also unique to 
within a change in sign of a state variable when the singular values of the 
balanced grammian are distinct (i.e. <r,- ■=£ <7y for t j). The proof of this fact 
for the weighted case is identical to that of the internally balanced case [Ref. 11] 
and hence is omitted. 

The algorithm for transforming a given asymptotically stable, minimal 
realization {F, G,H} n into a frequency weighted balanced realization, given state 
space realizations for the transfer functions of the asymptotically stable input 
and output weightings is summarized in Table DLL 

Thus it has been shown that there exists an essentially unique realization of 
the system obtained by a change of variables with a similarity transformation, T, 
such that the axes of the ellipsoids are the same (balanced) for the new state 
variables (say, z where x = Tz). For model reduction this balancing is the key 
idea. It provides a basis for the n dimensional x -state space where the direction 
of a given basis vector is as controllable with the weighted input as observable 
with the weighted output in a well defined sense. The lengths of the ellipsoid’s 
semi- axes provide a scalar measure, <r,- of how important a given basis vector 
direction is with respect to the weighted controllability and observability. 

The reduced order model is then obtained by neglecting the states which are 
weakly controllable/observable with the weighted input/output. In mathematical 
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Table ELI Frequency Weighted Balanced Realization Algorithm 

Given: F , G , H with R c fxjFjj < 0 

F it Gi t H it Di with R e jxfF,]] < 0 

F g , G 0 , H 0 , D 0 with R t [\[F 0 ]j < 0 

where G(s) = H{sI-F)~ l G 

W { (s) = H { (sI-F { )-'G { + D { 

W 0 (S) = H 0 (SI-F 0 )-' Gg+Dg 

Solve for U and Y from 



Solve for eigenvalues and eigenvectors of UY 

i.e. UY = TAT- 1 , A = diag {X, } , X, > X 2 > ••• > X n 
£ = diag{cr,} o { A y/kj i = 1, n 

A = T~ X FT B = T~ l G 

B = HT 
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terms this involves partitioning of the frequency weighted balanced realization, 

{A,B,C}: 



C = [C, C 2 ] 


The rth order reduced order model is then given by £,.(«) = C x {sl-A 
where the partitioning was done such that A u is r X r and <r f+ x < a r . 

Stability of the Frequency Weighted Reduced Order Model 

Consider the two special cases: 

Case 1: Wj(s ) == / and W # («) ^ I 
Case 2: U^s) ^ / and = / . 

By examining the partitioned form of the grammian equations we have: 

For Case 1: the reduced order model satisfies: 

-AjiEj + = 0 

For Case 2: the reduced order model satisfies: 

4- + C\C X = 0 

Since the full order model is assumed to be minimal, the balanced grammian 


is positive definite i.e.: 
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Case 1: AE 4- E A T + BB T — 0 , E > 0 
Case 2: A r E 4- EA + C T C = 0 , E > 0 

Now since Ej is just a leading partition of a positive definite matrix, it too is 
positive definite, thus for either case Ej > 0. 

It is well known, [Ref, 36] that the partitioned grammian equations (i.e. 
Lyapunov equations) together with Ej > 0 implies that /?e [X[>1 n]] < 0, 
furthermore if for 

Case 1: (An, B x ) is controllable then i?e [X[^4 n]] < 0 
Case 2: (A n , Cj) is observable then i?e[X[Anj] < 0 

Thus for either case, the reduced order model is guaranteed to be stable and 
generically (assuming controllability or observability of the reduced order model) 
the reduced order model is asymptotically stable. 

For the general case of non-unity weightings it is not known whether the 
reduced order model will be stable or not. As will be seen in Chapter IV, the 
general case will not be required for the remainder of this thesis anyway. 

Frequency Weighted Balanced Realization Error Bound 

The manipulations leading to the internally balanced (unity weightings) 
realization model reduction error bound were carried out for the frequency 
weighted (non-unity weightings) balanced realization. These manipulations are 
contained in Appendix D and result in 



- 105- 


£oc M II W # («)[G(.)-Gf r (.)]Wi(.)|| 00 < 2(1 + a)tr[Z 2 ) 

where unfortunately a simple bound for the positive quantity, a could not be 
found. The quantity, a is related to the H ^ norm of the transmission from the 
weighted input to the neglected states and/or from the neglected states to the 
weighted output. The premise of the model reduction is that these transmissions 
are small, hence a conjecture is that: a < 1 when E <*, < 1. 

Short of knowing something in closed form about a it is probably best to 
obtain the error, E ^ for the model reduction by direction computation. 

Model Reduction of Full Order Models With Unstable Poles 

As we have seen, both the internally and frequency weighted balanced 
realization model reduction techniques require that the full order model be 
asymptotically stable. In many cases model reduction is necessary when the full 
order model contains poles in the closed right half plane. Some work has been 
done in formally extending the internally balancing technique to include unstable 
poles [Ref. 42] however these results were deemed to be beyond the scope of this 
thesis. Rather the following scheme was used and is recommended. 

Given the transfer function of the full order model, G(s) it is always 
possible to perform the decomposition (partial fraction expansion) 

C( S ) = C,(») + C„(») 
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where G t (s) contains all the asymptotically stable poles of G(s ) and G ue (s) 
contains all the unstable poles of (7(s). 

Now since for control design, the reduced order model must have the same 
number of unstable poles as the full order model anyway, the reduced order 
model can be obtained by reducing the order of only the stable part, G g (s) and 
leaving the unstable part as is. That is apply the balancing model reduction 
technique to obtain G' #r (« ), a reduced order model of (?,(«), then the reduced 
order model, G r (s) of G(s) is given by 

G r (s) = G tr (s) + (?„,($) 

The error using this technique will be that due to the order reduction of the 
stable part only, which can be seen by 

J IlinWI «?(*) - <?,{«)] H'-(.)||„ 

= iitnwiG,w+ <?„(.)-<;„(*)- w;-(») il 
= II iv.(«) [«?,(.) -(?„(.)] w,(.)||„ . 

The unanswered question is: could a smaller E ^ be obtained by allowing 
the locations of the unstable poles to change but remain in the closed right half 
plane? A partial answer to the question is obvious for poles on the ju axis. In 
this case, the reduced order model must have poles at exactly the same locations 
on the ju> axis as the full order model or E ^ will not be bounded. A complete 
answer to the question was not obtained. 
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Model Reduction Concluding Remarks 

This chapter has developed a unique model reduction technique which is an 
attractive solution to the model reduction for control system design problem. 
The key idea is that of frequency dependent weightings. The choice of 
weightings will be discussed in the next chapter. 



- 108 - 


IV. CHOICE OF MODEL REDUCTION WEIGHTINGS 
FOR CONTROL DESIGN 

This chapter will develop two approaches for choosing model reduction 
weightings. The first approach assumes that it is the model of the compensator 
for which a reduced order model is desired and the second assumes that it is the 
model of the plant for which a reduced order model is desired. Both of these 
approaches result in a reduced order controller being designed which provides a 
stable closed loop system for the full order model of the plant. 

Motivation For The Choice of Weightings 

The two approaches for designing reduced order controllers are motivated by 
the stability robustness theorem of Chapter II which is repeated here for 
convenience. Let P(s) be a multivariable transfer function which is stable under 
unity feedback as shown in Fig. IV.l with 



Fig. IV.l Perturbed Closed Loop System 
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/ 


A(s) — 0. The result is that the perturbed closed loop system remains stable for 
all stable perturbations, A(s) which satisfy either 

||A( S )[/+ />(,))-' IUC1 
or 

in/ + /(sir 1 a(»)il<i 

and are referred to as input or output uncertainty tests respectively. This 
theorem is used to derive appropriate weightings for model reduction as part of 
control system design. 

Compensator Order Reduction 

Assume a compensator transfer function, K(s) has been designed for the 
transfer function of the system to be controlled, G(s) by some technique to meet 
performance and stability robustness specifications of the closed loop system. It 
is assumed that K(s) has an order large enough to warrant reducing its order to 
give A' f .(s), a reduced order approximation of K(s). Motivated by the stability 
robustness theorem, consider the block diagram in Fig. IV.2. For the nominal 



Fig. IV.2 Compensator Order Reduction As a Perturbation 
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loop (i.e. K r (s ) - K(s) — 0) the closed loop system is the result of the full order 
design. On the other hand the perturbed system corresponds to the closed loop 
system with the reduced order compensator and the full order model (this is the 
closed loop system for which stability is required). The stability of the closed 
loop system with the full order model and the reduced order compensator is 
guaranteed by 

E„ = ||H'.(»)[*(») -*,(•)! W i (.)|| 00 <l 

where to apply the theorem the uncertainty due to the model reduction is 
represented arbitrarily at the input or the output: 

for output uncertainty representation 
W,{a) ='[/+ <?(«)/£(* J]- 1 G[a) and W-(«) = / 

for input uncertainty representation 
w;.(«) = G(s)[7 + it:(s)G'(s ))- 1 and W,{a) = I . 

Actually part of the uncertainty can be put in both places but this generality was 
not used in this thesis. 

The reduced order controller, K r (s) is now determined by finding the 
reduced order model of K{s) from the frequency weighted balanced realization of 
K(s) with the weightings for either the input or output uncertainty 
representation. The algorithm for this approach is summarized in Table IV. 1. 

It can be readily verified that the two non-unity weightings are the same i.e. 
£($)[/ -I- K(s )G(s )] _1 = [I + G'(s)A^(5)] -1 G'(s). Thus the difference in the two 
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Table IV. 1 Compensator Order Reduction Algorithm 

• Design the full order compensator transfer function, K(s) to meet the design 
objectives for the plant transfer function, G(s). 

• Then choose input and output weightings by arbitrarily representing the 
model reduction uncertainty at either the input or the output: 



• Obtain the frequency weighted reduced order model, K r {s) of K{s) for 
iy-(s) and IT^s) such that E ^ < 1. Then the closed loop system will be 
guaranteed to be stable with the reduced order compensator K r (s). 
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uncertainty representations results in the weighting being used either as an input 
or an output weighting for model reduction of A(s). A realization of this 
weighting is given in Table IV.2. 

Note that for either case the input and output weightings have poles which 
are those of the closed loop system with the full order compensator and the full 
order model. This closed loop system will always be stable, hence, the weightings 
will always be stable. 

Some examples of the application of this approach are contained in Chapter 
V. This approach is quite straightforward and most of the remainder of this 
chapter will deal with the second approach. Advantages and disadvantages of 
the two approaches will be discussed and compared at the end of this chapter. 

Plant Order Reduction 

With this approach, the error due to the reduced order model of C(s) is 
associated with the perturbation term of the stability robustness theorem. The 
block diagram shown in Fig. IV.3 results. When the 
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Table IV.2 Realization of <?(«)[/ + K{s)G{s)]~ 1 = [I + G(s)A:(s)]- 1 G(s) 


let G(«) = C{8l-A)~ X B 

and H(a) = H{aI-F)~ l G 

then a realization of 

W{ s ) J (?(«)[/+ K{a)G{s j]' 1 = [/+ G(a)/f(s)] _1 G(s) 
is given by 

W(s) = CJsI-AJ-'B, 
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K r (S) 


Fig. IV.3 Plant Order Reduction As a Perturbation 

perturbation term, (?(s) - (7 r (s) is zero, the block diagram represents the result 
of a reduced order compensator design for the reduced order model, G r (s). That 
is, assume is designed to satisfy the design objectives for the reduced order 

model, G>(s). On the other hand, the perturbed system corresponds to the 
closed loop system with the reduced order compensator and the full order model 
(this is the closed loop system for which stability is required). 

The stability of the closed loop system with the full order model and the 
reduced order compensator is guaranteed by 

£00 = I! !*.(«)[<?(«) -G r (s)} w;-(«) IL < 1 

where to apply the theorem the uncertainty due to the model reduction is 
represented arbitrarily at the input or the output: 

for output uncertainty representation 
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W-(«) = AT r («) [J + <?(«)*:,(«) ] _1 and W,{s) = I 
for input uncertainty representation 

W { (s) = / and W 0 (s) = [I + K r (s)G r (s)}~ 1 K r (s) . 

Ag;ain, it can be readily verified that the two non-unity weightings are the 
same i.e. K r (s)[I 4- G' r (s)/C r (s) ] _1 = [/ + AT r («)G' r («) ] -1 /sT r (s). Thus again, 
the difference in the two uncertainty representations results in the weighting 
being used either as an input or an output weighting for model reduction of 
G(s). 

Note that for either case, the input and output weightings have poles which 
are those of the closed loop system with the reduced order compensator and the 
reduced order model. This closed loop system will always be stable, hence, the 
weightings will always be stable. 

Unfortunately the weightings for reducing the order of G(s) are not known a 
priori. The weightings depend on both the compensator and the reduced order 
model which are not known before the model reduction. 

Elimination of Compensator Dependence With Advanced Loop Shaping 

Advanced loop shaping, which was discussed in Chapter II is an LQG based 
design procedure, which for a large class of design problems, allows the designer 
to a priori specify a desired loop shape for one of the input or output loop 
transfer functions, K(s)G(s) or G(s)K(s) respectively. That is before actually 
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computing K(s), the input or output loop transfer function is known. 

Assuming that advanced loop shaping will be used to perform the reduced 
order controller design, it is of interest to make use of the a priori known loop 
transfer function to determine the a priori unknown model reduction weightings. 
Towards this end let 

#.(*) M c,(* )*,(») [/ + G r (.)K,(«)|- 1 

and 

Hii») J [/+ K,(s)GAs) ]-> K,(t)G r (o) 

then for input uncertainty: 

wjs) = \r + km 

= HM G r -‘(.) 

and for output uncertainty: 

MvU) = K t U)lI+ G r {,)KAs)}-' 

= C r -'(s) //„(,) . 

Note that H { (s) and H 0 (s) depend only on the loop transfer functions 
K r (s) G r {s) and G r (s) K r (s) respectively, which are known a priori (with 
advanced loop shaping). Thus either //,•(«) or H 0 (s) will also be known a priori. 


Thus for either the input or output uncertainty representation, the dependence of 
the weightings on the reduced order compensator has been eliminated. 
Unfortunately the weightings still depend on the reduced order model which is 
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not known a priori. 

Plant Order Reduction For Control System Design Problem Statement 

To simplify the further discussion assume the model reduction uncertainty is 
represented at the output. The dual results will be summarized later. Then the 
problem can be stated: 

Given: G'(s), H 0 (s), r 

Find: a reduced order model, G r (s) of (7(s) 

such that E x = || [G(s) - C r (,)l )#.(•) Hoc < 1 ■ 

Well the model reduction technique of Chapter HI addresses problems like 
this, however, it assumes that the weightings are known. In this case the 
weighting G f -1 (s) H 0 (s) depends on the result of the model reduction. 

Parameterization of the Solution 

Without loss of generality assume an n tk order, output normal minimal 
realization, {F, G, H} n of the stable portion, ^(a) of the plant, <?(«) is given. 
That is 

F + F T + H t H = 0 

where C,(s) = H(sI-F)~ l G and G(s) — G>(s) + (^(s) (^(s) has poles in 



- 118 - 


the open left half plane and G ue (s) has poles in the closed right half plane). Such 
a realization (or the input normal, dual) can be readily obtained by performing 
the steps in Table IV.3. 

For the sake of further discussion assume the input weighting, ^-(s) which 
solves the problem is known. Then let U be the weighted controllability 
grammian for the stable portion of the plant, (7,(5) and the weighting, VF,(s). 
This grammian, U can be computed as shown in Table ELI with a Lyapunov 
equation, however, note that by using Parseval’s Theorem, U is also given by 
00 

U _ _L J {jui-F)-' GW i {ju)Wf I {ju)G T {-juI-F T )- 1 du . 

2n -00 

This is introduced only for the sake of a more compact notation than that of 
Table III. 1 and computations are probably easier with the Lyapunov equation of 
Table IH. 1 . 

The next step of the model reduction is to balance the weighted 
controllability grammian with the identity observability grammian (i.e. F , G, H 
is an output normal realization). This balancing is accomplished by finding the 
eigenvalue decomposition of U i.e. 

U = WXW T 

where A == diag {X, }, \ x > X 2 > • • • > X n > 0, X,eR \fi and W T W = I. 
Then the transformation to a frequency weighted balanced realization is obtained 


by letting 
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Table IV. 3 Input and Output Normal Realizations 

Assume an arbitrary, minimal, asymptotically stable realization {A ,B , C} 
is given. 


Input Normal 

Output Normal 

solve for U from 

solve for Y from 

AU + UA T + BB T = 0 

A T Y + YA + C T C = 0 

perform the eigenvalue decomposition 

perform the eigenvalue decomposition 

U = VAV T , A = • diag {X..} > 0 

Y = VAV T , A = diag {X-} > 0 

X.eS Vi' and V T V = I 

X^Jt Vi and V T V = I 

let T = VA 1/2 

let T = VA" 1/2 

then 

then 

F T + F + GG T = 0 

F + F T + H T H = 0 


where F = T -1 AT G = T -1 B 


H 


CT 
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T = WA 1 / 4 . 

Note 

T-'UT- T = A" 1 / 4 W T WAW t WA - 1 ! 4 = A 1 / 2 
T t IT = A 1 / 4 W t IWA 1 ! a = A 1 / 2 

and thus = diag {X,- 1 / 2 } is the balanced grammian. The frequency weighted 
realization {A, B, C) of the stable portion G,(s) of the plant, G'(s) is given by 

A = T~ X FT B = T~ l G 
C — HT . 


Note that this frequency weighted balanced realization is equivalent to that 
which would have been obtained by applying the alternate algorithm for 
computing it, given in Table ELI. 

Now then assuming the order, r of the reduced order model is given, the 
reduced order model is obtained by partitioning 

A __ ^ll ^ 12 

a 21 a 22 

C = [ C x C 2 ] where A n is r X r 

The realization of the reduced order model, (7, r (s) of the stable portion, <7,(s) of 
the plant, (7(s) is given by 

G„(*) = C x {sl-A n )- X B x . 


B = 


Bx 

B 2 


Note that: 
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GJs) = B x 

= HW X E// 2 (s/ - Sf 1 / 2 W{FW X Z }/ 2 j 1 Ef 1 / 2 VFfG 

= Wj(s/ - WjFW x Y l WjG 

[E, 0 

where W — [W x W 2 ] and E = q ^ which are both partitioned 

2 

compatibly with A, B, C (i.e. W x has r columns and E 1 is r X r). Thus the 
reduced order model, <7, r (s) does not depend explicitly on E and hence only the 
eigenvectors of U corresponding to the r largest eigenvalues of U are necessary 
to compute the rth order reduced order model. 

The reduced order model of the plant is given by 

c,(«) = <?„(»)+ <?„(.) . 

Since the input weighting from the beginning of this discussion was assumed to 
be the correct weighting it must also be given by 

HK») - . 

The previous, hypothetical, discussion is summarized schematically in Fig. 
IV.4. Note that if either U, W, W x , G ir (a), <? r («) or were known the 

problem would be solved. Note also that the diagram also represents functions of 
these quantities onto themselves. For example, take any positive definite, n X n 
matrix U in as an input to the upper left hand block in Fig. rv.4 and regard the 
output of the lower left hand block, U 01lt in Fig. IV.4 as the output of the 
function, / [U] i.e. 


4 . 




Figure IV. 4 Schematic of the Inter-relationships for Plant Order Reduction Solution 
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That is f(U) is defined algorithmically as (assuming F , G , H, G vt (s), H 0 (s) and 
r are given): 

U = VF = [W 1 W 2 ) 

where W\ has r columns 

A n = ITfFtTj B x = Cj = W, 

<?,(»> = C, (»/->! i C„(.) 

»f(») - G r -‘(*)H.(*) 

00 

/(tf) = -f- / (y«/-F)-»C7W;.(;w)W;fli[yw) C r (-/w/-F r ) - 1 do; . 

2?r -00 

This series of steps is rather involved but nevertheless does define the function, 

f(U). 


Statement of the Plant Order Reduction Problem as a Fixed Point 
Problem 

The plant order reduction problem can now be stated: 

Given: F , G, H, G u ,{s), //„(«), r 
Find: U such that U — / ( t/) 

Thus the solution to the plant order reduction problem is equivalent to finding 
the stationary or fixed point of the function f{U) which maps n X n positive 
definite matrices onto n X n positive definite matrices. 
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Similar functions can be obtained for W , W j, G, r (s), G r (s) or Wf(«) as 
well. The functions involving the transfer functions: G, r (s), G r (s) or W[-(s) are 
not of interest, since specifying a transfer function typically involves specifying 
realizations which are not unique. The function involving is of interest, since 
it involves the fewest parameters as unknowns. 

The existence and uniqueness of stationary points for these functions is an 
open question. These issues were deemed to be beyond the scope of this thesis 
and the following solution technique was used successfully. 

Successive Approximation Solution 

The successive approximation method to obtain the solution of the equation 
x = f (x) is to guess an initial x say, x l and let 

*.-+i = /{*,) * = 1. 2, 

then under appropriate conditions 

x — lim X{ 

I — *00 

The appropriate conditions [Ref. 43] are that the solution exists and that the 
function, f (x) is a contraction i.e. the first derivative of the function as a matrix 
should have eigenvalues with magnitudes strictly less than unity in a 
neighborhood of the solution. 

This technique resulted in satisfactory solutions for all the examples tested 
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for this thesis (see Chapter V). However, the procedure did not always converge 
to a unique point. In several examples the procedure converged to a limiting set 
rather than a limiting point. That is a limit cycle was reached: 

/(*,+ ») = /(*,') ■•»! 

where k was the number of elements in the limiting set. In these cases, the value 
of was compared for the k elements of the limiting set and the one 
corresponding to the smallest E ^ was taken as an acceptable solution. 

Other Solution Techniques 

The problem of solving x == f (x) is also theoretically amenable to more 
advanced iterative solution procedures involving derivative(s) of f(x). These 
techniques were not used, however, since the computation of the gradient of this 
function was deemed to be beyond the scope of this thesis. 

Obtaining a Realization of G~ ][ (s)H{s) 

In order to carry out the successive approximation solution procedure, it is 
necessary to obtain a realization of G~ l (s)H(s) where realizations of the square 
tn X m transfer functions, G(s) and H(s) are given. Although a realization for 
G _1 (s) does not exist in general a realization for G~ l (s)H(s) often does exist. 
The existence depends on the pole-zero excess of both G(s) and H(s). Two 
special cases were used in this thesis and they will now be discussed. 
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G(s) With First Order Rolloff 

Assume {A, B, C} is a state space realization of G(s), i.e. 
(7(s) = C(sI-A)~ l B . Then first order rolloff is defined to mean that CB has 
full rank. The connection between this and a rolloff slope of negative one is 
obtained from the high frequency expansion of G(s): 

G(s) = C(sl-A )~ l B 

= C[s{I - - A) ]-* B 
s 

= C— (/ + —A + —■ A 2 + ...) B 

S S S 2 

= - CB + -V CAB + CA 2 B + • • • 

® s 2 s 3 

Thus if CB has full rank then all the singular values of G(ju) will have a slope 
of negative one for large w when plotted on a log-log scale versus ui. 

For a realization of G~ 1 (s)H(s) to exist in the full rank CB case, it is 
necessary that H(s) have a rolloff slope of negative one or less. This will always 
be true if H(s) has a state space triple (as opposed to quadruple) realization, say, 
{F , G, H } (i.e. zero feed thru or D term). 

It can be verified by direct computation that for any aeR : 

W( S ) = G-'(s)H(s) = CJsI-AJ-i B w + D w 
G{s) = C(sI-A)~ l B, CB full rank 
H{s ) = H(sI-F)~ l G 


where 
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A 

A + B{CB)~ l C(aI-A) -B{CB)~ l H{cd-F) 

B A 

B{CB)~'HG 


0 F 

G 


C. J \cBy'C(al-A) -(CB)->ff(a/-F)j D. ± [(CBy'Ho) . 

Thus {A„,, B w , C w , D w } is a realization of VK(s) = G'" 1 (s)//(s). However, 
this realization of W(s) is not minimal, in fact, it has m uncontrollable poles at 
s — a (these poles had to be introduced artificially to produce the realization for 
W(s)). The fact that these poles are uncontrollable can be verified by showing 
that the left eigenvectors of corresponding to the eigenvalues of A w at s — a 
are orthogonal to B w [Ref. 34]. These eigenvectors are given by 

V T = [c -//] 

and it can be verified that V T A w — aV T and V T B w — 0. 

Thus a non-minimal realization of IT(s) has been obtained. There are 
several ways of eliminating these m uncontrollable poles at 8 = a. in fact, a 
method using internally balanced realizations for obtaining a minimal realization 
was discussed in Chapter II. 

The realization of the dual, W(s) = i/(5)G _1 (s) and the realization of 
Mis) == G~ 1 (s)H(s) are summarized in Table IV.4. 


C(s) With Second Order Rolloff 

Assume {A, B, C} is a state space realization of G'(s), i.e. 



W(s) 

Note: 


W(s) 


Note: 


Table IV. 4 Realizations of G~^(s)H(s) and H(s)G~^(s) For G(s) With First Order Rolloff 
G(s) = C(sI-A) -1 B with CB full rank, H(s) = H(sI-F) -1 G, ae!R 


z\s)H(s) 

C (sI-A ) _1 B +D 
w v w' w w 



A + B(CB)" 1 C(aI-A) -B(CB) _1 H(aI-F) 


B(CB) -1 HG 

A £ 


B £ 


w 


w 



0 F 


G 


I 


V T A = aV T 
w 


V B w ■ 0 


w 

,T 


H ( s ) G~ 1 (s) 

C (sI-A ) -1 B +D 
w v w' w w 


w 


w 


^ [ (CB) _1 C(aI-A) 

£ [ C 


A 


A + (aI-A)B(CB) -1 C 
-(aI-F)G(CB) _1 C 


= [ HG(CB) _1 C 


~(CB) -1 H(aI-F) ] D w ^ [(CB) _1 HG] 

-H ] 

(al-A)B(CB) -1 
-(al-F)G(CB)' 1 

H] D w = [ HG(CB) -1 ] 


0 

F 


A, ,V = aV 
w 

C w V = 0 


00 
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C(s) = C(sl-A )~ X B . Then second order rolloff is defined to mean that CB = 0 
but CAB has full rank. Again the connection between this and a rolloff slope of 
negative two is obtained from the high frequency expansion of G(s): 

G(s) = - CB + -- CAB + ™ CA 2 B + • • • 

9 S 2 fi 3 

Thus if CB = 0 and CAB has full rank then all the singular values of G(ju) will 
have a slope of negative two for large oj when plotted on a log-log scale versus u. 

For a realization of C -1 (s)//(s ) to exist in the CB = 0 and full rank CAB 
case it is necessary that H(s) have a rolloff slope of negative two or less. Thus if 
{F, G, H } is a realization of H[s), i.e. H(s) = H(sI-F)~ l G , it is necessary that 
HG = 0. 

It can be verified by direct computation that for any a /3 and a* /?cR : 
W(s) = <?-*(«)#(«) = C W ( S I-AJ- 1 B u + D w 

where C(s) = C{sI-A)~ x B, CB — 0, CAB full rank 
H{s) = HisI-Fj-'G, HG=0 


n A 

W = 

C(al-A) [fil-A ) 

D, A 

CAB 

Hk & 

H{al-F) (/3/-F) 

Dk & 

HFG 

A A 

A-BD~ l C 9 BD~ [ H h 
0 F 

B . A 

BD-'D h 

G 

Ho 

[D;'C, Df'H h } 

D A 

w = 

] 
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Thus {^4^ , B w , C w , D w } is a realization of W(s) = G~ 1 (s)H(s). However, 
this realization of VK(s) is also not minimal, in fact, it has 2m uncontrollable 
poles: m at s — a and m at s = fi (these poles had to be introduced artificially 
to produce the realization for W(s)). The fact that these poles are uncontrollable 
can be verified by showing that the left eigenvectors of A w corresponding to the 
eigenvalues of A w at « = a and s = ft are orthogonal to B w [Ref. 34). These 
eigenvectors are given by 


f/r = 


r C(fil-A) -HI/3I-F)] 
[C(al-A) -H(al-F ) J 


and it can be verified that 


V T A, 


\al 0l v t 

10 fill v 


and 


v t b w = 0 . 

The realization of the dual, IT(s) = H(s)G~ 1 (s) and the realization of 
H 7 ( s ) = ■G~ 1 (s)H(8) are summarized in Table IV.5. 

Realizations of = G~ l (s )H(s) and = H(s)G~ l (s) have been 

given for two special cases which include a large class of transfer functions, G(s). 
In those cases not included it is still possible to obtain a realization of IT(s) by 
modifying the procedure used to obtain the results in Tables IV.4 and 5. 




i A _B O"" 1 C 

= C (sI-A )" B +D ft A 99 
w' w' w w A = i 

G.D^C 

n g 

al 0~J L 

= V 

0 BI 


Roll off 


a f Be ft 
D h £ HFG 


D. = HFG 
h 

(BI-A)B (al-A)B 


-(BI-F)G -(BI-F)G 
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Obtaining a Stable Weighting When G'(s) is Non-minimum Phase 

Often the reduced order model, G r (s) resulting from a balanced realization 
will be non-minimum phase (typically at high frequency) and in these cases the 
weighting, )H(s ) or W(s) = H(s)G~ l (s) will have unstable poles. 

This is because the transmission zeros of G r {s) will become poles of W(s). In 
this case it will be necessary to modify the transfer function of the weighting to 
make it stable. 

Consider the definition of the model reduction error criterion, E ^ 

•e» A. IKw [c(..)-c r (.)].« , i(*)IU 

then by definition of the H ^ norm 

E oo = snp a[w t (ju) j G(ju)- G r (ju ) j W,(/w)j 
then by definition of the maximum singular value 

e„ = supxig [wjju) [c(M - <?,(/«)] 

[C'l/w) - W?(ju) 

then since X max [v4.5] X max [fiA] \/A,B , 
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e , o = tupx'fi [[c(y«)-e f (y u )]H'.(y w )HV , (/«) 

[c? fl (y w )-G r iI (; W )] Mflyujw.yw) . 

Thus E oq depends only on the products W^*(/w) W^(/w) and Wf(ju) W 0 (ju>) and 
not on the weightings by themselves. 

Spectral Factorization 

Consider the product W(ju)W B (ju>) = IT(s) W^(-s) | ,_y w when VF(s) has 
some poles in the open right half plane. To facilitate the following disscussion let 
VT + (s) A IT(s) to denote the fact that VF + (s) has poles with positive real parts. 
It is well known [Ref. 44] that a transfer function VF_(s) always exists with poles 
in the open left half plane such that 

IT_(s)IT_ r (-s) = W + {s)Wl{-s) . 

This is called spectral factorization in the literature. Since E <*, depends only on 
products like above, the stable spectral factor can always be used in place of the 
unstable factor without changing E^. 

Computation of the Stable Spectral F actor 

It is well known [Refs. 1-4] that the minimum input (i.e. minimum of 
J — f™ u T (t)u(t)dt) control (i.e. u = -K e x) required to stabilize the unstable, 



- 134 - 


controllable system x — Ax + Bu just reflects the unstable poles of the system 
across the imaginary axis. That is the eigenvalues of .4 - BK C are the stable 
eigenvalues of .4 plus the stable reflections of the unstable eigenvalues of .4. The 
gain K e is computed by solving for P > 0 such that 

A T P + PA - PBB t P = 0 

and then K c — B T P. 

It is not so well known that these facts can be used to solve the spectral 
factorization problem in the multivariable case (it’s trivial in the SISO case). It 
has been shown in Reference 5 and it can be verified by direct computation that 
the given transfer function IT + (s) A C(sI~A)~ 1 B can be factored 

W + (s) = lT(44(s) 

where 

IT (s) J C(sl-A+ BK c )- l B 
A{s) J 1+ K c (sI-A)~ l B 

and where 

.4 T P + PA - PBB t P = 0 , P > 0 

and 

K c = B t P . 

It can also be verified that A(s) is a multivariable all pass, i.e. 
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A r (-s)A(s ) = A(s)A t (s) = I 

Thus we have 

= VK_( s )^(s)A 7 'H)^_ r (-s) 

= VK_(s)lT_ r (-s) 

which accomplishes the objective. 

These results have a dual for the other spectral factorization problem 
W’J (-s ) W + (s ) = (— s ) W_(s ) and both cases are summarized in Table IV.6. 

The results for plant order reduction are summarized in Table IV.7 for both 
representations of the model reduction uncertainty. 

Advantages and Disadvantages of the Two Approaches 
(Compensator vs. Plant Order Reduction) 

Compensator order reduction requires that the high order compensator be 
designed for the high order plant. This is clearly undesirable and is in fact the 
motivation for searching for reduced . order models. The advantage of this 
method is that implementing the algorithm for compensator order reduction is 
straightforward (i.e. the model reduction weightings are known a priori). Often a 
low order approximation to a given compensator is desired and thus this 
attractive, formal method for obtaining the approximation is a very useful result 


of this thesis. 
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Tabl e IV. 6 Soectral Factorization 


W + (s) = 

C(sI-A) _1 B 

W + (s)wJ(-s) = W_(s)W^(-s) 

wJ(-s)W + (s) = W^(-s)W_(s) 

Solve for P > 0 from 

Solve for P > 0 from 

A T P + PA - PBB T P = 0 

AP + PA T - PC T CP = 0 

let K = B T P 
c 

let K f = PC T 

then W + (s) = W_(s)A(s) 

W + (s) = A(s)W_(s) 

where 

where 

W_(s) = C(sI-A+BK c ) -1 B 

W_(s) = C(sI-A+K f C) _1 B 

A(s ) = I + K c (sI-A) -1 B 

A(s) = I + C(sl-A) -1 K f 


and in either case: A^(-s)A(s) = A(s)A^(-s) = I 




Table IV. 7 Summary of Plant Order Reduction Results 


Preliminary Assumptions 

G(s)=H(sI-F)“"*G + G u $(s) where Re[X[F]]<0 and G us (s) has unstable poles 

Use first steps of advanced loop 
shaping to obtain (Table II. 5) 

Regulator Loop Transfer Function 

= L.(s) 

Filter loop Transfer Function 

* L 0 (s) 

Then the indicated closed loop 
transfer function is known 

H.(s) = [I + l^s)]" 1 L - ( s ) 

H 0 (s) = L 0 (s)[I + L 0 (s)r n 

Assume realization { F , G , H} is 
(Table IV. 3) 

input normal: F+F^+GG^ = 0 

output normal : F+F T +H T H = 0 

Eigenvalue decomposition of the 
positive definite input to 
f ( ■ ) ( to be defined) 

Y = WAW T 

u = waw t 

Select r columns of W correspond- 
ing to the largest eigenvalues 

A = diagfA.. } where ^ > • • • W = [W-| W^] 

Compute the reduced order model 

G r (s) = HW ] (sI-w|fW-j ) _ 1 w|g + G(,g ( s ) 

Compute the weightings 
(Tables IV.4,5) 

W Q (s) = M i (s)G“ 1 (s) 

w^s) = G“ 1 (s)H 0 (s) 

Spectral Factorization (Table IV. 6) 

W Q (s) = A(s)W q _(s) 

W.(s) = W.Js)A(s) 

Solve for grammian (Table III.l) 

f (Y) = ^_/j-jtoI-F T )' 1 H T wJ_(ju3)W o _(ju.) 
H(ja)I-F) _1 d(o 

_co 

f (U) = 2¥/ m (jtoI " F) ’ lGW i- (ja,) 
W l | l _(ju.)G T (-ju)I-F T )‘ 1 du3 

Solve for stationary point of f(-) 

Y = f(Y) 

U = f(U) 

Use final steps of advanced loop 
shaping to obtain (Table II. 5) 

K^(s) such that K r (s)G r (s) = (s ) 

K r (s) such that G r (s)K r (s) 
= L 0 (s) 
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Plant order reduction on the other hand permits the compensator to be 
designed for the reduced order plant. This is clearly an advantage and satisfies 
the intentions of model reduction for control design. The disadvantage of this 
method is that the algorithm for implementing plant order reduction is not closed 
form (i.e. the model reduction weightings are not known a priori). In some cases 
it is not possible and in many cases it is not desirable to design a compensator for 
the high order plant and a reduced order model is required. This attractive, 
formal approach for obtaining this reduced order model is a very useful result of 
this thesis. 

The examples of the next chapter serve to illustrate the effectiveness of both 


of these methods. 
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V. EXAMPLES 

The purpose of this chapter is to illustrate by example both the model 
reduction and control design techniques discussed in the previous chapters. The 
intent is to use the examples to make specific points. Thus, typically only 
segments of the model reduction for control system design theory will appear in a 
given example. 

The chapter will illustrate three topics: model reduction, compensator order 
reduction and plant order reduction. The order of the full order model will vary 
from 2nd to 168th order. 

Model Reduction Examples 

The purpose of this section is to: verify some statements made in the 
section, “On the Question of Optimality” of Chapter LQ, illustrate graphically the 
effect of the weighting on the model reduction and to compare results for the 
balancing model reduction approach to a more classical model reduction 
approach. 

“Simple” H ^ Example 

The SISO second order transfer function for this example is given by 

2(s+2) 

(5 + l)(s + 3) 


G(s) == 
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It can be verified that the first order reduced order model of G(s ) obtained 
with the internally balanced realization technique is given by 


G r,ib( s ) — 


1.8320 
s + 1.4453 


and the smallest singular value of the balanced grammian, ct 2 is given by 
(7 2 = 0.032871. It can also be verified that the error magnitude, 
| G(jui) - G r ib {j oj) | is bounded 


| G(joj) - G r ; b (ju) | < | G(0) - <7 r> »(0) | = 2(0.032871) \/ W . 

It can also be verified that the first order reduced order model of (7(s) 
obtained with the Hankel norm technique is given by 


G rH {s) = — ~~ZT + 0.032871 
r ’" v s + 1.3028 


It is of interest to compute the error for the Hankel technique. It can be 
verified that the error is given by 

CW-I!,,W = -0.032871 . 

' r ' HK \ (i«+l)(i«+ 3)(;Vf 1.3028) 

Note that the error magnitude is independent of frequency! That is 

| G{ju) - G rH (jd)\ = 0.032871 \/w ! 

This is the smallest possible error as measured with the H ^ norm of 
G(s) - G r (s) when a D term is allowed [Ref. 14]. 

Thus this example is consistent with the comments made in the section: 
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“On the Question of Optimality” of Chapter HI. That is the reduced order 
model obtained with the Hankel technique (where D terms are required) has one 
half the error magnitude of the reduced order model obtained with the internally 
balanced realization technique (where D terms are not used). The error bound 
||C?(fi) - < 2<r[]T] 2 ] is also verified for this example. 


An Example to Illustrate the Effect of the Weightings 


A 46th order, 3 input and 3 output model representing a large space 
structure with many resonances was reduced to 28th order using both the 
internally balancing and frequency weighted balancing techniques. The transfer 
function for this example was of the form 


G(s) 


23 

E 


w 


l s + 2 + uf 


where uq = u> 2 = 0 (rigid body poles), six modes had damping ratios of f = 0.7 
and the rest had damping ratios of f = 0.005. A Bode plot for G(s ) (i.e. log 
a i \G(ju)\ vs. u, i — 1, 2, 3) is shown in Fig. V.l. 

For the purposes of illustration the model reduction weightings for the 
frequency weighted balancing case were arbitrarily chosen to be given by 

W -(«) = *.{') = W (*) = - — / 3 

(s+2) 2 


and a Bode plot of H^(s) is shown in Fig. Y.2. 





u. g. si grid 



Figure V. 2 Flexible Structure Model Reduction Weighting 
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Bode plots for a 28th order reduced order model, (^(s) using either 
balancing technique (weighted or unweighted) look very similar to that for (?(«) 
and thus are not shown. The only noticeable difference is that the peaks in the 
G(s) plot beyond cu = 30 do not appear in the plots for (?,.(«). 

Noticeable differences are however apparent in the Bode plots of the error 
G'(s) - ^(s). These plots are shown in Figs. V.3, 4 for the (?,.(«) obtained from 
the internally balanced realization and the (7 r (s) obtained from the frequency 
weighted balanced realization respectively. Notice that the singular values of the 
error are relatively flat for the unweighted reduced order model but that they are 
not flat for the weighted reduced order model. In fact, the singular values of the 
error are smaller near the peak of the weighting as expected. Note also, that this 
smaller error is obtained at the expense of larger error at higher frequencies 
where there is less weighting. 

The example clearly shows the effect of the weighting on the. model 
reduction. The example also illustrates the effectiveness of the frequency 
weighted balanced realization model reduction algorithm in tailoring the error 
with respect to frequency. 

Comparison of Balancing to Classical Model Reduction 

To facilitate the comparison of model reduction via balancing to more 
classical model reduction, an example was chosen for which a classical reduced 
order model is readily obtained. The example is a 7th order model of a spinning, 
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symmetric projectile. The 7 states are: 3 components of velocity, a, v, 'w, 3 
components of angular velocity, p, q, r and pitch angle, 9. The inputs are 2 
thrusters perpendicular to the axis of symmetry and the outputs are pitch rate, q 
and yaw rate, r. The coupling of the 3 states: forward velocity, u, roll rate, p 
and pitch angle, 9 with the other 4 states is typically very small for such a 
projectile and almost always ignored. This leaves a 4th order model. This 4th 
order model will be called the classically obtained reduced order model of the 7th 
order model. 

It is of interest to compare this reduced order model to that obtained by 
internally balancing. The 7th order model was balanced and the 4th order 
reduced order model was extracted. This reduced order model was for all 
practical intents and purposes identical to that obtained classically. Detailed 
data to substantiate this claim are given in Appendix E. 

The point here is that balancing is consistent with more classical model 
reduction when the model reduction is trivial. 

Compensator Order Reduction Examples 

The purpose of this section is to present some examples of compensator 
order reduction. Three relatively low order SISO compensators (less than 10th 
order) were designed and then reduced with the weighting technique discussed in 
Chapter IV. This technique proved to be effective in obtaining satisfactory (i.e. 
maintain closed loop stability) reduced order models of these compensators when 
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other techniques failed, such as: internally balancing (i.e. unity weightings), 

mode truncation and elimination of modes with small residues. The fact that the 
integral squared impulse response error criterion is not a good discriminator for 
model reduction for control design is also demonstrated. The compensator order 
reduction algorithm is not limited to LQG designs and was used successfully on a 
non-LQG design. 

Four Disk Example 

The system to be controlled is the subject of an experimental research 
project at Stanford (Ref. 45]. The system, indicated in Fig. V.5 consists of 



Fig. V.5 Four Disk System 


four disks (unity inertia) connected by a flexible wire (unity spring constant) with 
a motor for applying torques to the third disk and a sensor for measuring angular 
displacement of the first disk. 


The transfer function for this system with non-collocated actuators and 
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sensors is given by 


s 2 + 2c zs+ z 2 


G(s) = 



' > 

s+ a 


a 


4s" 


s 2 + 2faJi s+ c jf 




X 


s 2 + 2fw 2 s+ at 2 


wf 


s 2 -0.86s+ 6 2 


6 " 


s 2 + 2fw 3 s-f w 2 

«£ 


where the pole and zero data is given in Table V.l. 

Table V.l Four Disk System Pole and Zero Data 

u> l — 0.765 2 = 1 

ui 2 = 1.41 a = 4.84 

w 3 = 1.85 b = 5.65 

? = 0.02 

Note that the vibratory modes are assumed to have 2% damping and that the 
system is non-minimum phase. 

Note that the vibratory modes are assumed to have 2% damping and that the 
system is non-minimum phase. 

The performance and stability robustness requirements chosen for this 
example resulted in the constraints for the loop shape shown in Fig. V.6. The 
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relationship between these requirements and Fig. V.6 was discussed in 



Fig. V.6 Loop Shape Constraints For Four Disk System 

Chapter II and Fig. V.6 is just Fig. 11.13 for this example. Note, that the 
frequencies for which high gain is required are significantly less than the 
magnitude of the non-minimum phase zeros of C(s) (i.e. 0.07 « b) and thus 
the non-minimum phase zeros will have little impact on the design. 

A compensator was designed using the standard LQG loop shaping 
procedure of Chapter II (see Table H4). This compensator satisfied the low 
frequency constraint of Fig. V.6. However, because LQ loop transfer functions 
rolloff with a slope of negative one it could not satisfy the high frequency 
constraint. Hence, the compensator was modified by introducing an additional 
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lag to satisfy this constraint as well. Note that this additional lag does degrade 
loop properties in the crossover region due to the attendant phase loss of the lag. 
This degradation (45 deg. compared to 60 deg. of phase margin) was deemed 
acceptable. 

The compensator is given by 


*»(«) 


= 0.0216 


I s - f a 


f ' 

s 2 -f- 2fWjS + w 2 


1 - J 





s 2 + 2c zs+ z 2 ) 

S 2 + 1.9^5+ /3f 

^ 1 
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r y 

s 2 + 2 cw 2 s+ wf 


•> 

s 2 + 2$co$s + u) 2 
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S 2 + 1.2 S@2 S + P$ 


0.4/?3 s+ ft 2 

U+PA 

Pi 


Pi 
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where the pole and zero data is given in Table V.2. Note that the compensator 


Table V.2 Four Disk Compensator Pole and Zero Data 


/?! == 2.74 o = 0.0503 

P 2 = 2.87 
/? 3 = 2.99 

= 0.521 (additional lag pole) 

essentially inverts the plant as discussed in Chapter II. 

This 9th order compensator was reduced to 6th order with the compensator 



- 152 - 


order reduction algorithm given in Table IV. 1. The result is given by 


K 6 (s) = 0.0232 


s+ a r 





s + + cjj 




5 2 + 1.7 /? 4r s+/?| r 




4 r 


- x 


s 2 + 3$w 2r «4- w 2r 


«?r 


S 2 + 2C2S+ 2 2 ) 

s 2 +2tf 5 s+ fig 

* 1 

* I 


where the pole and zero data is given in Table V.3. 

Table V.3 Reduced Order Four Disk Compensator Pole and Zero Data 

/3 4r = 0.567 = ce T = 0.0515 = a 

h = 3.29 w 2r = 1.47 = w 2 

Note that the reduced order compensator essentially inverts only the low 
frequency behavior of the plant, i.e. the poles with | s | = u 1 , u> 2 and the zero 
with | s | = z. Results in Appendix E show that oj 1 , u 2 and z must be known to 
within 10% to guarantee closed loop stability. The zero at s — ~a r is nearly the 
same as the zero of the full order compensator at s — - a. The additional lag 
pole at s == -/? 4 of the full order compensator is approximated by a 2nd order 
pole with j s | = /? 4r = /? 4 . 
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This compensator is stable (closed loop poles given in Appendix E) with the 
full order model of the plant and satisfies the loop shape constraints of Fig. V.6. 
In addition the crossover properties of this compensator are essentially the same 
as those of the full order compensator (e.g. the compensator has a phase margin 
of 46 cleg, compared to 45 deg. for the full order compensator). 

It is of interest to note that reduced order models of K$(s) for any order 
obtained with the internally balanced realization of Kq(s) produce an unstable 
closed loop system. The same would be true for any unweighted model reduction 
method (e.g. residue technique, covariance cost, minimum integral squared 
impulse response error, ...). The eighth order reduced order model of A' 9 (s) 
obtained by using the internally balanced technique is discussed in Appendix E. 

It is also of interest to examine the integral squared impulse response error 
criterion for this example. This error criterion is defined 




1 1^(0112 

l|tf(0ll 2 


where for this example 

H,(t) = L-'[K,(s)-K 6 (s)](t) 

//(<) = L-‘|K,(s))(t) 

i.e. H e (t) is the impulse reponse of the error and H(t) is the impulse response of 
the full order compensator. The error criterion is arbitrarily normalized to make 
it dimensionless and thus have some meaning with respect to unity. 
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The value of E 2 was computed for K 6 (s), the 6th order reduced order model 
of Kq(s) obtained from the internally balanced realization of K$(s). These 
results are given in Table V.4. Note that while the reduced order model obtained 


Table V.4 Comparison of E 2 and Closed Loop Stability 


Model Reduction Technique 

E 2 

Closed Loop System Stability 

Frequency Weighted Balancing 
Internally Balancing 

0.956 

0.119 

stable 

unstable 


with the internally balancing technique has an E 2 error criterion eight times 
smaller than that of the reduced order model obtained with the frequency 
weighted balancing technique, it results in an unstable closed loop system, 
whereas the weighted reduced order model results in a stable closed loop system. 

These results clearly suggest that E 2 is not a good discriminator for model 
reduction for control system design. It is also clear that properly selected 
frequency dependent weightings are critical in obtaining a reduced order model 
for this example. 

This example is somewhat academic because there really isn’t any need to 
reduce the order of the compensator. However, the example has served to 
illustrate the procedure and make the following points: frequency dependent 
weightings are critical and E 2 is the wrong error criterion. 
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Robust Four Disk Design 

The above design for the four disk system, although satisfactory with respect 
to the constraints of Fig. V.6 and the crossover requirements (say, phase margin 
> 45 deg . ), it is not satisfactory with respect to more stringent stability 
robustness requirements. For instance, more robustness would be required if the 
disk inertias were significantly uncertain. This in turn would lead to uncertain 
natural frequencies Wj, w 2 , w 3 as well as overall gain. The above design based on 
plant inversion would be suspect in this case. 

In another study by this author the inertia of the fourth disk was 
taken to be uncertain but bounded between one and four. The pole-zero 
configurations for these two extremes are shown in Fig. V.7 



Fig. V.7 Four Disk System Pole-Zero Configurations 

Note that the inertia variation between the extremes causes the lowest frequency 
pole and zero to interchange locations in the s plane. Since the system exhibits 
such significant variation, a rigid body model was used for the following. 
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The intent of the previous study was to maximize performance 
subject to satisfactory robustness requirements (including robustness with respect 
to the uncertain 4th inertia). The loop shape requirements for this problem are 
shown in Fig. V.8. Note that the 



Fig. V.8 Four Disk Robust Design Loop Shape Requirements 


high frequency constraint of Fig. V.6 was modified. The loop is required to have 
more attenuation for frequencies greater than 0.5 compared to Fig. V.6, to satisfy 
the uncertain inertia requirement (more details are given in App. E). However, 
the loop is permitted to have less attenuation for frequencies just less than 0.5 
compared to Fig. V.6 to facilitate achieving maximum performance. Satisfactory 
crossover is also required, say phase margin should be 45 deg. or more. 
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The design technique of the previous study resulted in u 0 = 0.13 as the 
maximum performance subject to the above constraints. It is of interest to note 
that this performance is even higher than that of the previous design! An 8th 
order compensator which achieves this performance is given by 


K s [s) 


0.0639 


/ \ 


S 2 + lAlOjS-f of 

s+a 2 
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s 2 + 1.93o 2 s-f o 2 
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s 2 -b I.9303S-I- O3 
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s 2 + 0 . 52 o 3 s + O3 
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“3 
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9 

a 3 

- j 
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s 2 + 1.41 o 2 s+ af 


r \ 
s 2 + 0.52o 2 6'+ o 2 
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-n 

0 
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o 4 
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where the pole and zero data is given in Table V.5. The magnitude of the pole 

Table V.5 Four Disk Robust Design Pole and Zero Data 


Oj == 0.0850 a 3 — 0.500 

a 2 — 0.403 a 4 = 200 

at s == - o 4 was arbitrarily made very large to make the compensator strictly 
proper. 

This 8th order compensator was reduced to 5th order with the compensator 
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order reduction algorithm given in Table IV. 1. The result is given by 


A 5 (s) = 0.0638 


s 2 + 1.35o lr s+ a 2 r 


r 

fy 2 

) 

i 

' s +7i) 

f s +7 2 



7i 


72 


s 2 + 0.752a 3r 5+ a| r 


a 3r 


s“+ 0.552a 2r s + a 2r 


a 2r 


s+a 4r 


a 


ir 


where the pole and zero data is given in Table V.6. 

Table V.6 Reduced Order Four Disk Robust Design Pole and Zero Data 


a l r 

= 0.0872 = a, 

ct 4r 

= 32.1 

a 2r 

= 0.382 = a 2 

7i 

= 0.107 

a 3 r 

= 0.531 = a 3 

72 

= 0.310 


This compensator provides essentially the same closed loop properties as the 
full order compensator. That is the loop satisfies the constraints of Fig. V.8 and 
hence the closed loop system is stable for any inertia of the fourth disk between 
one and four (a root locus is shown in App. E). In addition, the compensator 
provides essentially the same performance and crossover properties as the full 
order compensator. 

Again this example did not really require model reduction, however, it shows 



- 159 - 


that the procedure of Table IV. 1 can be applied to a compensator designed with 
any method, not necessarily LQG designs. Also it demonstrates that 
performance and robustness properties are not compromised as a result of the 
reduction. 

Flexible Beam Example 

The system to be controlled is also the subject of an experimental research 
project at Stanford. The system, indicated in Fig. V.9 consists of a 
flexible beam with a motor for applying torques at one 


-;; v 

7777777r 

Fig. V.9 Flexible Beam System 

end and a sensor for measuring tip displacement at the other end. 

The transfer function for an 8th order model of this system with non- 
collocated actuators and sensors is given by 
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where the pole and zero data is given in Table V.7. Note that this system is also 
non-minimum phase. 

Table V.7 Flexible Beam System Pole and Zero Data 


Oj = 10.3 

u j = 11.3 

ft = 0.05 

o 2 = 11.7 

ui 2 = 22.0 

?2 = 0.02 

o 3 = 36.1 

tn 3 = 52.8 

<M 

o 

d 

II 

V/i 

a 4 = 37.6 




The performance and stability robustness requirements chosen for this 
example resulted in the constraints for the loop shape shown in Fig. V.10. Note 
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that the minimum phase zeros for this problem will seriously degrade the 
crossover properties for designs obtained by any technique. This is because of 
the proximity of the non-minimum phase zeros to the crossover region. 

A compensator was designed using the standard LQG loop shaping 
procedure of Chapter II (see Table 11.4). An additional lag was appended as in 
the LQG design for the four disk system. This compensator satisfied the 
constraints of Fig. V.10. The loop crossover properties were seriously degraded 
(only 24 deg. of phase margin and 5 dB of gain margin) due to the non-minimum 
phase zeros and the additional compensator lag. The compensator is given by 
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where the pole and zero data is given in Table V.8. Note that the compensator 
essentially inverts the 


Table V.8 Flexible Beam Compensator Pole and Zero Data 


0 X = 1.56 

Wj = 11.4 = U>! 

fi 

= 0.057 

= fi 

0 2 = 22.9 

uio — 22.4 = ui 2 

?2 

= 0.009 

= f 2 

0 3 = 32.2 

£j z = 52.8 = u> 3 

?3 

= 0.02 

= fe 


0 A = 53.4 
04 = 470. 

minimum phase behavior of the plant as discussed in Chapter II. 

This 9th order compensator was reduced to 6th order with the compensator 
order reduction algorithm given in Table IV. 1. The result is given by 
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where the pole and zero data is given in Table V. 9 . 


Table V.9 

Reduced Order Flexible Beam Compensator Pole and Zero Data 


°lr = 

8.66 = Oj 

Wlr = 11-6 = W, 

fir = 0-08 — fl 

hr = 

1.59 = /?! 

C^2 r 21.3 C^2 

?2r = 0.01 = ft 

7i = 

22.3 



72 = 

41.8 




Note that this compensator essentially inverts the low frequency, minimum phase 
behavior of the plant, i.e. the poles with | s | = u v u > 2 and the zero at s — - 07. 

This compensator is stable (closed loop poles given in App. E) with the full 
order model of the plant and satisfies the loop shape constraints. The crossover 
properties of this compensator were however degraded (only 17 deg of phase 
margin and 3 dB of gain margin). 
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It is of interest to compare this 6th order compensator to one that would be 
obtained by neglecting the third mode (i.e. the plant pole with | s | = w 3 ). This 
is referred to as modal truncation and is often done in practice. An LQG design 
was carried out for a 6th order plant model obtained by truncating the w 3 term 
of the partial fraction expansion of <7(5 ). The design was carried out such that 
the low frequency constraint of Fig. V.10 was satisfied to facilitate a meaningful 
comparison. The closed loop system with this compensator and the full order 
model was found to be unstable. Details of this design are given in App. E. 

Again this example does not really require model reduction, however, it 
shows that the procedure of Table IV.l can design a reduced order compensator 
that is closed loop stable with the full order model when the classical modal 
truncation technique fails. This is not meant to be an indictment of the modal 
truncation technique. It is a useful technique when used carefully. For example 
it is used successfully in a later example of this chapter. 

Plant Order Reduction Examples 

The purpose of this section is to present some examples of the plant order 
reduction algorithm discussed in Chapter IV (see Table IV.7). A couple simple 
SISO (less than 3rd order) examples as well as a MIMO (6th order) example will 
be used to demonstrate the successive approximation solution procedure. This 
technique was found to be acceptable for these examples. The section concludes 
with a significant design/demonstration example, utilizing a 168th order, MIMO 
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non-collocated model of a flexible spacecraft. The model reduction and design 
techniques of this thesis are shown to be very effective in performing this design. 


Simple Low Order Examples 

The plant order reduction algorithm was exercised with the “simple H ^ 
example” (i.e. the first example of this chapter) with transfer function given by 


G(s) 


2(3+2) 
(«+ l)(«+3) 


The closed loop transfer function was chosen to be 


H(s) = 


0.5 

5 + 0.5 


An initial weighting was determined from 

VK(s) = G-'is )H(s) 

and the successive approximation algorithm (Table IV.7) converged after five 
iterations. The resulting reduced order model is given by 


G r (s) = 


1.78 

8 + 1.36 


0 5 

The compensator, AT(s) — G'~ 1 (s) — 1 — is guaranteed to produce a stable closed 

8 

loop system with (?(s) since == 0.024 < 1. 
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Third Order Example 

The plant order reduction algorithm was also exercised with the 3rd order 
transfer function: 

= (»+ 0 . 8 )(.+ 2 ) 

( s + 1.5) (s 2 -f 1.4s + 1) 

For this example the closed loop transfer function was chosen to be 


u initial weighting was determined from 

\V(s) = G~ l (s) H{s) 

and after 3 and 5 iterations the successive approximation algorithm (Table IV. 7) 
converged for the 2nd and 1st order reduced order models respectively. These 
reduced order models are given by 


Go(«) = 


0.822(«-F 1.14) 
s 2 + 1.34s + 1.07 


G,(s) = 


1.16 

(a + 0.819) 


The compensators, K r (s) = G r 1 (s ) — r = 1, 2 are guaranteed to produce 

s 

stable closed loop systems with G{s ) since E ^ — 0.0085, 0.249 < 1 for r = 1, 2 
respectively. 


The only significance of these two examples is that the successive 


approximation solution procedure was successful. 
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C5A Example 

The successive approximation solution procedure was also demonstrated 
with a more realistic example. The plant transfer function for this example is a 
6th order representation of the C5A aircraft wing root bending and torsion 
moments measured in in-lbs due to aileron and elevator commands measured in 
radians [Ref. 46]. The MEMO system with two inputs and two outputs has the 
poles and transmission zeros given in Table V.10. A Bode plot of the singular 
values of the 2X2 6th order plant transfer function, G'(s) is shown in Fig. V.ll 
and an output normal realization of the transfer function is given in Table V.ll. 


Fourth Order Reduced Order Model 

To obtain a 4th order reduced order model, the following output closed loop 
system was used 


H 0 (s) = 


+ 1 


• An initial input weighting was determined from 

Wi( s ) ---- G-\s)HAs) . 

The successive approximation algorithm converged after 5 iterations. The 
resulting 2X2 4th order reduced order model poles and transmission zeros are 
given in Table V.12. A realization of this reduced ■ 
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Table V.10 C5A Poles and Transmission Zeros 

Poles 


s 2 ; a; 

II 
1— ‘ 

II 

0.6 

short period mode 

s 2 ; u> 

= 5.5, ? = 

0.09 

flexible mode 

s 

= -6.0 


aileron actuator 

s 

= -7.5 


elevator actuator 


Transmission Zeros 

s 2 ; u — 8.1, f = 0.69 
s 2 ; w = 11.3, f = 0.44 

where the notation s 2 ; w = u, f = f means the complex conjugate pair 

s — j 




Table V.ll C5A Transfer Function Realization 


Far + Gu G{s) = H{sI-F)~ x G 

Hx F + F T + H t H = 0 


r aileron command! __ [Bending momentl 

Lelevator commandj ^ LTorsion momentj 


-1.192e-01 
-4.412e-01 
-5.366e + 00 
7.003e-01 
-9.3l5e-01 
. 2. 980 e -02 


5.806e-01 4.758e+ 00 -1.464e+00 


-4.414e-02 
5.039e-01 
-8.856e-01 
-3.954 e -01 
-2.697e-01 


-1.014e-01 
-9.381 e -01 
1.491e-01 
-1.598e-01 
2.673e-02 


1.343e + 00 
-2. 174 e 4- 00 
-1.232 e + 00 
-4.563 e -01 
-4.245 e -01 


2.060e+00 1.640e + 00 

-4.941e-01 -5.637e-01 

4.632 e + 00 3.238e+00 

4.452 e + 00 5.533e + 00 

-6.579 e+ 00 -2.592 e + 00 
-4.385 e -01 -7.364 e + 00. 


G = 


-2.577 e+ 08 2.985e + 08' 
1. 865 e 4- 08 2.345 e + 08 

-2.491 e+ 07 -8.587 e + 07 
-1.875c + 07 -2.817e + 07 
-1.139e + 07 -1.851 e+ 07 
-3.218e + 06 -2.683e + 06. 


r-4.682e-01 

L-1.386e-01 


-2.971 e-01 
2.364e-03 

3.010e + 00 
-2.024c +00 


-1.356e + 00 
1.932 c -01 

2.784 e +001 
2.642 e+OOj 


-1.538e + 00 
-3.131e-01 
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Table V.12 

C5A 4th Order Reduced Order Model Poles and Transmission Zeros 
Poles 


s 2 ; u> == 1.25, 

f = 

0.511 

short period mode 

s 2 ; w = 5.30, 

? = 

0.0925 

flexible mode 

Transmission Zeros 



s 2 ; u -=' 10.5, 

f = 

0.395 



order model is given in Table V.13. Clearly the algorithm has just eliminated 

Table V.13 C5A 4th Order Reduced Order Model Realization 


G r (s) = C{sI-A)-'B 


A 


-1.003e-03 
-3.387 e+ 00 
-5.1466+00 
-2.574e+01 


1.140e+ 00 
-8.392e-03 
9.648e+ 00 
2.352 e+ 01 


2.218e-01 
-1.280e+ 00 
-1.320e-01 
5.203e+ 00 


1.170e-01 
-3.100e-01 
-5.580e-01 
-2.114e+ 00, 


B 


2.759e+ 05 
-8.895e+ 07 
1.1956+08 
.3.5126+08 


-7.092 e + 07' 
4.467 e+ 07 
3.905 e+ 07 
7.355e+ 06 . 


C 


M.413e-02 

L-2.261e-01 


3.680 e -02 
-3.854e-01 


1.251e-01 

-5.397e-01 


.998e-0ll 

,998e-02j 


the actuator dynamics, leaving the short period and flex mode dynamics 
essentially unchanged. 


Second Order Reduced Order Model 


To obtain a 2nd order reduced order model, the following output closed loop 


transfer function was used 
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H 0 (s) 


0.1 


s + 0.1 


An initial input weighting was determined from 

»;(») = <?-'(») • 

The successive approximation algorithm did not converge to a single point in this 
case. Rather, the algorithm resulted in a limit cycle containing three points. A 
graphical representation of the limit cycle is shown in Fig. V.12a, b where the 
balanced grammian’s singular values are plotted versus iteration. Clearly after 
the 10th iteration, the successive approximation algorithm cycles between three 
points. 

As discussed in Chapter IV, when this happens, the solution is the reduced 
order model of the limit set corresponding to the minimum error, E Xl . Recall 
that the error criterion in this case is defined 

E„ = ||[C(.)-C r (.)| C r - 1 (*)/U*)ll«, • 

Bode plots of the maximum singular value of \G{ju) - G r {ju)\ G~ l (j cj) H 0 {ju) 
versus u are shown in Figs. 13a, b, c for the three points of the limit set 
respectively. Recall, from the definition of the H x norm that is just the 
peak value of these plots. Thus Fig. V.13c corresponds to the reduced order 
model which is the solution by definition. Note that E ^ > 1 here and thus a 
control design based on this reduced order model is not guaranteed to be 


successful. 




Figure V.12a The 3 Largest Sin 









- 178 - 


This reduced order model has no transmission zeros and a 2nd order pair of 
poles with magnitude of 0.755 and a damping ratio of 0.728. Thus the 2nd order 
reduced order model retains the short period dynamics. A realization of this 
reduced order model is given in Table V.14. 


Table V.14 


C5A 2nd Order Reduced Order Model Realization 


G r (s) = C(sJ-A )~ l B 

, _ f-5.554e-01 1.142e-01 1 

* ~ L-2.344e+00 -5.438e-0lJ 

n _ r~1.186e+ 07 1.026e-f 07 1 
° ~ Lo.073e+07 1.700e+08j 

r _ r3.498e+00 1.851 e+00l 
° “ L-8.419e-01 2.896e-01 J 


This example has illustrated the use of the plant order reduction algorithm 
on a realistic example. The results obtained are consistent with intuition. That 
is for the 4th order reduced order model, throwing out the actuator dynamics 
seems reasonable. The fact that a compensator designed with the 2nd order 
reduced order model is not guaranteed to be successful is also reasonable. The 
procedure for dealing with non-convergence of the successive approximation 
solution procedure was also illustrated. 


ACOSS EXAMPLE 


The purpose of this example was to demonstrate all the steps of the plant 
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order reduction and reduced order controller design algorithm on a realistic 
example problem where model reduction was necessary. The example chosen was 
the ACOSS II large space structure [Ref. 47). The spacecraft, as illustrated in 
Fig. V.14, consists of an equipment section, solar panels and an optical structure 
to support the optical hardware for the surveillance mission of the satellite. 

The control problem is to provide acceptable line of sight errors in the face 
of the disturbance environment. Line of sight error is the focal plane x-y 
position error of the image of the optical system of this surveillance satellite. 
The disturbance environment consists of external disturbances: solar, gravity 
gradient, aerodynamic and thermal and internal disturbances: imperfectly 

balanced rotating machinery on board the equipment section (cryogenic coolers 
for the mirrors and control moment gyros). 

Sensors and Actuators 

Three actuators for solving the control problem are located near the center 
of the equipment section (Node 44). They are control moment gyros (CMG) and 
provide torque inputs about the x, y and z axes. The sensors for solving the 
control problem are located on the optical structure near the focal plane (Node 
II). They are rate-integrating gyros and provide angular position measurements 
about the x, y and z axes. Note that the sensors and actuators are non- 


collocated. 



primary 



Figure V.14 ACOSS II Large Space Structure 
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Full Order Model 

A finite element model of the spacecraft was developed at Draper Labs [Ref. 
48] and is given by 


Gm( s ) = 


84 

£ 



s 2 + 2 + uf 


where oj { , b { , c,- and jy i — 1, 2, 84 are given in Appendix G. A Bode plot 

of the singular values of G 84 [jw) is shown in Fig. V.15. That is G(s) is the 
3X3 transfer function relating the CMG inputs to the gyro outputs. The first 
six of the modes are rigid body modes and thus u>,- = 0 * =1, 2, ..., 6. Six of 

the remaining modes have damping ratios of == 0.7 (see App. G) and the rest 
have damping ratios of f,- = 0.005. The higher damping ratios reflect the 
presence of passive isolators connecting the equipment section to the optical 
structure. The lighter damping is just that due to structural damping. 


Loop Shape Constraints 

The loop shape constraints for this control problem are shown in Fig. V.16 
and would apply for either the input or the output loop shapes. 





- 183- 



Fig. V.16 Loop Shape Constraints 

The low frequency portion of the figure reflects the requirements for maintaining 
acceptable line of sight error in the presence of the disturbances. The high 
frequency constraint is due to sensor noise and high frequency modeling errors. 
The section of negative one slope is due to non-minimum phase zeros and will be 
discussed shortly. The crossover requirements are that the minimum singular 
values of the return difference and inverse-return difference for either the input or 
the output remain greater than 0.5. 

Constraints Due To Actuator Dynamics 

Actuator dynamics impose additional loop requirements which are expressed 
in terms of the input inverse-return difference. This requirement is due to the 
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fact that the design will be done neglecting the actuator dynamics. 
Requirements for closed loop stability with respect to neglected actuator 
dynamics can be generated with the stability robustness theorem of Chapter II. 

Consider the block diagram of the actuator dynamics shown in Fig. V.17. 

ACTUATOR 
OUTPUT 


ACTUATOR 

COMMAND 



Fig. V.17 Actuator Dynamics 

The actuator dynamics are nominally modeled as an identity (i.e. A ac ,(s) = 0 
nominally) and A act (s) accounts for the actual dynamics of the actuator. The 
point is that although A act (s) is not known exactly it can be bounded. The 
bound for this problem is given by 

< UctM 

and the scalar function of frequency, l act (w) is plotted in Fig. V18. To fix ideas 
this representation includes 


G act (s) = 


a 

s + a 


I 



a 


2 


G a Js) = 


s 2 + 1.4as + a 2 


I 
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log e act (co) 



Fig. V.18 Actuator Uncertainty Bound 


and 




q(s-f a) 
s 2 + 1.4as+ a 2 


I 


where a > 200 r/sec as potential actuator dynamics. 

To guarantee closed loop stability for any of the actuator dynamics included 
in the above representation it is necessary that 

c[/ + ^K{ju)G{ju)Y ] > I act M \/w 

Thus, this additional requirement is imposed on the input inverse-return 
difference. 

Motivation For Using Methodology of Chapter rV 

A compensator was designed using LQG and a model of only the rigid body 
portion of the transfer function <S' 84 (s), i-e. Deglecting the flexibility effects. This 
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design satisfied the performance imposed high gain constraint of Fig. V.16. 
However, the closed loop system was unstable when the flexible effects were 
added. Clearly the design process must take into account the flexible modes. 
But, there are 78 2nd order flexible modes to take account of ! Surely not all 78 
are important for the control design, but which ones of the 78 are critical? The 
plant order reduction methodology of Chapter IV addresses this question. 

Preliminary Model Simplification 

The three rigid body modes corresponding to translation are neither 
controllable nor observable with the actuators and sensors and thus were 
discarded (see App. G). The modes of the model with > 100 r/sec and 
1 1 6, | | • 1 1 c,- 1 1 =0 were also discarded (see App. G), Note that the modes beyond 
100 r/sec are sufficiently far beyond the crossover region and thus it should not 
be necessary to include them in the model. The closed loop stability of the 
system in the presence of these discarded modes, will however, have to be 
checked, after completing the design. The justification for eliminating the modes 
with small residues and those with magnitudes greater than 100 r/sec (roughly 
two decades greater than the desired bandwidth) is that this model reduction is 
trivial and the formal balancing would eliminate these modes anyway. 

The preliminary simplification resulted in a model with 29 2nd order modes. 
The transfer function for this model will be denoted by ( 729 ( 5 ). A Bode plot of 
the singular values of Gog(ju) is shown in Fig. V.19. Note that this plot is 
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essentially the same as that for C 84 (s) except that peaks are missing for oj > 100 
r/sec. 

Minimum Phase Approximation of G 2 g(s) 

The 58th order reduced order model, G 2 9 ( s ) h 35 1® transmission zeros in the 
right half plane. The non-minimum phase zeros all have magnitudes greater than 
6 r/sec and all but one are very close to the juj axis. It is desirable to 
approximate the system with a minimum phase model, G(s). This will satisfy 
the technical requirements for using the LQG loop shaping design procedure 
discussed in Chapter D. 

This is accomplished by factoring 

<?»(«) = G(s)P(s) 

where G(s) is minimum phase and P(s) is a MIMO all pass, i.e. 
P(s)P T (-s) = P T (-s)P(s) = I. In addition P(s) was chosen such that 
lim [P(s)-7] =0. This can always be done and a general procedure for 

B —-00 

performing this factorization is discussed in Appendix F. This of course, 
introduces additional uncertainty into the design process. The stability of the 
closed loop system with G 2 g(s) can however be guaranteed by satisfying 


e[/ + [/OHOf/w)]- 1 ] > s[p{ju) - /] 



tf w P •*■« -*-* z) “O 0.1 


G29(s> 



Figure V.19 Bode Plot of the Singular Values of G 2 g(s) 
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A Bode plot of the maximum singular value of P(ju)-I is shown in Fig. V.20 

This leads to the portion of the mid frequency loop shape constraint with -1 
slope. In other words, the loop must be rolled off for frequencies greater than the 
magnitude of the smallest non-minimum phase zero. 

Stable and Unstable Decomposition of (^(s) 

The transfer function G(s) can be expanded into stable and unstable (rigid 
body poles for this example) parts: 

<?(•) = <?.(«) + <?„(«) • 

That is 



where R~ ] is the full rank 3X3 inertia matrix for the center of mass of the 
structure about the x, y and z axes and (?„(«) represents the flexible 
characteristics of the satellite. 

Preliminary Design Steps 

A filter loop transfer function was designed using the standard LQG loop 
shaping; procedure (2nd column of Table II. 4) to satisfy the high gain constraint 
of Fig. V.16. A satisfactory high gain characteristic was obtained using only the 
rigid body dynamics, i.e. ("^(s). In other words, for this step of the design, the 

/ 



of eic-nv =j tj a* 


Haxinunun Singular Value of P<jw)-I 
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process noise driving the flexible mode states was taken to be zero. 

To be more precise let G ul (s) = C(sI-A)~ l B then the scalar p was chosen 
such that Wj(s) — C(sI-A)~ l T satisfied the high gain constraints of Fig. V.16 
where T — pB . Then the filter ricatti equation was solved for Pj > 0: 

AP f + p } a t + rr r - Pj c t cp, = o 

and the filter gain, Kj = PjC T was computed. This resulted in the filter loop 
transfer function Lj(s) — C{sI-A)~ l Kj and a Bode plot of the singular values 
of Lf(ju) is shown in Fig. V.21. Note that as discussed in Chapter II this 
transfer function also satisfies the high gain loop shape constraint, since 
W f {s) = CisI-AY'T did. 

Plant Order Reduction 

The filter loop transfer function was then used to obtain the output closed 
loop transfer function 

H g (s) = L 0 (s)[I+ L 0 (s) j- 1 
= C(sI-A+K f C)- l K f 

An initial input weighting was obtained from 

w i{*) = G~ t l (s)H c {s) 7"^° j Q0 where 

the lag at s = -100 (actually three lags since VF t -(s) is a 3 X 3 matrix transfer 
function) was needed such that the initial weighting would be proper. 



3WBP3 


C*inv<sI-A>*Kf 



Figure V . 21 Bode Plot of the Singular Values of the Filter Loop Transfer Function 
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The successive approximation solution procedure resulted in a satisfactory 

4 

32nd order reduced oder model, G r (s) of the 58th order (7(s). The iterations 
had not converged after 32 (32nd order and 32 iterations is just a coincidence) 
iterations but the values of were changing very little and hence the process 
was stopped. The error plot, i.e. a\\G[ju )- G r (j<jj)\ G~ l (joj) H(joj) } versus u> 
is shown in Fig. V.22. Note that E ^ is less than one and thus a satisfactory 
design based on (7 r (s) is expected. A Bode plot of the singular values of G r (ju>) 
is shown in Fig. V.23. 

Recall that the model reduction algorithm produces a reduced order model, 
G jr (s) of G'Js) and that G r (s) is obtained from 

G r (s) — C? er (s) + £„„(«) 

A realization of G ?r (s) is given by 

G <r ( 5 ) == C^sI-AnY 1 B l . 


A 


A 


A 0 

0 A u 


B A 


B 


C J \C C x 


Letting 
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we have 

G r {s) = C{sI-A)- l B 


Final Design Steps 

To complete the design of the compensator, i£(s) the filter loop transfer 
function Lj{s) — C{sl-A)~ l Kj is recovered with the remaining steps of Table 
n.4. That is the regulator Riccati equation is solved for P c > 0: 

A T P c + P C A + q 2 C T C - P c BB t P c = 0 
and the regulator gain is obtained from 

K c = B T P. e . 


The scalar q w’as chosen large enough to satisfy the loop shape constraints of Fig. 
V.16. 

Finally the compensator transfer function is given by 

K{s ) = K e {si -A + BK C + K f C)' 1 K f 

where 




A Bode plot of the singular values of K(ju) is shown in Fig. V.24. Note that the 
overall shape is that of the familiar lead compensator. Note also that 
considerable notch filtering is going on (i.e. plant inversion). The sensitivity of 
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this plant inversion is an issue. Hence an analysis is presented in Appendix E for 
how accurately the resonant frequencies must be known to guarantee closed loop 
stability. 

Verification of Stability 

Well the compensator is certainly closed loop stable with the reduced order 
model. G r (s) (it is an LQG design). However the stability with the actual 
system <7 84 (s) must still be ascertained. The stability of the compensator, K(s) 
with the 58th order non-minimum phase model, G 29 (s) was verified with an 
eigenvalue check (results are tabulated in App. E). Finally the stability of the 
compensator, K{s) with G 84 (s) was verified by checking that 

ll[c M («) - <?»(«)] *(«)[/ + CmMatw]' 1 iu < 1 . 

This is just the stability robustness theorem of Chapter II applied to this 
problem. A plot of (f[\G Bi (jtju) - G 2 g(;w)] + G 2 Q{j<^)K{joj)]~ l ] versus u> 

is shown in Fig. V.25. Note that its peak value is less than one and thus the 
closed loop system is stable with the full order model, Gg^s). 

It remains to verify that the neglected actuator dynamics do not cause an 
instability. This can also be verified with the stability robustness theorem of 
Chapter II. In this case the following inequality must be satisfied 


a[l + [KUv)G84Uu)]~ 1 } > UctM Vu; 



s < 
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By comparing the plot of jjj/ + [K(ju)G 84 (ju>)]~ 1 j versus w shown in Fig.' V.30 

to the plot of l aet {ui) shown in Fig. V.18, it can be verified that this condition is 
satisfied. 

In conclusion with regard to stability, the 32nd order compensator is closed 
loop stable with the 168th order full order model (7 84 (s) the presence of the 
actuator uncertainty. 

Verification of the Other Loop Requirements 

Bode plots for the input loop transfer function, /tT(s)G' 84 (s) and the output 
loop transfer function, G 84 (s)K(s) are shown in Figs. V.26 and 27 respectively. 
By comparing these plots to the constraints of Fig. V.16 it can be seen that the 
low and high frequency constraints are satisfied but the mid frequency constraint 
with negative one slope is not. Recall that this constraint was imposed for 
stability requirements due to the minimum phase approximation used for design. 
Since stability was checked and found to be satisfactory (i.e. stable closed loop 
system) violation of this constraint is of no concern. 

Bode plots of the input return difference, I + JC(s)G(s ) and the output 
return difference I + G(s)K(s) are shown in Figs. V.28 and 29 respectively. 
Recall that the loop shape constraints required that the minimum singular value 
of these return differences be greater than 0.5 for all frequencies. It can be seen 
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from the plots that while the output return difference satisfies this constraint, the 
input return difference is in violation of the constraint for a narrow band of 
frequencies near 8 rad/sec. 

Similar remarks apply to the input inverse-return difference, 
/■ + [AT(s ) )]— 1 and the output inverse-return difference 1 + [G'(s)/tT(s)]“ 1 
■with Bode plots shown in Figs. V.30 and 31 respectively. 

Note that the input loop properties are not as good as the output loop 
properties. This is a consequence of the LQG loop shaping design methodology. 
That is the desirable properties of the output loop were designed formally (i.e. 
2nd column of Table II. 4 and Table IV.7) but the input loop properties are just a 
consequence of the design. 

Step Responses 

Time histories (see App. E) were computed for step commands and were 
found to compare favorably with expected results based on the performance 
specification. 

Conclusions For ACOSS Design 

The plant order reduction and reduced order compensator design 
methodology of Chapter IV has been shown to be effective on a realistic, high 


order design example. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 
FOR FURTHER RESEARCH 

Conclusions 

The major conclusion of this thesis is that model reduction and control 
design, which have long been treated tacitly as independent design steps, can be 
successfully interrelated through a formal design process developed in the thesis. 

The formal process combines frequency weighted model reduction with 
modern linear quadratic gaussian loop shaping design. The model reduction 
approach involves a frequency weighted error between the full and reduced order 
model. An important point is that the size of this error is measured by the peak 
value of its frequency response. This criterion is important because it enables the 
size of the model reduction error and closed loop stability to be related by 
appropriate choice of frequency weightings. 

Internally balanced realizations were investigated and found to have 
desirable properties as a model reduction tool. A new theoretical result for 
bounding the peak value of the error frequency response was obtained. The 
balancing mechanics was extended to include frequency dependent weightings. 
Appropriate weightings for control system design were derived. 

Several example problems demonstrated the effectiveness of the model 
reduction technique for obtaining successful control designs. The examples 
illustrate that the frequency weightings and peak value error criterion are critical 
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for non-trivial model reduction problems. For trivial (i.e. near pole-zero 
cancellations, small residues, modes with frequencies far beyond control 
bandwidth) model reduction problems any reasonable model reduction method 
will give acceptable results. 

Recommendations For Further Research 

There are always some open questions left at the end of a research endeavor. 
Some of the open questions of this thesis and the author’s comments regarding 
them are: 

1. Are reduced order models obtained from frequency weighted balanced 
realizations asymptotically stable if both the input and output weightings 
are non-unity? The author believes that they are stable (possibly a weak 
condition on the weightings is required, say no transmission zeros on the jui 
axis). 

2. Can a simple error bound for frequency weighted balanced realizations be 
proved? The author’s opinion is that this is hopeless. 

3. Would a scalar transfer function times an identity used as an input 
weighting produce the same reduced order model if it were used as an 
output weighting? Intuitively the answer is yes and a proof should be 


possible. 
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4. Can the theoretical properties of the plant order reduction fixed point 
problem be established (existence and uniqueness of solutions)? The author’s 
opinion is that solutions do exist and probably more than one in many cases. 

5. What are the similarities and differences of plant and compensator order 
reduction in terms of the resulting control system design? It would be 
pleasing if the results were the same but the author’s opinion is that they 
won’t be the same. 

6. Can the Hankel results be generalized to include the frequency dependent 
weightings? Certainly a weighted Hankel model reduction can be done. The 
interesting question is: can an error bound be proved? 

7. Can the plant order reduction problem be solved with something other than 
the linear quadratic gaussian synthesis procedure? Combining the model 
reduction ideas with the recent results in H ^ optimization is a fruitful area 
for several Ph.D. dissertations. 

8. Can better numerical algorithms for solving the model reduction problem be 
obtained? Any progress in this area would be of great value to the 


practicing control engineer. 
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9. Can model reduction weightings be chosen for purposes other than 
guaranteeing closed loop stability? Certainly weightings can be chosen for 
other purposes e.g. to satisfy performance objectives. 
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APPENDEX A 

SINGULAR VALUES 

The singular values of a complex, n X m matrix, A are denoted by (7,- [A] 
(or just cr,- when context permits) and are defined to be the k largest, 
nonnegative square roots of the eigenvalues of A n A (or equivalently AA n ) where 
k == min (n,m). That is 

a ; \A\M^[A B A\ ( = 1 , 2 ,...,* 

and the ordering is such that cr { > a i+ j. 


The maximum and minimum singular values ofA, denoted by cr[A ] and q\A\ 
respectively are equivalently given by 


a[A] 

z{A] 


max 

xj^Q 

min 

xj^O 


[ Ax | 


= max 


Ax 


I 1*1 
I Aar I 


11 * 11=1 11*1 
I Aar I i 


mm . . 

11*11 = 1 11*1 


If A" 1 exists, the minimum singular value is also equivalently given by 




1 

afA- 1 ) 


By the definition of singular values it is also clear that 

a(aA] = | a | a[A] 


for any complex scalar, a. 

The complex, n X m matrix, A can be decomposed in terms of its singular 
values, < 7 ,- as follows 
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k 

A == UZV H = £j <TiU { vP 

1=1 


where 


U = [«j u 2 • • • «„] 
F = [ Vl v 2 • • • v m \ 



[Ej 0] n < m 


and E t = diag {cr v er 2 , cr*}. 


The columns of U and V are right eigenvectors of AA H and A n A, 
respectively and are known as the left and right singular vectors of the matrix A . 
In addition the matrices U and V are unitary, that is 

U H U = UU H = I 
V H V = VV H = I . 


Several standard inequalities and theorems involving singular values are 


stated in Table A.l. 
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Table A.l Singular Values: Some Inequalities and Theorems 

Inequalities: 

. a[A+ B\ < a[A\ + <?[B\ 
a[AB\ < cr[A\ a[B\ 

<z[A] -1 < a.[I+ A\ < s[A\ + 1 

z{a\ < IM^II < 


Theorems: 

<l[A], s[B\ 0 => s[AB\ > s[A\s[B\ 
3[E] < s[A] => s[A+ E\ > 0 
ff[A\ < 1 <l[I+A] > l-Jf[A] 
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APPENDEX B 

STABILITY ROBUSTNESS CONDITIONS 

This appendix will give proofs that either of 

1) <L [/ + G(ju)\ > ?[A(/w)] \fu 

2) <l\1 + G-'Uu)} > U[G- l (ju)A(ju)} \/u 

3) sl{I 4- G~ x {juj)\ > 5 r [A( < /w)G _1 (ya;)] \fu 


implies that there exists an L[s) and an R(s) such that 


||L(s)A(s)/?- 1 (s)|| 0O < 1 


||/?( 5 )[/+ G(s)\~' L-'(*)\\oo < 1 . 

Proof for 1): 

z[/ 4- G{ju)\ > 5[A(/w)| Vw 

** W+ w < 1 

<=$> ff([/ 4- C(;w)] _1 ] cr[A(;w)] <1 \fu 

=$> ||[/+ G'(a)] -1 A(a)|| 0O < 1 

||[/+ G(s)}~' A(«)IL < 1 

L(s) — A _1 (s) and /?(a) = I 


let 


Q.E.D. 
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Proof for 2): 

<l[I + G _, (yw)| > a[G~\juj)A(ju)] Vw 

<=> — — ~rr . 7 cr[G- l {ju)A(ju)\ < i \/w 

s[I + G *(;w)] 

a[[I + a[G~ l (ju)A(ju)\ <1 \fu 

=> cf[[ I + G _1 (;w)] _1 G~ l (ju!)A(ju)\ <1 \fu> 
a[[I + G(;o;)]- , A(ya;)j <1 \/u 

||[/ + G'(*)]- 1 A(3)|| 0O < 1 

let L(s) = A- 1 (s) and R(s) = I Q.E.D. 

Proof for 3): 

£Z[/ + C? _l (;w)] > ?[A(/w)G -1 (/u/)] \/w 

^(ywJG-^yw)] — ; ~T—~7 < 1 Vw 

a[l + G *(;wj 

«S==> a[A {j(jj)G~ l (ju)\a\[I + G _I (yw)] _1 ] < 1 Vw 
=> <j{A(ju)G-'(ju)[I + G-^yw)]- 1 ] < 1 Vw 

<=$► cf[A(ju)[I + G(yw)]- 1 ] <1 \/u 

l|A(s)[/ + G( 5 )]- , || 00 < 1 


let L(s) = I and R(s) = A(s) Q.E.D. 
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APPENDEX C 

DETAILS OF BALANCING THE CONTROLLABILITY AND 
OBSERVABILITY GRAMM LA NS WITH EIGENVECTORS 

OF THEIR PRODUCT 

Lemma Cl: Product of Positive Semi-Definite Matrices is Similar to a 
Positive Semi-Definite Matrix 

Let A and B be positive semi-definite matrices, i.e. A ,B >0. Since A is 
positive semi-definite it has the singular value decomposition: 

A = VKV T 

where 


and 


Let 


At 0 
0 0 


A, > 0 


vv 


T __ 


V T BV = 


*^11 ^12 
s n s 22 


where 5 n has the same dimensions as Aj. Note that: 


B > 0 => S > 0 => 5 n > 0 


Hence 5 n has the singular value decomposition: 
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S u = UZU T 


where 


and 


Let 


E 

U 


E t 0 

0 OJ 

Wt U 2 ] 


uu 


T _ 


/ . 


^11 $12 


U 2 

sj 2 s 22 ' 


X 


C(X) = [uj X T ] 


Note that: 

S >0 => C(X) >0 MX . 
From the singular value decomposition of 5 jj it can be seen that 

U?S U = 0 


and thus 

C(X) = u[s 12 x + x t sJ 2 u 2 + x r s 22 x . 

From this expression for C(X) and the fact that C(X) > 0 for any X it can 
be seen that U 2 S X2 — 0 or else it would be possible to find an X such that C(X) 
was not positive semi-definite. To be precise, say U 2 S l2 ^ 0, then there exists a 
z such that 8 A z T U 2 S l2 Sf 2 U 2 z > 0. Next let 7 A z T U 2 S x2 S 22 Sj 2 U 2 z, then 
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S > 0, implies S 2 2 > 0, implies 7 > 0. Then let 


7 

1 


7 7^ 0 
7 = 0 


Finally let X = -eS^(/ 2 , then 


2 r C(X)z 


— 2e<5 -f € 2 7 


v 


_£ 

7 

-2<5 


7^0 
7 = 0 


In either case 2 r C(X)z < 0 and hence if 6bJS 12 7^ 0 it is always possible to 
choose X such that C(X ) is not positive semi-definite. Since C(X) > 0 for \fX 
it follows that UoS l2 — 0. 


Let 


T A V 



-Ufii'UfSa 

I 


then it can be verified that 


T' X ABT 


A'/'SiiA, 1 / 2 o' 
0 0 


This completes the proof since 5 U > 0 and > 0 imply that T X ABT > 0. 
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Lemma C2: Nonzero Off-Diagonal Elements of Two Symmetric 

Matrices Whose Product is Diagonal Correspond to Repeated Elements 
of Their Diagonal Product 

Let A with elements, a- and B with elements, 6, ; - be nXn symmetric 
matrices, i.e. A=A T and B=B T where AB == A = diag {X,}. Then it will be 
proved that a,y(X,— X ; ) = 0 and 6 |; (X,-X ; ) = 0 for \f i,j. 

Firstly note that it is sufficient to prove the result for A since 
BA = B t A t = {AB) t = A . 

Next partition A, B and A compatibly: 


A 


A 11 A 12 

A 12 A 22 


B = 


B li 

B \2 5 22 


A = 


A, 0 
0 A 2 


where A n , B n , A t are rXr with 1 < r < n but r is otherwise arbitrary. The 
equation AB = A is then equivalent to the following four equations: 


*4 it + 

A 12 &12 

= Ai 

( 1 , 1 ) 

An B\2 + 

A 12 B 02 

= 0 

(1,2) 

A 12 B\\ + 

A 22 B i2 

= 0 

(2,1) 

A 12 B 12 + 

A 22 B 22 

— 

(2,2) 


Multiplying the transpose of equation (1,2) by A 22 on the left and equation 
(2,1) by A a on the right and then subtracting the two equations results in 


A 22& 22^12 — Ai 2 B n A n 


• 7 * 

Using this result to substitute for in the equation obtained by 
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*p 

multiplying the transpose of equation (1,1) by A x2 on the left results in 

A 22 B 2 * A i2 12 B 12 A 12 == -Aj^Aj 

{^22 B 22 + A 12 B 12 J A 12 — A 12 Ai 

^■2 A 12 — .-A 5 A 1 

where equation (2,2) was used to substitute for the term in brackets. Examining 
the final equation for all r between 1 and n implies that 

«,;( X < " x y) = 0 V *'»; 

Q.E.D. 

Lemma C3: Transformation Always Exists Which Diagonalizes Both 
Grammians, U and Y 

Lemma Cl shows that a transformation, T 1 always exists such that 
T x 1 UYT l = A = diag {X, }. In fact is any eigenvector matrix of UY. Let 

A A rfUTl 7 

and 

B M TjYT x 

then 

AB = T x 1 UT x t T{YT 1 . 

Noting that A = A T and B = B T , Lemma C2 shows that the elements of A 
and B, a {} - and b {] - respectively satisfy 
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a, 7 (X,. - X y ) = 0 

and 

hi (\ - >>,) = 0 . 

Clearly if X,- X ; - for » =p^ ; then Tj is a transformation with satisfies 

Lemma C3. Otherwise without loss of generality,, order the elements of A such 
that 

Xj = X 2 — ■ • • = X Ai > X* i+ j = X * i+2 = • - • = ^ki+k, > 

^ ^ k x + k 2 + ...+ k,.,+ l " ^k x + k 2 + ...+ k t 

That is A has k { elements equal to \ kl+ ki+ + k for i=l, 2, q and 
? 

n = t, A:,-. 

1=1 

Then Lemma C2 shows that 

A = block diag {A,} 

B = block diag {5,-} 

A = block diag jx* i+ >+ k> I k% 
where A { and B { are k { X k { . Then let Vy be any eigenvector matrix for A } -, i.e. 

v i lA i v i = dia ? !«*,+*,+ ...+ *,.,+ .• 

The fact that Aj can be diagonalized is a consequence of the symmetry of Aj 
which is a consequence of the symmetry of A . Another consequence of the 
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symmetry of Aj is that Vj can be chosen such that V- = Vj T . 

Next let 

T 2 = block diag { 7,-} 

then let 

A = T 2 l A T 2 t 

and B = r 2 r B T 2 . 

Then A — diag {a, } and 

AB = T 2 l A T 2 t T'£ B T 2 
= T 2 l AB T 2 
= To 1 A T 2 
= A 

* a * 

but since A and A are diagonal B must also be diagonal, i.e. B = diag {/?,}. In 
conclusion is a transformation which satisfies Lemma C3. 

Q.E.D. 

Choice of Scaling Such That Controllable and Observable Portions 
of Balanced Diagonal Grammains Are Equal 

It has been shown in Lemma C3 that there exists a transformation, T such 


that 
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T~ l UT~ T = diag {a,-} 
T t YT = diag {/?,•} 
T-'UYT = A = diag {X, } 


and X,- = a * =1, 2, N. Without loss of generality assume the X,- 

have been ordered such that 


^1 ^ ^2 ^ ^ > ^n+l 


^n+! 



Let 


and 




V 


1 


Pi ¥= 0 

a,- =0 or Pi — 0 


D — diag {</,} 


and 

T b = TD 


then 



T b - l UT b - T 


T b T YT b 


E 


D~ x T- x UT- t D~ l 
D~ l diag {a, } D~ l 



block diag {E, E y } 


DT t YTD 
D diag {P{}D 
diag {/?,- rf, 2 } 
block diag {E, E y } 


diag |\/X7, y^ 2 , •••, \/X„ | 

dia g | a n+l. a n+2> -» a N j 
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APPENDEX D 

CONJECTURE REGARDING E x BOUND FOR FREQUENCY 
WEIGHTED BALANCED REALIZATIONS 

The objective of this appendix is to carry out the manipulations leading to 
the conjecture about the E ^ bound for frequency weighted balanced realizations. 
The notation and results of Chapter III for the proof of the unweighted error 
bound for internally balanced realizations will be used. 

F rom Chapter III we have that 

£(«)-<?,(*) = C( a JA-V )£?(*) 

Let E{s) A VF ( ,(a)[C(s) - G r (s )j H^s) where ^(a) and are the specified 

output and input model reduction weightings. Letting 

BJs) A B(.)WM 

and 

c.w J »v.(.)a(*) 

it can be seen that 

£(») = . 


From the definition of the maximum singular value we have 
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= X m(x [&~\ju)BJju)Bj(-ju)&~ T Hu)Cj(-ju)CJju)} 

The fact that the realization is balanced with respect to the model reduction 
weightings: ITJs) and V^-(s) is then used to derive expressions for 

B w {ju) B T (-ju) and Cj(-ju) C w (ju). This requires an enormous amount of 
algebra which will only be sketched here. 

The basic idea is the same as that for the proof of the bound for the 
unweighted case given in Chapter ID. That is the two grammian equations: 

A { U + UA? + B { B? -■= 0 
AjY + YA 0 + CjC 0 === 0 

are partitioned into six equations each, by substituting 

_ . -4-n A l2 B x Hi _ 

Af — - A 21 A 22 B 2 ^^ B{ == 

0 ()” Ff 

_ A 12 0 

A 0 = -4 21 A 22 0 

G.C X G t C 2 F 0 


_ Sj o y 13 
y = o e 2 r 23 

XL y'L y M . 

The result for B w (s) will be discussed and the result for C w (3) will follow by 



C, = [z),C, d,c 2 j/„] 


’ b , d - 

b 2 d { 

G { 



duality. Using the definition of B v {s) and B(a) to expand 5 w (s) Bj(-s) results 
in terms involving B^D^D^Bj , B , BiH { U zz , B 2 D i D- r Bl , B 2 Hi f/ 33 and 
GG t . Substitutions are made for these six terms using the six equations for the 
weighted controllability grammian. 

This results in an expression involving S 2 and U 22 (the E 1; U zv U zz terms 
exactly cancel). This final expression and its dual for Cj(-a)C ,(«) are given by 

BJs)Bjt-s) = A(s)£ 2 + S 2 A r (-s)+ NM 
Cj(-s)CJs) = A T (-»)£ 2 + E 2 A(») + NJa) 

where 

NM M A(a) Uj 2 {-sI-Ffr'Hi T B T {M+ BUWM-Fir'U^l-') 

NM & A T {-<,)Y a UI-F.)-'G.C(‘)+ C r MOj{-st-Fj)- l Yl A(») . 

With these expressions for B v (a)Bj(-s) and Cj(-a)C w (s ) we have 

*[E {}<*>) 1 = [[ S 2+ A- 1 (;w)E 2 A r (-;'w) + 

js 2 + A- r (-yw)E 2 A(;'w) 

+ A - T (-j to)N 0 (j u>) j j . 

If N { (ju>) = N 0 {ju) — 0 then with the same steps as in Chapter HI it can be 


shown that 
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however, N { (ju) and N 0 (joj) are not zero in general. A bound for the general 
case could not be found but a conjecture is that 

£00 = sup &[E{jw)\ < 2(1+ a) <r[E 2 ] 

U) 

where the conjecture is that 

a < 1 when E ^ < 1 
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APPENDIX E 
SUBSTANTIATING DATA 


Spinning Projectile Data 

A seventh order realization of the rigid body dynamics of a spinning 
projectile is given by 

x — Ax + Bu 
y — C'x 

where 

x T = [« v w p q r 9] 

* T = I F, F,] 


y T = I? r 


The matrices. A. B, C are tabulated in Table E.l. 


The units of the velocity components, u, v, w are measured in terms of the 


magnitude of the velocity, V ( = 1053 ft /sec). That is — which is dimensionless 

and — and which represent sideslip and angle of attack in degrees 

V 7T V ir 

respectively. The units of the components of the angular velocity p, q and r are 
measured in degrees/sec, the units of the pitch angle, 9 are measured in degrees 
and the units of F y and F. are measured in pounds (lbs). 



Table E.l A, B, C Matrices For Spinning Projectile 



-4.757p-OP 

8. HI 7 p-04 

P.50 4 P-05 

0. OOOp+OO 

-4.01 Pp-04 

-8.31 Op-04 

-4.6162-04 


-5.77/2-01 

-1.554p-0| • 

-! . 3/6a-0? 

->.?85e-07 

0. OOOe+OO 

-0.853e-OI 

-5.350P-04 


- 1 .?n4e+no 

1.4102-0? - 

-1 .551 2-01 

-o. 30/2-07 

0.C)f<6p-f'1 

?.74«e-o? 

1 .4162-0? 

A = 

-3.855e+OP 

.3 . p08p-0I - 

-1.6 40p-01 

-5. 44/2-03 

0. OOOp+OO 

0. OOOp+OO 

O.OOOp+OO 


8.621e+0? 

-4 .036e+'V> 

6. 360p+n? 

-1 .4 30e-m. 

-2 . P00e-0! 

- ! . p /8e+0P 

- 4 . 06Pe— 04 


1 . /?Qp+ CM 

-6.433P+O0 - 

-7.5/6p+ oo 

-P.O(V)p-oi 

1 .2 78 p+ OP 

-?.OO0e_M1 

-4 .06?e-04 


0. OOOo+OO 

0. OOOe+OO 

O.OOOp+OO 

0. OOOp+OJ 

1 . 0002+ 00 

0. OOOp+OO 

O.OOOp+OO 


0. 000 p+00 

O.OOOp+OO 







1 .7012-0? 

0. OOOe+OO 







O.OOOp+OO 

1.701 e-0? 






B = 

O. OOOe+OO 

O.OOOp+OO 







O.OOOp+OO 

4 . 584 p+OI 







-4 . 584e+0| 

0. OOOo+OO 







o. nnop+on 

o. rvm p+no 






0 = 

o.fionp+on 

n. ooop+oo 

o.ooop+no 

0. OOOo+OO 

I . 0002 + no 

0. OOOe+OO 

0. OOOp+OO 


| o.nnop+m 

o. ooop+oo 

O.OOOp+OO 

O. 0002+00 

O. MOOp+OO 

1 # OfVtp+OO 

0.0002+00 
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Fourth order reduced order models were obtained with the internally 
balancing technique and the “classical” technique of neglecting the u, p and 9 
states. These two reduced order models are compared below. To facilitate the 
comparison of the two reduced order models and the full order model let the 
transfer function relating u(s) and y(s) be parameterized as shown in Table E.2. 
Then the differences between the two reduced order models can be seen by 
examining the data in Table E.3. 

Note that all the parameters differ by less than a few percentage for the two 
reduced order models. The only exception is the precession damping for which 
the internally balanced reduced order model is in closer agreement with the full 
order model than the classical reduced oder model. The authors opinion is that 
for all practical intents and purposes the two reduced order models are identical. 

Four Disk Example 

Sensitivity With Respect To Knowledge of uq, w 2 and z 

The fact that the compensator, K 6 (s) involves pole-zero cancellation near 
the jjj axis is a concern. The poles of G(s) with | s | = uq, u> 2 and the zero of 
G(s) with |s j =£ are cancelled. Hence, an analysis was performed to assess 
the stability of the closed loop system in the face of simultaneous ± 10% 
uncertainty in these three critical parameters, uq, w 2 and z (sensitivity to other 
parameters would be less severe). The analysis was performed assuming that for 
any uncertainty the DC gain of G(s) was known exactly. 



Table E.2 Parameterization of Spinning Projectile Transfer Function 


G ( s ) 


1 

< s2+ v P s H;) (s2+2 w + “n> 


9'|'j(s)(s +2c-j i w^s+u)-| i ) 


-9 21 (s)(s+a 21 )(s-b 21 )(s+c 21 ) 


For Full Order Model: 


let d(s) 

= (s+5.2e-3)[s 2 +5.8e-2s+l.le-3] 

then g-j i (s ) 

= K 11 s(s+4.5e-3)(s-2.3e-2)/d(s) 

9-| 2 ( ^ ) 

= K 12 s(s+5.5e-3)(s+4.5e-2)/d(s) 

9 21 ( s ) 

= K 21 (s+5.2e-3)(s 2 +6.1e-2s+l .le-3)/d(s) 

9 22 ( s ) 

= K 22 (s+5.6e-3)(s-6.1e-2)(s+0.1)/d(s) 


9l 2 ( s )( 2)( s “b-j2)( 2 ^ 


2 2 
H22 kS J vS + 2C22 U) 22 S+w 22 / 


For Reduced Order Model : 
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Table E.3 Comparison of Full Order Model, Internally 

Balanced and "Classical" Reduced Order Models 


Description 

Parameter 

Full Order 
Model 

Internally 

9a.laneed 

Classical 

X Difference In 
Reduced firmer Models 

precession 

CP 

0.0045631 

0.0045631 

0.005601 7 

-PP.76 

ooles 

wp 

5.1790 

5.1780 

5 . 1 860 

-0.15 

nutation 

{n 

0.002 78? 4 

n. 0027824 

0.nn27816 

0.03 

poles 

<*»n 

1 PP. Vi 

1 PP. V 

! PP.65 

0.01 

Fy to q 

K 1 1 

5«50.p 

585Q.P 

5850.? 

0 

numerator 

{11 

0 . 708 1 3 

n. 70741 

0.711 71 

-0.61 


<e| 1 

0.3?365 

0.P66I 1 

0.P67R2 

-0.64 


K1 P 

45.94? 

■ 45.94? 

45.94? 

0 

Fz to q 

a12 

0.30231 

0.37564 

0.30344 

-3.64 

numerator 

bl 2 

24.052 

24.055 

24.082 

-0.11 


ClP 

P5.1I0 

P5.3P2 

25.346 

-0.00 


K?l 

45.94? 

45.84? 

45.84? 

0 

Fy to r 

a2l 

O.3052Q 

0.30935 

0.30483 

0.98 

numerator 

bPI 

P5 . 000 

?5 . 000 

?5 . 000 

0 


c2l 

25.06Q 

25.360 

25.383 

-0.06 

Fz to r 

K?2 

5950.5 

5850.6 

5850.6 

0 

numerator 

(P2 

n . 60610 

0.71724 

n.7on56 

2. ->3 

MI MO 


0.34905 

0.3P859 

0.3477? 

-5.93 

Transmls- 

{ 

O.OQ89* 

0. 00056 

0.00017 

0,04 

slon zeros 

<1 1 

0.3030? 

0.39750 

0.3038? 

-1.35 


The x it ** or-’nc® was commuted by subtraction the classical osrw» ter 
fror> the Internally halsnceri oarameter, dtvldlnn by the Internally 
balancer) parameter and then multiolytna by inn. 
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Consider the multiplicative perturbation implied by the uncertainty in , 'jJo 


and z 


A( S ) = 


G(SJ CJj, CUo. -) ~ G(sj Gj, UJo, *) 
G(s; d>i, d> 2 , z) 

s 2 + 2c zs 4- z 


fi 

UJr) 

2 f 

[Z Cjy 

(I/O 

l 


j 2 + 2 $zs + 


O - « „ O 

•$“ 4“ 4- 


’ , 
0 _ . A 0 
S~ 4“ 2fCJoS *4* GJn 

•$“ + 

< J 


S w 4“ 2^2^ 4“ GJo 

k J 


where G'(s ; wq, oj 2 . z ) represents a possible true system and G(s ; t . Go z ) 
represents the nominal system. Note that the other poles and zeros cancel out of 
the expression for A(s). Also note that A(0) = 0. 

Let /(vj) be a bound for the magnitude of the perturbation, i.e. 


/(vj) = max 

0.9 < — < 1.1 

a 

( *« = <i/j, / 


A (jco) 


The bound, l(^j) is plotted versus in Fig. E.l. 

The stability robustness theorem of Chapter II (see also App. B) guarantees 
stability of the closed loop system in the presence of the uncertainty if 


T i-! 

1 -f j^G'fya;; Gj, Go, z) K§(ju)^ 


> /(c 


Va; 


This condition can be seen to be satisfied by examining Fig. E.l where the left 
and right hand sides of the inequality are plotted versus u. Thus the closed loop 




FREQUENCY (RADIANS/SECONO) 

Fig. E. 1 Stability Robustness Test For Four Disk System With K g (s) 
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system is stable in the face of simultaneous ± 10% uncertainty in co h w 2 and 2 . 

Closed Loop Poles 

The closed loop poles for the control system with G(s) and 7F 6 (s) are given 
in Table E.4. Note that all have negative real parts which is guaranteed by the 
fact that E <*, < 0.2 (i.e. less than unity) for the weighted model reduction of 
K(s) to K 8 (s). 

Results With Internally Balanced Model Reduction 

An eighth order reduced order model of K(s) was found using unity 
weighting or the internally balanced method. This compensator produces an 
unstable system with closed loop poles shown in Table E.5.. The fact that unity 
weighting produced an unstable closed loop system is easily seen by comparing 
the Bode gain plots of the three compensators K(s) or Kg[s), A'g(s) and K 8 (s) 
shown i^Fig. E.2. While K s (s) (unity weighting) captures the high gain portion 
of K(s) nicely it sacrifices low frequency model fidelity, whereas, K s (s) (non- 
unity weighting) with only six poles matches the low frequency behavior and 
produces a stable closed loop system. 

Robust Four Disk Design 

For zero damping the transfer function of the four disk system as a function 
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Table E.4 Four Disk System 

Closed Loop Poles for G(s) and K g (s) 

(Non-unity model reduction weighting) 


real imaqinary 

1 -7.30749-0? 6 . 44 1 5e-02 

2 -7.30749-02 —6 .441 5e-0? 

3 -4.0331 9-0 1 1 .36819-01 

4 -4.03319-01 -1.8681e-01 

5 -1.46R9e-02 7.6457e-01 

6 — 1 ,4689e— 02 —7.64579—01 

7 -2 . OOOOe-OP 9.99R0e-01 

8 -2 . oo00e-02 -9.09809-01 

9 -?.9580e-02 1.40829+00 

10 -2.9580e-02 -l.4082e+no 

11 -4.24719-02 1.84709+on 

12 -4.2471 e-02 -l.3470e+00 

13 -6.51269-02 3.28959+on 

14 -6.51269-02 -3'.0395e+00 


f r equen cy 

damping 

9.741 2e-02 

7.501 6e-01 

9.74129-02 

7.501 6e-oi 

4. 44479-01 

9.0739e-0l 

4 . 44479-01 

9 . O739e-O1 

7.6471 e-01 

1 .9209e-02 

7.6471 e-01 

1 .92099-02 

1 .0000 e+00 

2. 00009-0-5 

1 . ooooe+on 

2 . O 000 e-o? 

1 .40859+00 

2. 1 001 9-02 

1 .40859+00 

2. I 001 e-02 

1 .3475 9+00 

2 .2988 9—02 

1 .84759+00 

2.298Re-0? 

3.2901 9+00 

1 .97949-02 

3.2901 e+00 

I . 9794 e-02 


Table E.5 Four Disk System 

Closed Loop Poles for G(s) and Kg(s) 

(Unity model reduction weighting) 


real Imaginary frequency damping 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 
1 3 

14 

15 
1 6 


1 .4037e-()2 
1 . 4037e-02 
-1 .27469-02 
-1 .2746 e-02 
-1 .9501 e-02 
-I .9501 e-O? 
-2. 7707 e-02 
-2. 7707 e-02 
-3 . 687?e-02 
—3 ,6872e— 02 
-1 .4981 e+00 
-1 .4981 e+00 
-5 . 9795e-01 
-5 . 9795e— 01 
-3.0948 9+00 
-3. 09 48 e+00 


1 .67859-01 
-.1 . 6785e-0 1 
7 . 6926e-01 
-7. 69269-01 

1 .001 2 e+00 
-1 .0012 9+00 

1 . 4088e+00 
-1 . 4088 9+00 
1 .34969+00 

-1 .8496 9+00 

2. 1 8739+00 
-2. 1873e+oo 

2,93 55 e+ 00 
-2.9355e+00 
1 . 9086 e+00 
-1 . 9086e+00 


1 .68439-01 
1 .68439-01 
7. 6 93 7 e-01 
7.69379-01 
1 . 001 4 e+00 
1 .0014 e+00 
1 .4091 e+.OO 
1 . 4091 p+ oo 
I ,8500e+00 

1 .8500 e+00 
9.651 2p+00 
2 .6512e+00 
2 . Q953 e+00 
2. 9958 e+ 00 
3.6360 9+00 
3. 6360 e+00 


-8.33339-02 
-8. 3338 e-02 
1 .6567 e-02 

1 .65679-02 
1 .94749-02 

1 . Q 4 74 9— 02 

1 .96639-02 
1 .96639-02 
1 . 9Q32e-02 
1 .9932e-02 
5.^5089-01 
5 ,6508e-oi 
1 . 0960e-oi 
I . 9960e-01 
8.51159-01 

8.5.1 1 5 e-01 
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of the inertia of the fourth disk is given by 


G 


1 

, s: T 

l r 


(s 2 + -T 2 ) 


s 2 js 6 + (5-t- -t 2 )^ 4 “t* (6+4 t)s 2 + (1+37 2 )j 


where is the inertia of the fourth disk. A partial fraction expansion of 

r 

G'o(*: “T) is S iven b y 

r 


Gofs ; 




i ?. 


3 

E ^ , 

( -=0 8 + Wf. 


where the residues, /?, i = 0, 1, 2, 3 and w,- » =1, 2, 3 are functions of the 

inertia of the fourth disk and = 0. The constant damping ratio f = 0.02 was 
introduced to produce the transfer function 


O, «; A 

r 


3 

= E 


% 


i=o 5 + 2fu/ t -s 4* oJf 


Consider the multiplicative perturbation given by 


A(») 


r 


<?(•) 


G(s] 


where the nominal model, C(s) was chosen to be given by 


G(s) 


0.2 


Let l(u>) be a bound for the magnitude of the perturbation, i.e. 
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l(u) == max |A(;w)| 

1 <" 7<4 

r 

This bound is plotted versus oj in Fig. E.3 and resulted in the high frequency 
constraint of Fig. V.8. 

The stability robustness theorem of Chapter II (see also App. B) guarantees 
stability of the closed loop system in the presence of the uncertain but bounded 
inertia of the fourth disk if 

1 1 + |G(ju)/f 5 (jw)|-'| > Ifu) Mu 

where K$(s) is the compensator of Chapter V. This condition can be seen to be 
satisfied by examining Fig. E.3 where the left and right hand sides of the 
inequality are plotted versus oj. 

It may also be of interest to examine the closed loop pole locations as a 
function of the inertia of the fourth disk. A root locus plot is shown in Fig. E.4 

for 1 < < 4. Note that all the poles are in the left half plane as guaranteed 

r 

by the stability robustness theorem. Also note that a closed loop pole at 
s = - a ir = -32.1 is not shown. 

Flexible Beam Example 
Closed Loop Poles 

The closed loop poles for the control system with G(s) and A 6 (s) are given 
in Table E.6. Note that all have negative real parts. 





Figure E.4 Root Locus For Robust Four Disk Design 
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Table E.6 Flexible Beam Closed Loop Poles for G(s) and K g (s) 
(Non-unity model reduction weighting) 


rea 1 

imaginary 

frequency 

damping 

1 -2. 85 33 e+00 

0. 0000e+.00 

2.8533e+00 

1 . 0000e+00 

2 -9.9263e-01 

4.6585e+00 

4 .7631 e+00 

2 .0840e-01 

3 -9. 9 26 3 e-01 

-4. 6585e+00 

4.7631 e+00 

2 . 0840e— 01 

4 -6. 35 22 e+00 

0.0000e+00 

6.3522 e+00 

1 . 0000e+00 

5-1.1 824e+00 

1.1811 e+01 

1 . 1870e+0I 

9.961 Oe-OP 

6 -1.1824e+00 

-1.1811 e+01 

1 . 1 870e+01 

9.961 0e-02 

7 — 9. 8649 e— 02 

2. 1721 e+01 

2.1721 e+01 

4.541 7e-03 

8 -9 . 8649 e— 02 

-2.1 721 e+01 

2.1721 e+01 

4.541 7e-03 

9 -2.9323e+00 

2.1 774e+01 

2.1971 e+01 

1 .3347e-Ol 

.10 -2.9323e+00 

-2.1 774e+01 

P.1971 e+01 

1 . 3347e-0| 

11 -1 . 1634e+0O 

4.21 74e+oi 

4.21 90e+01 

2 . 7575e-02 

12-1 . 1 634 e+00 

-4.21 74 e+01 

4. PI 90 e+01 

P. 75 75 e- 02 

13 -1 . 9634e+ 00 

5.2381 e+01 

5 . P4 1 7e+01 

. 3.7458e-02 

.14 -1 .9634 e+00 

-5.2381 e+01 

5.2417 e+01 

3 . 7458e-02 


Table E.7 Flexible Beam Closed Loop Poles for G(s) and K mt (s) 

(LQG design based on reduced order model of G(s) 
obtained with mode truncation) 


real imaginary frequency damping 


1 -3.4733e+00 2 .4961 e+OO 4.2772e+00 8.1206e-01 

2 -3. 4733 e+OO -2,4961 e+OO 4.2772e+GO 8.1206e-01 

3-1.0l43e+0l 3.4765e+00 1 .07?2e+0l 9.4598e-0l 

4 -1 .0143'e+01 -r-3.4765e+Od 1„0722e+O1 9.4598e-01 

5 -2. 1 691 e+.00 1.1628e+0! 1.18?9e+0l ! .8338e-01 

6 — 2.1691e+00 — 1.1628e+01 1.1829e+01 1 .8338e-01 

7 -7.3621 e-01 2.2054e+0! 2.2066e+0! 3.3364e-02 

8 -7.3621 e-01 -2.2054.e+01 2 .2066 e+01 3.3364 e -02 

9 -4.0673 e+00 3. 6747 e+01 3.697! e+01 1.1 001 e-01 

10 -4.O673e+00 -3.6747e+01 3.697le+01 l.moie-01 

11 1.51 50e+01 6 . 0890 e+01 6.2746e+01 -2.4 145e-01 

1? 1.5l50e+0l -6.0890e+0l 6.2746e+01 -2.4145e~01 

13 -2.0670e+02 2.0346e+02 2.9003e+02 7.l267e-oi 

14 -2 . 0670e+02 -2.0346e+02 2.9003e+02 7.l267e-0l 
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Results With Modal Truncation 

Another sixth order compensator, K mt (s) was obtained by truncating the 
second order mode with | s \ — o» 3 in the partial fraction expansion of (7(.s) and 
doing a LQG loop shaping design with this sixth order reduced order model. 
This compensator produces an unstable system with closed loop poles shown in 
Table E.7. The differences between the two designs involving K & (s) and K mt (s) 
are readily seen from the loop gain plots for the two designs shown in Fig. E.5. 
While both loops have the same low frequency performance characteristics, the 
peak at s = ju> 3 is less than one for /C 6 (s) and greater than ten for K mt (s). 
Note that the unstable pole has a magnitude of approximately w 3 . 

ACOSS Example 

Sensitivity With Respect To Knowledge of Resonant Frequencies 

The transmission zeros of the compensator K{s ) are given in Table E.8. 
Comparing these zeros with the poles of G 8i (s) it can be seen that the 
compensation involves pole-zero cancellation near the ju axis which is a concern. 
An analysis was performed to assess the stability of the closed loop system in the 
face of simultaneous uncertainty in the critical resonant frequencies given in 
Table E.9 (sensitivity to other resonant frequencies would be less severe). The 
analysis was performed assuming that for any uncertainty the DC gain of C 84 (s) 


was known exactly. 
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1 - 

2 - 

3 - 

4 - 

5 - 

6 - 

7 - 

8 - 
9 - 

10 - 
I I - 
12 - 

13 - 

14 - 

15 - 

16 - 

17 - 

18 - 

19 - 

20 - 
21 - 
22 - 

23 - 

24 - 

25 - 

26 - 

27 - 

28 - 
29 - 


Table E.8 Transmission Zeros of K(s) 


real 

2 ,9816e— 02 
2.981 6e-02 
5 ,5829e-01 
5 . 9824e— 01 
8.5396 e-01 
5 .8649 e-01 
5 .86 4 9 e-01 
I . 2038 e-02 
1 .2038 e-0 2 
4. 42 55 e-02 
4.42559-02 
3. 75 44 e-02 
3 . 75 44e— 02 
9. 0013 e-02 
9.001 3.e-02 
6. 88 5 7 e-02 
6. 8857e-0? 
1 .03 16 e-01 
1 . 03 1 6 e-01 
1 . 0340e— 01 
1 . 0340 e-01 
I .071 2e-01 
1 .07 12 e-01 
I .2568e-01 
1 .2568 e-01 
3.41 54e-01 
3.41 54e-oi 
4..1 1 36e— 04 
4.11 36e— 04 


imaginary 

8 ,8096e— 02 
-8 , 8096e-02 
0. OOOOe+OO 
0. 0000e+00 
0. 0000 e+oo 
7.9556e-01 
-7 . 9556e-01 
-2.091 4e+00 
2.091 4 e+ 00 
-7.7002 e+00 
7. 7 002 e + 00 
-8 . 1 589e+00 
8. 1589 e+00 
1 .0757e+01 
-1 . 0757e+0 1 
1 .4851e+01 
-1 . 4851, e+01 
- 1 .877o ft +oi 
1 . 8770e+01 
-1 . 9954e+01 
1 . 9954 e+OI 
2.1 283e+01 
-2.1 283o+01 
5. 3 099 e+OI 
-5.3099e+01 
-7.1 2 33 e+ 0 1 
7. 12 33 e+oi 
8.5839e+0l 
-8.5839e+01 


frequency 

9.3005 e-02 
9. 3 005 e-02 
5.5829e-oi 
5 . 9824 e-01 
8. 5 3 96 e-01 
9. 3839 e-01 
9.9838 e-01 
2 ..091 5 e+00 
2 .09 15 e+00 
7.7004 e+00 
7. 7004 e+00 
8.1590a+ 00 
8. 1 590e+00 
1 .0/57 e+01 
1 .0757e+0l 
1 .4852o+01 
1 ,4852e+01 
1 .3 770 e+01 
I .3 770 e+01 
1 .9955 e+01 
1 .9955 e+01 
2.1 283e+0l 
2. 1283 e+01 
5. 3099 e+01 
5 . 3099 e+oi 
7. I 234e+01 
7 .1234o+01 
8.5839e+01 
8 .58 39 e+01 


damping 

3.2059e-01 
3 . 2059e-01 
1 . OOOOe+OO 
1 . OOOOe+OO 
1 . OOOOe+OO 
5 . 9339e-01 
5 .9 339e-oi 
5 . 7559e-03 
5 .7559e-o3. 
5 . 7471 e-03' 
5.7471 e-03 
4.6016e-03 
4.601 6e-03 
8. 36 8 Oe- 03 
8. 3680e-03 
4.6363e-03 
4. 63 63 e-03 
5.4960e-03 
5.4960e-03 

5 . 18 16e-03 
5 . 1 31 6e-o3 
5. 03 29 e-03 
5. 03 29 e-03 
2. 3 668 e-03 
2. 3668 e-03 
4.7946e-03 
4 . 7946o-n3 
4. 79 22o-06 
4 . 7922e-06 
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Table E.9 Critical Resonant Frequencies and Their Uncertainty 


Mode 


% Uncertainty 

10 

2.09 

0.5 

14 

7.69 

0.2 

15 

8.17 

0.2 

17 

10.8 

1.0 

21 

14.9 

2.0 

22 

18.8 

2.0 

23 

20.0 

2.0 

24 

21.3 

2.0 

29 

52.5 

20. 

30 

53.9 

20. 

33 

71.3 

20. 

36 

85.3 

20. 

37 

86.2 

i 

20. 






- 257 - 


Recall that G' 84 (s) can be written 


G g 4 ( s ;ft) 


R 


84 C; b; 


+ s 


Ul; 


1=7 w r .+ 4- w," 


where ft denotes the actual values of the thirteen resonant frequencies given in 
Table E.9. Let the nominal system be given by 


Ul; 


R 84 c i 

Gsiis: n) = 4 + £ — — ~ 

5“ 1=7 ujf s~ 4- 2f,-w,'S + U)- 


where ft denotes the nominal values of the thirteen resonant frequencies given in 
Table E.9. Note that in general u> t - ^ Cj ■ for w,-fft and ui ; = u i for uq^ft. 

Consider the additive perturbation implied by the uncertainty in uqe ft 


A(s) — G 84 (s;ft)- G' 84 (s;ft) 


I 8.t. 
ui <n 


C; b • 


Ul: 


s[{oj?-Cj?)s + 2$ i u { Cj i (uj r Cj { )} 
(.s 2 + 2$ { (jJ i s+uJ?){s 2 + 2<; { u> i s+ ui?) 


Note that the resonant frequencies not included in ft do not appear in the 
expression for A(s). Also note that A(0) — 0. 

The stability robustness theorem of Chapter II guarantees stability of the 
closed loop system in the presence of the uncertainty if 


max !cr[A( 7 vj)/v(yw)[/ 4- 
mn ( L L 


Gs A {joj](l)K(juj) 



< 1 


\fuJ 


where ft represents any ft such that 
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U; 

1 - P{ < — < 1 + Pi w,efi 

and the percentage uncertainties, 100 p { are given in Table E.9. This condition 
can be seen to be satisfied by examining Fig. E.6 where the left hand side is 
plotted versus u and can be seen to be less than unity for any frequency. Thus 
the closed loop system is stable in the face of simultaneous percent uncertainty 
given in Table E.9 for the thirteen critical resonant frequencies. 

Closed Loop Poles 

The closed loop poles for the control system with G 2 q{s) anc * K{ s ) are given 
in Table E.10. Note that all have negative real parts. 

Step Responses 

Consider the closed loop system shown in Fig. E.7 where y is a 3 X 1 vector 
containing 



Fig. E.7 ACOSS Closed Loop System 
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Table E.10 ACOSS Closed Loop Poles for G 2 g( 


rea 1 


imaginary frequency 


1 -2.981 6e— 0? 

2 -2. 98.1 6 e-02 

3 -6. 1 360e-01 

4 -6.1 360e-0t 

5 -4.0492e— 01 

6 —4 ..0492 e— 01 

7 -6.4070e-0l 

8 -6.4O70e-01 

9 -9. 3599e-01 

10 -9 . 3599e-01 

11 -7.2532e-01 

12 -7.2532e-01 

13 -1.3664e-02 

14 -1 . 3664e-02 

15 -6. 8568e-03 

16 -6.8568 e— 03 

17 -1 .0646e-02 

18 -1 ..0646 e-02 

19 -I .3(5396+00 

20 -1 . 3639e+00 

21 -I .9 663 e+00 

22 -1 . 9663e+00 

23 -3. 3905 e+00 

24 -3 .3905 e+00 

25 -3.6342 e+00 

26 -3 .63 42 e+00 

27 -2 . 5629e+00 

28 -2. 5629 e+00 

29 -7 .9609 e-01 

30 -7.9609e-0! 

31 -3 . 4522 e+00 

32 -3. 45226+00 

33 -5.2997e— 01 

34 -5. 2 997 e-01 

35 -3.7437e-02 

36 -3.7437e-02 

37 -5 . 32 00 e-02 

38 -5 . 3200 e-02 

39 -5..1 I 40e-02 

40 -5.1140e-02 

41 -2.01 23e-01 

42 -2.ni23 e-01. 

43 -8. 6668 e+00 

44 -8. 6668 e+00 

45 -1 .69QQe+00 

46 -I .6999 e+00 

47 -7 .81 80e— 02 

48 -7.81 80e-02 

49 -1 .2081 e+01 

50 -1 .2081 e+01 

51 — 9. 6233 e— 02 

52 -9 . 62 33 e-02 


8.81 00e-02 
-8 o 8100e-02 
3. 0984e-01 
-3. 0984 e-01 
7.9540e-01 
-7„9540e-01 
6.9728e-01 
-6.9728e-01 
6.1517 e-01 
-6. 1 51 7e-0t 
.1 „7343e+00 
-.1 .7343 e+00 
1 . 9930e+00 
-1 .09306+00 
2.0257e+00 
-2. 02 5 7 e+00 
2. 0926e+00 
-2. 0926 e+00 
1 .9261 e+00 
-1 .9261 e+00 
I .991 6e+00 
-1 . 991 6 e+00 
1 .2301 e+00 
-1 .2301 e+00 
4. 1695e-01 
-4. 1 695e-0 1 
4. 8492 e+ OO 
-4.8492e+00 
6. 709 7 e+00 
-6. 7097e+00 
6. 08 5 4 e+00 
-6.08546+00 
7. 4876e+00 
-7.48 76 e+00 
7 . 6858e+00 
-7. 6 8 58 e+00 
8.1 980e+00 
-8.1 98Oe+0O 
1 . 0824e+01 
-1 .08246+01 
1 . 1 2 77 e+01 
-1 .1 277e+01 
7. 35 10 e+00 
-7.351 Oe+OO 
1 .4005 e+01 
- 1 . 4005e+0 1 
1 .4847e+oi 
-I . 4847e+0 1 
8 . 0556 e+oo 
-8.9556e+00 
1 . 3733e+01 
-1 .87836+01 


9, 3 009 e-02 
9.3009 e-02 
6. 8739 e-01 
6 .8739e-01 
8.92546-01 
8. 9254 e-01 
9.46946-01 
9.4694e— 01 
1 .1201 e+00 
l . I 201 e+00 
1 .8 798 e+00 
1 .87986+00 
1 . 9930e+00 
1 . 09 30 e+OO 
2.0257e+00 
2. 02 5 7 e+ 00 
2. 0927 e+00 
2.0927e+00 
2.3601 e+00 
2.3601 e+00 
2.7988 e+00 
2 .7988 e+ 00 
3. 6 06 7 e+ 00 
3.6067e+00 
3. 6 580 e+OO 
3 . 6580e+00 
5.4848 e+00 
5 . 4848e+00 
6 . 7568 e+ OO 
6. 7568e+00 
6. 996 4 e+00 
*.9964 e+00 
7.5063e+00 
7.50636+00 
7. 6850 e+00 
7.6859e+00 
8.1981 e+00 
8.1 981 e+oo 
1 .0824 e+01 
1 .0824 6+01 
1 .1279 6+01 
1 .12706+01 
1 .1364 e+OI 
I . 1 364e+0I 
1 .4107 6+01 
I .41 07e+01 
1 . 4847 e+01 
| . 4847e+01 
1 .5039 6+01 
I . 5039e+01 
1 .3783e+01 
I .8783 6+01 


) and K(s) 


damping 

3 . 2058e-01 
3.2058e-01 
8.9265e-01 
8.92656-01 
4.5367e-01 
4 .5367e-Ol 
6. 7660e— 01 
6. 7660e-01 
8 . 3567e~01 
8.3567e-01 
3.8584e-01 
3.8584e-01 
6.8558e-03 
6 ,8558e-03 
3.3849e-03 
3.38496-03 
5 . 0874 e- 03 
5.0874 6-03 
5. 7790e-01 
5 . 7790e-01 
7.0258 e-01 
7. 02 58 e-01 
9.4005 e-01 
9. 4005 e-01 
9.93486-01 
9. 9348 e-01 
4 .6 72 7 e-01 
4.6727e-01 
1 .I782e-01 
1 .17826-01 
4.9343 e— oi 
4.93436-01 
7 . 06 03e-02 
7.06 03e-02 
4 . 8709e-03 
4.8709e-O3 
6. 4893e-03 
6 .4893e-03 
4.7248e-03 
4.7248e-93 
1 .7341 e-00 
1 .7841 3 ~02 
7 .6262 p- 01 
7 . 6262e-0l 
1 . 20506-01 
1 .?050e-0| 

5 . 2658 6—03 
5 . 2658e-03 
8 . 03356-01 
8 . 0335e-01 
5.123 4 e -03 
5.1 234 p-03 
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Table E.10 ACOSS Closed Loop Poles for G£g(s) and K(s) (continued) 


53 - 7. 82 38 e— 01 

54 —7.. 8238e-OI 

55 -9. 7238e— 02 

56 -9. 72 38 e-02 

57 -! .0795e-0l 

58 -1 . 0795 e-01 

59 -l.9858e+00 

60 -1 .9858e+00 

61 -3 . 7302 e+00 

62 -3. 7302 e+00 

63 -1 . 64 70 e- 0.1 

64 —I . 6470e-01 

65 -2.5546e— 01 

66 -2.5546e-0! 

67 -2.5671 e-01 

68 -2.5671 e-0! 

69 -2 ,3095e— 01 

70 -2. 3095e-01 

71 -2.7437e-oi 

72 -2.7437e-01 

73 -3.681 4e-0 I 

74 -3 .681 4e-01 

75 -3.5867e-OI 

76 -3. 5867 e-01 

77 -3.61 35e— 0! 

78 -3.6135e-01 

79 -4 . 0078 e— 0.1 

80 -4. 0078 e-01 

81 -4. 4628 e— 01 

82 -4.4628e-01 

83 -6. 7891 e-02 

84 -6.7891 e-02 

85 -4.1536e-01 

86 -4. 1536e-01 

87 -4.4261e-01 

88 -4.4261 e-0 I 

89 -4.92890-01 

90 -4.92890-01 


1 .3999e+01 
-1 .8999 0+01 
1 . 9954e+0f 
-1 .9954 0+01 
2.1 283e+0l 
-2 . 1 283e+01 
2. 1922e+01 
-2.1922,0+01 
2.5085 0+01 
-2.2085e+01 
3.3051 e+01 
-3. 105 1 e+01 
5. 0998 e+01 
-5. 0998e+0! 

5.25280+01 
-5.25280+01 
5.31030+01 
-5.31 03e+0! 

5. 3848 e+01 
-5.38480+01 
7. 1 237e+01 
-7.12370+01 
7. 1 286e+01 
-7.12860+01 
7.op42e+0I 
-7.22420+01 
7. 9957e+01 
-7.99570+01 
8.53420+01 
-8 . 5342e+01 
8.5837e+01 
-8.5837e+0.1 
8.61 720+01 
-8.61 72e+01 
8.3972 0+0.1 
-8.3972e+01 
9.33450+01 
-9 . 3345e+01 


1 .901 5 e+01 
1 .90150+01 
1 . 9955e+0! 
1 .99550+01 
2 . 1 283 e+01 
2.12830+01 
2.201 2 e+01 
2 .20.120+01 
2.2398e+01 
2.2398e+01 
3.30520+01 
3.30520+01 
5.0999 e+ol 
5. 0999 e+01 
5.25280+01 
5 . 2528e+0! 
5. 3103 e+01 
5.31 03 e+01 
5.3849e+01 
5.3849 e+01 
7. 1 2380+01 
7.1238 0+01 
7. 1 287e+01 
7.1287 0+01 
7. 2?43e+01 
7.22430+01 
7. 9953e+01 
7.9958 0+01 
8.53430+01 
8.5343e+01 
8.58370+01 
8 . 5337e+01 
8.61 73 0+01 
8.61 730+01 
3 .89730+01 
8 . 8973 e+01 
9.3346 0+01 
9.8346e+01 


4. 11 46 e-02 
4.11 46e-02 
4.8730e-03 
4.37300-03 
5.0723e-03 
5.0723e-03 
9.021 3 e-02 
9.021 3e-o2 
1 . 6654e-01 
1 . 6654e— 01 
4 .98 300-03 
4.9830e-03 
5 . 0091 e-03 
5.0091 e-03 
4 .83 71 e-03 
4.8371 e-03 
4 . 34900-03 
4.34900-03 
5.0951 e-03 
5.0951 e-03 
5. 1673e-03 
5 . 16 780-03 
5.03130-03 
5.03130-03 
5. no I 9 e-03 
5 . 00190-03 
5.01 24e-03 
5. 0124 e-03 
5.22930-0 3 
5 . 2293e-03 
7 .90920-04 
7.90926-04 
4 .82 01 e-03 
4.8201 e-03 
4.97460-03 
4.9746e-03 
5.01180-03 
5.011 80-03 
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Figure E.6 Stability Robustness Test For ACOSS 
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the Node 11 gyro outputs for attitude about the x, y, z axes measured in radians 
and y c is the commanded value of y. 

Based on the low frequency performance specification of Fig. V.16, the 
design of K(s ) was carried out such that 

J/(«) = C(sI-A+K f CT l K f y c (s) 

where .4, Kf, C are defined in the “Preliminary Design Steps” section of 
Chapter V and are based on rigid body characteristics of the space structure. 
The poles of this rigid body closed loop transfer function are given in Table E. 11. 

It is of interest to compare the step responses of the actual closed loop 
system with A'(s) and G' 2 g(s) to the rigid body response, i.e. that of 
C(sl-A + Kf C)~ x Kj . The x, y, z gyro responses are shown in Figs. E.8a, b for 
an x command, i.e. yj = [1 0 0], The x, y, z gyro responses are shown in Figs. 
E.9a. b for a y command, i.e. yj = [0 1 0], The x, y, z gyro responses are 
shown in Figs. E.lOa, b for a 2 command, i.e. yj =[0 0 1], 

Note that the actual responses are basically the same as the rigid body 
responses. The actual response is a little slower, has a little more overshoot and 
takes a little longer to settle out compared to the rigid body responses. Note also 
that the off axis coupling Figs. E.8b, E.9b, E.lOb is ahvays less than 8^0 of the 


command and less than 2 c c of the command after 6 sec. 
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Table E.ll Poles of the Rigid Body Closed Loop Transfer Function 


real Imaginary frequency damping 

1 -6. 346 5 e- 01 6.3465e-01 8.9753e-01 7.07lte~01 

2 “6. 3 465 e -01 ^6.3465e-01 8.9753e-01 7.07lle-0l 

3 -6 . 6 044 e— 01 6.6044.e-01 9.3400.e-01 7.07l1e-01 

4 — 6. 6044e-01 -6.6044e-01 9.3400e-0J 7.071 1e-01 

5 -9.8923e-01 9.8923e-01 1.3990e+00 7.0711e-01 

6 -9.8923e-01 -9.8923e-01 1.3990e+.00 7.071 le-01 
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APPENDDC F 

TRANSFER FUNCTION FACTORIZATION 
INVOLVING TRANSMISSION ZEROS 

An nth order MIMO transfer function, G(s) with q < n transmission zeros 
at s = Zi, z 2 , ■■■, z q can always by factored in the two ways: 

G(s) = &i(s)Pi(s) = /%(*)£,(«) 

where G,(s) and G g ($) have transmission zeros at s = -z h -z 2 , -z q and 
P,(s) and P 0 (s) are MIMO all pass transfer functions with transmission zeros at 
s — -r j, Zo, ..., z g and poles at s = -z lt -z 2 , -z q . MIMO all pass transfer 
functions have the property that 

PMPft-») = pA-°)p<{°) = i 
P.(=)Pj{-» I = p. t H)p.(*) = / • 

These factorizations are not unique. Consider a unitary transformation, i.e. 
a matrix V such that V T V — I. Let 

<?,'(«) = PM = v ; PM vfVt = i 

PA*) = PA*)''? GA*) = v,ga*) vjv . = i 

then 

G(s) = GMPiW = P.(«)G # («) 

for any unitary matrices V i and V g . Thus the factorizations are only unique up 
to a unitary transformation. 
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These factorizations can be expressed in state space terms [Ref. 31]. Let 
{A,B ,C,D} be a realization of <2(s) i.e. <?($) — C(sI-A)~ l B + D. The 
factorizations are expressed in terms of the generalized eigenvectors for the 
transmission zeros generalized eigenvalue problem. The results for a single 
( q = 1) real zero at s = z and a complex pair ( q — 2) of zeros at s = z and 
s = z* are given in Table F.l. More than one zero or one complex pair can be 
handled by repeated application of the formulas in Table F.l. 

These factorizations are particularly useful for obtaining a minimum phase 
approximation of a non-minimum phase transfer function. The factorizations are 
applied for the non-minimum phase zeros of G(s) resulting in a minimum phase 
approximation G,-(s)V'T or V 0 G 0 (s) for any V,-, V 0 such that — I and 

VjV 0 — I. Important measures of the approximation error are defined by 

Ad.) J [c,(»K r ]"‘ [<?(«) 

A.(«) J [<?(»)- V,c„(.)] [v.O.U)]' 1 = P,(s)V?-l . 

Typically the unitary transformations, V,- and V 0 are chosen such that 
A,(0) = A„(0) = 0. This is accomplished by letting 

V( = P,- r (0) 

V. =p.{ 0) . 



- 272 - 


Table F.1 

Transfer Function Factorization 
Involving Zeros In State Space Terms 


Input or Output 
Factorization 

G(s) 

= P(s)G(s) 


General ized 

n T n 



Eigenvalue 

x T W T 

II 

= 0 

Problem 

L J 




G(s) = G(s)P(s) 



1 

1 

IB 

1 

1 


= 0 


Eigenvector 
Normal ization 


w w = 1 


G(s) 


G(s) = C(sl-A) -1 B + D 


G(s) = C(sI-A) -1 B + D 


Im[z] = 0 Special Case 


B, C 


C - C - 2zwx 


B = B - 2zxw T 


P(s) 


P(s) = I 


2z ww T 


s + z 


Im[z] f 0 Special Case 


Definitions 


a = Re[z] B = w T w y = (} - ■— - ^1 - 


v-1 


B, C 


= C-4ayRe[wx H --y-w*x H 


B = B-4ayRe 


xw H -— x*w H 
z 


'] 


P(s) 


P(s) = I — 5 — — ~ \sRe[ww H - — w w H l+Re zww H -aBw w H l 

- 2 +2as+|z | 2 ( L z J L j 


Note: The expressions for B and y in Reference 31 have typographical 

errors. The above expressions are correct. 
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AP PEND IX G 


ACOSS H MODEL DATA 


The ACOSS II transfer function is given by 


: 6 7 


c i °i 


84 

<?(*) = E -r , 

i-i a‘ + 2f,ui,s + btf 


where is a 3 X 1 vector corresponding to x,y and z torque inputs at node 44, 
q is a 3 X 1 vector corresponding to x, y and z attitude outputs at node 11, u { 
is the frequency and ft is the damping ratio for the ith mode. The data for w,-, b { 
and ft is tabulated in Table G.l. The damping ratios are 0.005 for all modes 
except: 7, 8, 11, 12, 13, 16 which have damping ratios equal to 0.7. 

The three rigid body modes corresponding to translation were eliminated by 
letting 


Br = block row [t,- 7 *] 
C'r A block col [ft] 


*= 1 , 2 , ..., 6 
i = l, 2, ..., 6 


Then 


ft V 


R 84 

<?(«) = -4 + £ — ^ 2 

* .-=7 5 " + 2 ft-UftS + OJf 


where 


R J QB* 


Note that R 1 is the 3 X 3 inertia matrix for the center of mass of the structure 
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Table G.l ACOSS H Data 


Mode Freq Torque Incuts at Mode 44 Attitude Outputs at Mode 11 

* Y z x v 7 . 



rad/sec 

1/Mm 

1/Mm 

1/Mm rad 

rad 

rad 

I 

0. OOe+OO 

1 .38 9-05 

-3.849-04 

-4. 00 9-04 1.38O-05 

-3.349-04 -4.009-04 

2 

o.ooe+oo 

-4 .60e-o4 

-4.50 e -o5 

2.2Qe-04 -4.409-04 

-4.509-05 

2.200-04 

3 

0. OOe+OO 

4.83e-P5 

7.459-05 

-I.459-04 4 .8 3e-/i5 

3.450-05 -1.450-/14 

4 

0. OOe+OO 

1 . 3 I e-04 

-2.2O0-O4 

1.21 e— 07 1.71 e-04 

_7.70o_04 

1 .2 10-03 

5 

o. noe+oo 

- 7. 03 e—o4 

4 , 549-/15 

4.029-05 -7.079-/14 

4 . 54o— 05 

4 .029-/15 

4 

0.00a* 00 

-3.01 e-05 

-7.349-04 

-2.019-04 -4.019-05 

-7.849-04 -7.0! P-04 

7 

0. 14e-01 

2 .8 3e-08 

1 ,54e-03 

-8.889-05 — 1 . 75°— '14 

— 7,8 4 9—04 -| ,74e— 05 

8 

1 . 45p+ 00 

-7.-32 9—03 

-2.77—05 

-5,409-04 -1 . 709-/15 

2.749-04 

1 .159-03 

0 

1 .OOe+OO 

4 , 74e-04 

-1 .04e-I 1 

-5.4] e - II 3 . 58e-04 

1 . 1 40-04 -5.81 e-OR 

in 

2.0Oe+no 

2 .529-04 

1 . 739-00 

5.449-00 2.27e-o4 

2 ,50e-n7 

4.03e-07 

1 1 

2 . 78e+00 

-I .439-03 

-7.OO9-1 0 

-3.329-00 -0.31 9-04 

-3 ,45e-04 

1 .449-07 

12 

3.63e+no 

-o.23e-o4 

-5.4 3e— /1 7 

1.21 e-04 -7.009-04 

-4.079-04 -3 . 08e-04 

1 3 

3 ,45e+on 

2,0*19-04 

-1 .20— 04 

5.179-04 -5.55“- n 4 

-3 .479-04 

1 .05 0-/14 

1 4 

7.4Qe+00 

-3.07e-03 

5 .54e-03 

7.709-05 -3 .48e-04 

1 . 1 8 e-05 

1 , OQe-05 

15 

8.1 7e+oo 

2.14 9-04 

4 , 709-05 

-4,810-07 -0.049-05 

8. 51 e-o5 -2.32 9-04 

I 4 

8 . 47e+00 

-0.859-03 

4.059-00 

-0.709-07 0.4)9-05 

-5.079-04 

2.459-04 

17 

1 . 08e+oi 

-3.Q4e-o7 

-3. 229-07 

1.41 e -04 -4.42e-/14 

1 .179-03 -7.47 e -03 

18 

I . 1 4e+OI 

-?. 709-05 

-4.22 9-1/1 

5 ,32o_oo -7.440-07 

1 .539-/17 

2 . 709-/17 

10 

1 .! 4a+0I 

-1 ,03e-05 

I .509-00 

-2.519-08 7 . 87e-04 

Q.08e-07 -o.54e-07 

20 

1 . 1 Qe+OI 

1 .049-12 

3.409-14 

5.449-14 -8. 27e-15 

0. IO9-14 -3.429-14 

21 

1 , 4 0 e + o 1 

1 .559-07 

-2.7! 9-04 

— 2.48o— 04 2 , 40 p— 04 

-7.70o-03 -7.049-04 

22 

1 ,P9e.+oi 

-1 ,72e-o2 

-7. 74—08 

-4.049-08 1.819-05 

-8. 05 e-04 

3.379-04 

23 

2 , oOe+m 

-5.259-08 

-4.1 4 9-/15 

1,079-02 -2.549-05 

4.370-04 -2.209-04 

24 

2. 1 3e+0! 

5.43e-08 

-1 . 209-02 

—2.819—05 — 0 .ORo— 04 

7.01 o-05 

2.479-05 

25 

3.24e+oi 

1 . 55e-1 2 

-2.14 — 15 

7.049-14 4.07e-14 

4.450—14 -3.74e-l4 

2*i 

3.31 o+O! 

-4 .45e—o7 

1 .479-/10 

-7.459-10 -1.759-04 

-7. 18 0-04 

1 .059-/14 

27 

4.05e+0I 

0.32e-l 4 

5.54e_| 7 

-l.04e-!4 2 . 1 I e- 1 4 

-8.1 Oe-I 8 

1 . 07e-t 8 

23 

5.1oe+oi 

-0.01 e-o7 

3 . 25 e-04 

-7.420-04 1,249-07 

-i? ,?4e-o3 

0. 77e-04 

20 

5. 25e+0I 

-7.829-03 

-3.759-1 0 

2.279-00 -7.059-05 

8.11 0-04 -7,439-04 

30 

5.309+01 

1 . 309-04 

-1 .839-00 

7.24e-o8 -1.78e-07 

Q.1'3e-04 -4 . 00e-04 

31 

5.540+01 

-1 .079-05 

3.4/19-12 

2.7QO-11 3.779-04 

-1 ,400-04 

1 .1 79-46 

32 

5 ,54e+0l 

2. 00 e-05 

! ,27e- ! 1 

5. 75e- 1 1 d.7Re_04 

-7.419-04 

0.309-77 

3 3 

7.1 3e+ni 

1 .459-03 

-5.02e-/i8 

l.o 7 0-/17 -5.550-/14 

1.709-05 -1.509-/15 

34 

7. 22e+0 I 

I .70 9-05 

7.009-07 

-a.05°-04 -7.009-04 

3.500-04 -4.1 59-05 

35 

8. ooe+oi 

-3.94 e _o5 

4 .84e—o7 

-1 .7Q0-/14 1 .020—03 

-1 .749-07 -| .709-07 

3* 

fl.53e+01 

?. 72 o_o7 

1 . 729-/12 

-1 .529-07 2.059-/15 

-1 .049-/14 —4 . QOp— 05 

37 

8 . 42e+0 1 

4 . 1 49-04 

-I . 24e-03 

-4. ROe-04 -Q. 029-05 

-1.1 09-03 -4.809-04 

38 

8 .qoo+ni 

2 .109-05 

1 . 4q 9-/14 

-1 . 1 2 0— OS -1 .040-/17 

1 . 779-07 

2 .270-/14 

30 

0.839+01 

5.43 e-05 

| .429-04 

1 . 27 0— 04 -7 . 449-03 

-.7.540-07 -5.470-04 

40 

1 .01 9+0? 

-2 ,85e-o4 

5 . 779—05 

-2 . 2oe— 04 -2.080-01 

-7.489-04 -4.709-04 

41 

1 .04e+02 

1 .079-05 

-7.479-/14 

-4.109-07 -1.019-07 

1.05 0—07 -7.449-05 

42 

1 .059 + 02 

2. 70e-05 

3. IOe-04 

1 . 77e-04 -7. 1 4a —04 

-4.779-05 -7.740-04 

43 

1 ,08»+02 

3 . 4 59-/14 

-4.309-05 

4.079—04 7,779-/14 

-3 ,50 0-/15 -7.409-04 

44 

I . 1 2e+02 

2 .75 9-05 

-1 . 1 79-04 

4. 77e-07 -1 . 78e-07 

-1 .31 0-05 -1 .01 0-04 

43 

1 . 20e+o2 

3. 05 e-O.4 

2. 11 e-05 

-1 .21 e-05 -2.859-/17 

3.440-07 

| .429-05 

40 

1 .40e+0> 

4. 7O0-04 

4 .279-04 

-4. 21 0-/1 7 -2.74.0-/14 

5.009-04 -1,100-/14 

47 

1 . 5 3e+Q2 

4.00e-04 

-1 . ! 3 e-05 

2. 1 7 e-04 -2. 779-04 

1.01 0-07 -7.470-04 

48 

1 .439+02 

2.289-07 

-7.049-/14 

0. 709-07 2.149-05 

2.87 0—04 -5.739-05 

40 

1 . 449+02 

-4.47p-04 

-4 , 45 e-OR 

—4,449—07 4.049—04 

-3 . 40 0—05 

7.309-04 

50 

1 . 46e+o2 

-7.4i e -l4 

3 .24 — 17 

-5.359-17 -4.129-14 

-0.41 e- 1 7 

1 . 77 A- 1 4 

51 

1 ,40p+/)2 

2.309-05 

7.07“-o4 

7.449-00 -1 .749-03 

-1 .750-04 

3 .2 39-/15 
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Table G.l ACOSS II Data (continued) 


Mode Freq Torque Tnouts et Mode 44 Attitude Outputs et Mode || 

* y 7 x / t 

red/sec 1/Mm 1/Mm 1/M*« red red red 

53 I . 75e+03 -0.32e-08 3.63e-06 -7.70e-06 I . 1 4o-03 -3.I8e-04 -4. 75o-04 

53 1.75e+02 2.729-06 -8.45e-o7 -5.57e-o7 -P.jQo-na 1 . 049-04 - 4 .IO 0-04 

54 I . 7d=+03 1.059-04 -3.87<»-06 7.64»-07 3.70e-04 o.||o-04 -1 ,ooo-04 

55 1.SOe+o2 -1 ,84e-00 5 .OO 9 -OP -!.48e-06 -5. Ode-07 r,03o- 0R Q.4Ro_Q8 

56 1.03e+n2 -3.64e-n7 -6.03e-o7 4 . 6 O 0-03 -?.32e-o5 3 . 370-04 1 . 749-05 

57 2.069+07 ! .50=-05 4.55e-07 -Q.78e-07 5 . 479-03 -o.34e-03 | , eqe-03 

58 2 . 11 e+o? -1 ,66e-06 -7.41e-o7 6.5l9- n S -1.40e-o4 7 . 059-04 -3.26e-04 

5Q 2 . I 7e + 0? - 1 . 00e-05 7.64e-07 -I .50o-06 7.ddo-03 4.57e-o7 - 3 . 340-03 

60 2 . 36e+02 -1 .33e-07 - 2 .I 89 -O 8 d.56 e -0R -4.73e-04 3.07 e -05 - 1 . 679-04 

61 2.40e+02 4.440-07 -1 .Rdo-n7 -7.1 eo -07 5 . 400-03 1 . 01 »-oa 7.72e-o4 

62 3.55e+07 -I .31 e-07 -3.36e-07 l.70o_07 |.76o_04 -|.43e-04 3.00e-Q4 

63 2.5de+np -3.QRe-06 7.45e-o8 6 .00 e — 07 - 4 . 759-03 4 . 819-07 - 1 . 439-07 

64 2 . d I e+03 3 ,Q4o_07 -3.05o-07 R.|7<*-08 I.Ho-06 -| ,57e-o4 -5.18e-o5 

65 3.oi e +o2 -2.40e-08 -I ,63e-08 5 . 349 -nR R.70e-05 -|.85e-04 -I. I Re-04 

66 3.21e+o2 -2 .1 7o_o7 -1 .58«-oa -4. 460-00 1 . 080-04 5 . 400-05 - 1 , 710-05 

67 3.d7e+02 -I ,74e-07 -2.d2e-07 -I.549-0R -1.589-03 -4.47o_04 -Q.llo-04 

68 3 . 079+02 - 6 . 449-07 7.11e-oS -4.oq 0 -oq - 0 . 705-04 3 . 300-04 -5.709-04 

60 4.549+07 d,75o_| q -3.6|o_|fl -5.76o-l7 -4.d|o-20 4.3le-l7 -4,d4e-17 

70 4.55e+02 I.O’e-13 -I,78e-I3 4,6|«-08 -6.47e-0Q -7.569-08 1 . | 4e-08 

71 5 . 7de+o2 — 1 ,07e— o7 -4.60o_od 3.67e-oq 4 . 310-05 - 3 . 340-04 7.5de-o4 

72 5 . 709+03 8. I 3e-00 -1.349-07 6 . 3 R 9 -OR 3.20e-05 -3.44o-04 4 . 530-04 

73 6 . 029 +op 5 . O7o_io _i .33o_iQ -2.370-15 5.2?e-17 7 . 4 3 e -1 7 -4 . OO 0 - 1 8 

74 7.00o+07 -3.47e-l4 —1.17 e — | 3 -3.5 Ro-OR 7.41 o— od 1 .47e-o7 -2.04o_0R 

75 d.47e+07 I ,85e— I 3 3.53e-0R -f .d6e-07 -3.34e-03 -3. 36e-03 -l.5le-03 

76 0 . 510+03 - 4 . 3Q=>_OR i,4flo-od i.44e-os 3.24o-o3 2 ,19e-o3 7 . 710—03 

77 I.OIe+03 —7 ,8do— 03 — | . 4R o— OR | . 74 o-Od -5.57o-04 -5..R5e-05 -3.34e-04 

73 l.opo+oe i.4R e -OR - 3 . 050-08 6 . 4 O 9 -IO -I . 045-04 o.00e-05 -3.35e-05 

/Q I.IQ 0+03 1.15o_oq -4.60O-00 7.04e-OR 1.97o-07 1 . 779-07 2.78e-04 

80 I.Ide+03 5. 77e-00 3.76e-IO -Q. 059-00 -3.710-03 -l.06o-03 -6.0|o-04 

81 1 . 310+03 4.odo — 00 — 1 ,o3o — 08 4 . 339-08 2.7R=-03 3.30e-o7 1 ,65e-o4 

82 l.7| =+03 — 3 , 03 o— Od —4 , 40 o — | 0 I .75 “-OR 4. 579-03 ?.72=-07 l.l8e-07 

87 1.559+03 - 1 . 779-10 _3.37e-00 4,559-00 1.049-07 r. 479-05 l.2le-04 

84 1 . 55e+o? 7 . 779-00 7,77p_n - 5 . 940-10 -i.opo-oR -8.73e-05 -1 .58“-04 
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about the x, y and z axes. 


Finally a minimal realization of 


R 

— with 6 states is 
s~ 


obtained as 

4 = C{s!-A)-'B 


where 

A 
C 

The modes with small residues: 18, 19, 20, 25, 27, 31, 32 and the modes 
with frequencies greater than 100 r/sec: 40, 41, ..., 84 were also eliminated as 
part of the preliminary simplification to form (^(s). 


0 / 3 
0 0 
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16. Abstract 


In a number of applications areas, the control engineer is faced with controlling a 
physical system for which an analytical model can be derived in the form of a very large 
number of coupled, first order, linear, time invariant differential equations. This high 
order analytical model is an input to the controller design process for any design tech- 
nique . 

This thesis develops an approach and a technique for effectively obtaining reduced 
order mathematical models of a given large order model for the purposes of synthesis, analy- 
sis and implementation of control systems. This approach involves the use of an error cri- 
terion which is the H-inf inity norm of a frequency weighted error between the full and 
reduced order models. The weightings are chosen to take into account the purpose for which 
the reduced order model is intended. 

A previously unknown error bound in the H-infinity norm for reduced order models 
obtained from internally balanced realizations was obtained. This motivated further deve- 
lopment of the balancing technique to include the frequency dependent weightings. This 
resulted in the frequency weighted balanced realization and a new model reduction technique. 

Two new approaches to designing reduced order controllers were also developed. The 
first involves reducing the order of a high order controller with an appropriate weighting. 
The second involves linear-quadratic-Gaussian synthesis based on a reduced order model 
obtained with an appropriate weighting. 


Several numerical examples are used to illustrate the theoretical developments of the 
thesis. These examples include aircraft and large space structure problems. The examples 
clearly illustrate the usefulness of this research for practical problems . 
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